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THERMAL NEUTRON ABSORPTION CROSS SECTION OF Xe™! 


H. R. FICKEL? AND R. H. TOMLINSON 


ABSTRACT 


The effective neutron absorption cross section of Xe™ has been measured 
with a mass spectrometer by observing the variation in the Cs'/Cs!*7 fission 
yield ratio obtained at various thermal neutron fluxes. Values of 3.15+0.06 
megabarns and 3.27+0.11 megabarns have been determined for neutron tem- 
peratures of 120°C and 137°C respectively. 


INTRODUCTION 


No fission product requires more consideration in reactor operations than 
Xe'® since it is formed in relatively high yield and possesses an abnormally 
high thermal neutron absorption cross section. It is important to obtain the 
best possible value for its cross section as well as for its absolute yield. Measure- 
ments of the total cross section of Xe! have been made at Oak Ridge National 
Laboratories by Bernstein (1956) using a crystal spectrometer, and by Smith 
(unpublished) using a fast chopper. Differences of about 10% exist in results 
of the two methods, and at present there is no basis for favoring the data of 
either one or the other. 

Westcott (1958) has tabulated both sets of data along with calculated values 
of the effective absorption cross section. The values depend on the choice 
between the possible values of g = 5/8 and g = 3/8 for the statistical weight 
factor and hence the total cross-section measurements leave considerable 
uncertainty in the value of the effective absorption cross section. 

The effective neutron cross section has been measured by Petruska et al. 
(1955) and Ivanov et al. (1957) using mass spectrometric techniques. Since 
Cs" is the daughter of the short-lived Xe" nuclide, it is possible to evaluate 
the neutron absorption cross section of Xe! by comparison of the ratios of 
Cs!*/Cs'*7 from fission product samples obtained at high and low neutron 
fluxes. Values of the neutron temperature were not determined for either 
of the effective neutron absorption cross-section measurements made by 
Petruska et al. (1955) or Ivanov et al. (1957) and hence comparison of them 
either with each other or with the total cross-section values is not possible. 
In the present work the effective cross section of Xe'*® has been measured at 


'1Manuscript received December 30, 1958. 
Contribution from Chemistry Department, McMaster University, Hamilton, Ontario. 


*Holder of Canadian Industries Limited Fellowship. 
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two temperatures and compared with the total cross-section values of both 
Bernstein and Smith. 


EXPERIMENTAL 

Block assemblies containing five equally spaced holes were prepared using 
superpure magnesium. Pins of uranium—aluminum alloy (0.108 cm diameter 
and 0.965 cm in length containing 1 mg U**5) crimped in superpure aluminum 
tubing were located in positions on either side of the central hole. Cobalt 
flux monitors (.005 in. diameter and 1 cm in length weighing approximately 
1.15 mg) sheathed in aluminum tubing were placed in the outside positions. 
Two such assemblies were simultaneously installed at positions S-3-5 and 
S-6-3 in the graphite reflector of the N.R.X. reactor, Chalk River, and 
irradiated for 87.83 hours without shutdowns or serious power fluctuations. 
This time at the flux of these positions was more than that required to ensure 
2% precision in the measured values of the effective neutron cross sections. 
During the irradiation the reactor power was constant within 0.5%, 8-hour 
power checks showing random fluctuations about 40.4 megawatts, not ex- 
ceeding 0.3 megawatt at any time. 

The U**-Al alloy pins from the two samples A and B were dissolved 
separately in vycor crucibles using HCI and HNO;. The aluminum sheaths 
were dissolved along with the pins in order to include the recoil fission frag- 
ments. Each solution was evaporated to dryness and the residue redissolved 
in 10 ml of deionized distilled water. The aluminum was precipitated as 
Al(OH); using NH,OH. The ammonium salts in the supernatant liquid were 
destroyed by repeated evaporations with aqua regia. The resulting sample 
was then transferred to the ion source filament of a 90° sector mass spectro- 
meter as described by Inghram (1953). 

The Co® content of the individual cobalt flux monitors after irradiation 
was compared with standards of known activity by means of a well-type 
ionization chamber filled with argon at 30-atmosphere pressure. The ion 
currents from the ionization chamber were determined with a vibrating-reed 
electrometer. The cobalt standards were calibrated both by coincidence 
and 47 counting techniques. 


RESULTS 
The ratio of the number of atoms of Cs'> found after a neutron irradiation 
of U*** to the total number of atoms of mass 135 produced may be determined 
by solution of the following differential equations: 
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where y'** is the per cent cumulative yield of mass 135 chain, and ¢ is the 
2200 meters/second neutron flux, o and ) are the effective neutron cross section 
and decay constant for the nuclides indicated by the appropriate subscripts 
and superscripts. 

In these equations it is assumed that there is a 5% primary yield of Xe! 
which is the amount estimated from the Glendenin hypothesis of equal charge 
displacement (Glendenin and Steinberg 1955). It is also assumed that the 
remaining 95% of this mass chain is produced directly as I'** since the half-lives 
of the percursors are short compared with the irradiation times for which 
solutions to these equations are required. 

Equation (1) may be derived from these differential equations for the 
conditions that all the I'** and Xe" present at the end of the irradiation have 
decayed to Cs!*> and that there is a negligible depletion of U*® during the 
irradiation. 

ae Sg 0.95(e"""'—1)  0.95(e*'— | 
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The experimental values of Cs'**/Cs'*? divided by the fission yield ratio of 
1.043+0.002 previously obtained by Petruska et al. (1955) give the quantity 
Cs5/ Y"% for each of the four simultaneously irradiated samples. The values 
of 54015 shown in Table I have been calculated from the ratios of Cs'*5/ Y1%6 
using equation (1) with values of 53\!° = 0.1034 hr, 54A° = 0.752 hr-', and 
t = 87.83 hours. 


TABLE I 


Neutron absorption rates in Xe! 





Cs'85/Cs!57 corrected for 








Sample decay 71/2 = 30 years Cals / Ys ssa, ro! 
A—Top 0.4868+0.0021 0.4667 +0 .0022 0.1308+0.0014 
Bottom 0.4793 +0 .0027 0.4595 +0 .0028 0.1355+0.0018 
B—Top 0.7649+0.0057 0.7334+0.0057 0.03801 +0.0011 
Bottom 0.7632+0.0039 0.7317 +0.0040 0.03830 +0 .00075 








The value of the flux ¢ for the positions of each of the cobalt monitors can 
be obtained from its disintegration rate dCo®/dt with the following relation: 


dCo” 1 
at NGOS St 


where ) is the decay constant of Co®(4.18 X 10~° sec), @ the effective absorp- 
tion cross section, and S the self-shielding correction. 

In Table II are shown the values of the disintegration rates and weights of 
the cobalt monitors. The effective absorption cross sections of the cobalt shown 
in column 6 of Table II were calculated (assuming 2200 meters/second absorp- 
tion cross section of 36.5 barns) in the manner described by Westcott (1958). 
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The values of the neutron temperature and the epithermal indices required 
for this evaluation have been determined by Bigham (1958) in connection 
with this present work. The flux depression or self-shielding in the cobalt 
was calculated as described by Petruska et al. (1955). The final column of 
Table II gives the calculated values of the flux using the irradiation time of 
87.83 hours. 


TABLE II 
Cobalt monitors data 


Self- | Weight 





Disintegration Neutron* Epithermal 
serve Co®, rate, temp., G, index, ¢X10%, 
Sample position mg counts/sec X 107 2 barns* r n/cm? sect 








A—Top S-3-5 1.163 0.6904+0.0018 120 36.33 0.018 1.141+0.003 


Bottom 1.151 0.8123+0.0027 1.238+0.005 
B—Top S-6-3. 1.134 0.1585+40.0012 137 36.89 0.0098 0.3265+0.0022 
Bottom 1.163 0.1591+0.0010 0.3286+0.0015 





*Epithermal indices defined by Westcott (1958) and neutron temperatures both measured by Bigham (1958). 
tCorrected by factor of 0.977 for self-shielding (Case et al. 1953). 


In order to obtain the absorption cross section of Xe'*® from s540!**@ given 
in Table I, it was necessary to know the flux in each uranium sample. The 
flux at each of the uranium positions was interpolated from the values at the 
cobalt positions given in Table II assuming a linear flux gradient between 
the monitors in each block assembly. These values corrected for a 0.988 
self-shielding factor in the uranium pins are given in Table III. 


TABLE III 


Effective neutron absorption cross section of Xe!5 





Flux (@) 10%, oxe!5, Average @x~!*5, 





Self-serve | Neutron 
Sample position temp., °C n/cm? sec megabarns megabarns 
A—Top S-3-5 120 1.151+0.0040 3.154+0.034 3.15+0.06 
Bottom 1.198+0.0044 3.138+0.042 
B—Top S-6-3 137 0.3243+0.0022 3.264+0.095 3.27+0.11 
Bottom 0.3244+0.0019 3.278+0.068 





The effective neutron absorption cross sections in Table III were obtained 
by dividing 540° given in column 4, Table I, by these fluxes. 


DISCUSSION 

The precision of the values for the effective absorption cross section of 
Xe® shown in Table III includes only the standard deviation of the obser- 
vations carried out in this work. No estimate of the absolute accuracy of the 
values is possible without including the accuracy of the 36.5-barn value 
used for the 2200 meters/second absorption cross section of the Co** and 
the 9.20-hour half-life used for Xe'*®. To a lesser extent values given in Table 
III] would depend on the accuracy of other quantities such as the resonance 
integral for Co®* and the half-life of I'*°. The precision of the cross sections 
obtained from the samples in assembly B are not as good as those from 
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assembly A. This arises from the relatively smaller amount of capture by the 
Xe! at the lower flux rather than from the precision of the measurements. 

Previous determinations of the effective absorption cross section of Xe!*® 
have been made by essentially the same method as in this work. Ivanov et al. 
(1957) have found a value of 3.2+1.0 megabarns and Petruska et al. (1955) 
3.47 megabarns. Direct comparison of these with each other or with the values 
given in Table III is difficult as neutron temperatures were not measured. 
The precision of the 3.47-megabarn value obtained by Petruska et al. was 
difficult to estimate in view of flux variations during the irradiation, a factor 
which was not inherent in the present work. Although the neutron tempera- 
ture was not determined by Petruska et a/., the irradiation position in the 
N.R.X. reactor for this sample was similar to that of assembly B. If the 
neutron temperature at the time of this irradiation was essentially the same as 
for the irradiation in the present work, the agreement of the values is very 
good. 

The total cross-section measurements of Smith (unpublished) and Bernstein 
(1952) have been used by Westcott (1958) to evaluate the effective neutron 
absorption cross section of Xe at various temperatures. The calculated 
values for temperatures corresponding to this work are given in Table IV. 


TABLE IV 


Effective absorption cross sections of Xe! calculated from total cross sections 


éxe'*§ (g = 3/8) X105, barns oxe'5 (g = 5/8) X10, barns 














120°C 137° C 120° C 137° C 
Smith 29.78 29.96 33.78 33.98 
26.95 27.12 30.19 30.45 


Bernstein 


These values are not strictly comparable with those obtained in this work 
since they represent the effective absorption for a pure Maxwellian spectrum 
of neutrons. The contribution of absorption by epithermal neutrons cannot 
be estimated until the resonance integral for Xe" is known. The values com- 
parable with this work would be slightly less than those given in Table IV, the 
greater reduction applying to the 120° values for which the epithermal flux is 
greater. Comparison of the possible values for the effective cross section of 
Xe! given in Table 1V with the values obtained in this work given in Table 
III suggests, however, that the values calculated from the data of Bernstein 
with the statistical weight factor of g = 5/8 are most consistent with the 
present work. 
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A *=-* TRANSITION IN CrH! 


BENGT KLEMAN? AND ULLA UHLER 


ABSTRACT 


The infrared band system of CrH has been photographed at high dispersion. 
Rotational analyses of the 0-0 and 0-1 bands are given. The transition is shown 
to be a ®2+—8E*, The upper *Z-state shows an extended perturbation. Rotational 
and vibrational constants for the two states are given. The application of theo- 
retical expressions for the fine structure splitting in °Z-states is discussed and 
approximate values of the splitting constants in the observed states have been 
derived. The corresponding system of the isotope molecule CrD was also photo- 
graphed and a vibrational analysis has been carried out. Finally the electronic 
configurations of the molecular states and the nature of the perturbing state are 
discussed. The ultraviolet systems of CrH are briefly mentioned. 


INTRODUCTION 


In the band spectrum of CrH one infrared and two ultraviolet systems are 
known (Kleman and Liljeqvist 1955; Gaydon and Pearse 1937). These systems 
have now been photographed in emission from a King furnace with the 21-ft 
grating at the National Research Council of Canada, Ottawa. As expected the 
infrared system shows a simpler structure than the ultraviolet ones. A rota- 
tional analysis of the 0-0 band of the infrared system was undertaken and it 
is shown in the present paper to be due to a ®2+—*Z* transition. It was not 
possible at the time, due to a temporary lack of suitable plates, to obtain 
useful plates of the 0-1 band. This band was later photographed with a plane 
grating of somewhat smaller dimensions at the Physics Department, Univer- 
sity of Stockholm. The accuracy of the measurements of the lines in the 0-1 
band is slightly less than in the 0-O band. Thus the analysis of the 0-1 band 
does not contribute essentially to the study of the fine structure of the ®°2-terms 
other than to serve as a useful check on the assignment of the lines in the 
perturbed regions of the 0-0 band. 

There have been comparatively few studies of the fine structure splittings 
in high multiplicity 2-terms. The most thoroughly studied ones are probably 
the 7Z-term in MnH (Nevin 1942, 1945) and the ‘4-terms in O.+ (Nevin 
1938, 1940) and GeH (Kleman and Werhagen 1953). 

In a 2-X transition generally only the lines in the main branches can be 
observed. In the CrH bands no satellite branches were observed and conse- 
quently the differences between the sublevels of a rotational state cannot 
be directly measured. Such a transition does not seem specially suited for a 
study of the fine structure. In the present case, however, the splittings are 
rather large. This makes it possible to determine the splitting constants in an 
indirect way. Experimentally one obtains differences (4:F-values) between 
rotational levels within one and the same substate. Comparing the observed 
differences with the corresponding theoretical expressions one obtains approxi- 
mate values of the splitting constants. 

1Manuscript received January 19, 1959. 

Contribution from the Department of Physics, University of Stockholm, Stockholm, Sweden. 

2National Research Council Postdoctorate Fellow, 1954-1956. 
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EXPERIMENTAL 


The CrH bands were photographed in emission from a King furnace. The 
furnace was filled with hydrogen to a pressure of about 1/3 atm. A few grams 
of chromium were placed in a separate graphite tube in order to prevent the 
heating tube from being too quickly destroyed. The temperature of the furnace 
was kept in the range of 2150° to 2250° C. The 0-0 band of the red system was 
photographed at Ottawa in the first order of a 21-ft concave grating at a dis- 
persion of 2.5 A/mm. The plates used were Kodak I-N, I-M, and I-Z and 
the exposure time ranged between 10 minutes and 1 hour. In the third order of 
the same spectrograph the ultraviolet systems of CrH were photographed on 
Kodak III-O plates with an exposure time of about 5 minutes. In order to 
study the vibrational isotopic effect of the red system the furnace was filled 
with a mixture of helium and deuterium to a pressure of 1/3 atm and the region 
of the 0-0, 1-0, and 2-0 bands was photographed on I-N plates with an 
exposure time of 10 minutes. As reference lines for the spectra mentioned 
above high-order lines from an iron arc spectrum were used. 

The 0-1 band was photographed in the first order of a plane grating in 
Littrow mounting at a dispersion of about 4 A/mm. The pressure was 2/3 atm 
and the exposure time 20 minutes on sensitized Kodak I-Z plates. First-order 
lines of the arc spectra of barium, strontium, and chromium served as 
comparison spectrum. 


Structure of the Bands 

As can be seen from a comparison between Figs. la and 10, the structure of 
the infrared system of CrH is simpler than that of the ultraviolet one. A 
thorough study of the infrared 0-0 band showed 12 branches of approximately 
the same intensity. No strong lines remained after picking out these branches. 
The branches form six pairs of associated R- and P-branches. The straight- 
forward assumption, which is also supported by the further analysis, is that 
they constitute a °2—*S transition. As is apparent from the intensity formulae 
for °2—*Y transitions given by Premaswarup (1953) only the 12 main branches 
will be of any appreciable strength. In the present case we have not been able 
to trace any of the satellite branches. All the branches show one big pertur- 
bation and some branches, in addition, show smaller ones. The perturbations 
originate in the upper state. 

The exposures of the 0-1 band are not as good as those of the 0-0 band. 
Here, too, six R- and six P-branches are found. The agreement of the com- 
bination differences and the presence of the same perturbation pattern confirm 
the analysis of the 0-0 band. 

The CrD plates show a rather dense structure. The systems of both CrH 
and CrD are present and we were only able to measure the CrD heads of the 
0-0, 1-0, and 2-0 bands. 

The ultraviolet bands have, as is seen from Fig. la, a very complicated 
structure. No attempt has been made to analyze these bands rotationally. 


Vibrational Analysis of CrH and CrD 
As only two bands have been analyzed rotationally, the vibrational constants 
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band heads 





Fic. 1. CrH bands. (a) Upper band. The 3600-A system. (b) The 0-0 band of the infrared system. Part 
of the analysis is shown in the figure. 
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had to be derived from measurements of the band heads. The heads measured 
are listed in Table I. 











TABLE I 
Band heads 
CrH CrD 
Band 1, A Int. cm! Band , A Int. cm 
0-1 9968.5 3 10029 
0-0 8611.1 10 11610 0-0 8626.6 10 11589 
1-0 7641.7 4 13082 1-0 7901.8 3 12652 
2-0 6889 .6 1 14511 2-0 7302.9 1 13689 





Using the simple isotope formulae the w, value of CrD can be calculated 
from the vibrational constants of the hydride. The value obtained is w/ = 1083 
cm~! in good agreement with the observed value 1089 cm. In the lower 
electronic state it is only possible to obtain the interval between the first two 
vibrational levels. 

CrH: AG,’ = 1581 cm“. 


The vibrational constants are summarized in Table V. 


Rotational Analysis 
The 0-0 and 0-1 bands were analyzed. The wave numbers of the lines of 


the 12 main branches in each band are given in Tables II and III. 

As mentioned before, the measurements of the lines in the 0-1 band are less 
accurate than those of the 0-0 band and only the combination differences of 
the 0-0 band were used for the determination of the rotational constants of the 
upper state. 

A general procedure for obtaining the rotational energy levels of a =-state 
has been developed by van Vleck (1929). It leads in the case of a *Z-state to a 


secular determinant of the sixth order. The energy levels W,... Wg. are the 
roots of the following equation: 

a-W aq 0 0 0 0 

ay b—W by 0 0 0 

0 b, c-W ¢g 0 0 

0 0 Cy c—-W b, 0 zs Q 

0 0 0 b, b- W ay 

0 0 0 0 a a—W 
where 
a = BIJ(J+1) —(15/4)]+10e; 6 = BIJ(J+1)+(17/4)]—2e; 

c = BIJ(J+1)+(33/4)]—86¢ 

a, = V5BIJ(J +1) —(15/4)}; b) = 2V 2BlU(J +1) — (3/4); cp = 3B +3). 


The magnetic interaction between V and S can to a first approximation be 
taken into account by adding the expression (y/2)[J(J+1)—N(V+1) 
—§(S+1)] where y is a constant. 
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The splitting constants do not occur in the expression 
6 
> AF,(J). 
y=] 


Thus by plotting 
€ 
be airs) |/ (J+4) 


against (J+}3)* the rotational constants B and D can be derived. The values 
obtained are given in Table V. 


Determination of the Splitting Constants 

The secular determinant of the sixth order mentioned above reduces to two 
of the third order. A first-order development of the roots of these equations 
leads to the following expressions for the energy levels (Rao 1952). 


F\(N) = B | ww+1)-10 ew where N = J-—5/2 
F,(N) = nl (N+1)+A (x9) N = J-8/2 
FAN) = bi ail stay r3— =) N = J+1/2 
F,\(N) = Bl v4u+a(1 ty N = J+3/2 
F,(N) = aly N(N+1)—100 ae N = J+5/2 


The parameter \ is equal to the spin splitting constant e divided by B 
(A = «/B). 

In order to study the influence of the neglected high-order terms of \ the 
equations have been solved numerically for the four values \ = +1.0, +0.1. 
The result, compared with the first-order solution is seen in Fig. 2. The dis- 
crepancies between the two solutions are small for \ = +0.1, but quite im- 
portant for \ = +1.0. It is to be observed that the exact solutions are not, in 
contrast to the approximate ones, symmetric with respect to a change of sign 
of \. A rough calculation showed that the values of \ to be expected in the 
present case are of the magnitude +0.1 and +0.01 in the upper and lower 
state respectively. Thus the approximate energy expressions of the rotational 
sublevels can be assumed to be sufficiently accurate for the determination of the 
splitting constants in the present case. 

In the absence of satellite branches it is impossible to obtain experimental 
values of differences between any two of the sublevels in a rotational state. 
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The only observational data convenient for the determination of the splitting 
constants are the second differences: 

















AFi(N) = Bl aw+)— aya | +5y 
AsF:(N) = | s+ +4) ~ayp wen | +3y 
aA) =e [4v+4)- sere 5] ¥Y 
wri = s[sren arctan] 
Set) = Blt + aye es | =a 
AaFa(N) = B[ w+ +a a | ~5y 


In the expressions above, the term due to JN, S interaction 
E YN, +3 (8N—5), +5 (V—-8), —5 (N+9), —5 (@N+8), and 


oY 
“a w+1)| 


in F\(N), Fo(.N), F3(N), PCV), Ps(V), and Fe(V) respectively has been added. 
The constants \ and y are derived graphically by a linear representation of 
expressions of the type 


AF m(N)—A2F,(N) (where m= 1,2...6,”=1,2...6and m#n). 


Every combination of m and m has been examined. The scattering of the 
individual determinations is considerable, however, and the accuracy in the 
y’’-, \”’-, and \’-values obtained is roughly 30%. The magnitude of - is much 
larger than that of y’’ which reduces the siadsi error to about 5%. For the 
upper state some of the curves were slightly bent indicating that the approxi- 
mate expressions are not quite sufficient. Mean values of the constants derived 
from the diagrams are given in Table V. 


The Perturbations 

As mentioned before, large irregularities in all main branches have been 
found. 

The displacements of the lines from their expected positions, which in some 
cases amount to 7 cm, are explained as perturbations in all six sublevels 
of v’ = 0. This is checked by the analysis of the 0-1 band which shows exactly 
the same perturbation picture as the 0-0 band. The perturbation maximum 
occurs at N = 20, 21 in the levels F.-F, and at N = 19, 20 in the F;-level. 
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In addition to the main perturbations smaller ones have been found at N = 9, 
10 in the F;, Fs, and F; levels. The perturbations are illustrated in Fig. 3. 





10.5 20.5 t 


Fic. 3. The perturbations in the level v’ = 0. The unperturbed parts of the sublevel energy 
functions have been schematically illustrated as straight lines. The magnitude of the pertur- 
bations has been represented on a scale given on the ordinate axis. The separations of the lines 
in the figure are not significant. 


A ®d-state necessarily belongs to Hund’s case 6. The occurrence of the 
perturbations at the same N-number in all the sublevels indicates that the 
perturbing state also belongs to Hund’s case b and that its multiplicity must 
be at least six. For reasons which are given below the perturbing term is 
assumed to be an °2-state. 

Budé and Kovacs have treated the case of perturbations between 2-terms 
of equal multiplicity (Kov4cs 1938) and perturbations between 2-terms of 
different multiplicity (Bud6é and Kovacs 1938). A necessary condition among 
others for a perturbation is that the J-values are equal in the perturbed and the 
perturbing states. In the case of multiplet terms the N-values must be equal 
or differ by an even number. 

In the present case it is natural to assume that the large perturbation at 
N = 21 occurs when both J and N are equal for the perturbed and the per- 
turbing state. The smaller perturbations at V = 10, which have only been 
partially traced, may represent the case, when the V-numbers differ by 2 units. 
From these assumptions and assuming the rotational energy to be roughly 
BN(N+1), we are now able to calculate an approximate B*-value of the 
unknown perturbing state. The value obtained is B* ~ 6.0 cm™. 








qr 
~ 
N“N 


KLEMAN AND UHLER: *-8 TRANSITION 


Expected Electronic States 

In general the low-lying electronic states in diatomic hydrides can be pre- 
dicted on the basis of the electron configurations in the ‘‘united atom’’. When, 
however, there are several electron configurations in the ‘‘united atom” 
giving rise to electronic states of approximately equal energy, there will be 
difficulties in predicting the order of the states in the corresponding hydride 
molecule. This is the case for CrH as well as for the neighboring molecules as 
for example MnH. 

It therefore does not seem feasible to predict the expected states simply by 
trying to consider the relative order of excitation of all possible electron con- 
figurations. Instead we proceed in the following manner. First, accepting the 
experimental evidence of the two low-lying ®Z-states, we will try to relate 
these to the proper electron configurations, and then discuss the reason for the 
assumption that the lowest of these observed states is the ground state of the 
molecule (X ®2), and finally we will compare the result with the related mole- 
cule MnH. We only need now to consider those electron configurations which 
can give rise to sextet and higher states and shall assume that of all molecular 
states arising from one electron configuration the one with the highest spin 
will be lowest and that, when there are several with the same spin, the one 
with the highest A-value will lie lowest. In Table IV those electron configura- 
tions which can give rise to low-lying sextet and octet states are shown. 


TABLE IV 





Lowest 
KL 3de 3dx 3d6 4s0 Apo Apr state* Other states 
CrH KL @ r? 62 o 6>+... Quartets, doublets 
KL_ a? x & o 7 Quartets, doublets 
KL a? n 5 7 (S11)... Quartets, doublets 
KL a x & ¢ o 8>+... Sextets, quartets, 
doublets 
MnH KL a am oe o o 7+... Quintets, triplets, sing- 
lets 
KL @ ad 6? o 7 TR... Quintets, triplets, sing- 
lets 











*This is the lowest predicted state of the given electron configuration and is (with the exception of the (*I)-state) 
assumed to be identical with an observed state. 


The two observed *2-states in CrH are interpreted as being the *2*-states 
derived from the electron configurations (3da)*(3dx)?(3d6)?(4sc) and (3da)? 
(3dm)?(3d5)?(4pc) respectively. Both these states are the only sextet states 
which can be derived from these electron configurations. Although it cannot 
be rigorously proved from the present observations, there are many indica- 
tions in favor of the assumption that the observed X *+-state is the ground 
state of the CrH molecule. From the experimental point of view the strength 
of the A *S+—X ®E+ system, considering the observed behavior of other 
hydride spectra in a low-energy light source of this type supports the above 
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assumption. The ground state of Mn is a 8S-state and this has to be correlated 
with the ground state 7S+2S of Cr+H via the lowest *2+-state in the CrH 
molecule. By the assumption that this *2+ is the ground state we avoid a two- 
fold-crossing of this potential curve by another presumed ground state curve 
in the molecule. The above assumptions also look very reasonable in view of 
the electron configurations which one is forced to assume in MnH to explain the 
two lowest observed states 72+ and "II, viz. (3dc)?(3dm)?(3d5)?(4sc)(4po0) and 
(3do)?(3dx)?(3d5)?(4s0)(4pr) (Herzberg 1957). From the appearance of the 
perturbation in the A *2+-state it is evident that the perturbing state, which 
has to be a sextet or an octet, is very near to a case b state. Except for 2-states 
pure case b states are very rare in molecules containing a relatively heavy 
atom like Cr and, therefore, the perturbing state is in all likelihood a =-state. 
If this is true, it has to be the *2+-state indicated in Fig. 4. The perturbation is 
rather weak which also fits with the relatively weak interaction between 





Fic. 4. Potential curves of CrH. 


states of different multiplicities. The *2+-state can only be derived by breaking 
up the (3dc)? electron pair, and it seems likely that it corresponds to the 
electron configuration (3da)(3d7)?(3d6)?(4sc) (4pc). 

The lowest octet term in the united atom Mn is an *P° term, which via the 
perturbing ®2*+ is to be correlated with 7S(Cr) + 2S(H). The upper A *2*-state 
of the infrared system may be correlated with the 5S(Cr) + 2S(H) level of the 
separated atoms. The suggested term scheme is given in Fig. 4. 

An extrapolation of the vibrational levels of the upper *2+-state (A ®2) in 
CrH and CrD points to a dissociation energy D for this state of around 23000 
cm~'. This gives for the dissociation energy of the ground state (X ®2+) D?’ 
= 26000 cm~. The shapes of the potential curves of A *£+ and X ®+ (Fig. 4) 
have been based on these figures. 
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The Ultraviolet Systems 

The complex structure of the ultraviolet bands rules out the possibility that 
they are of the 2-2 type. The strength of the 3600-A system probably implies 
a transition which is in accordance with the selection rules for an allowed 
transition. Assuming the ground state to be involved in this transition, the 
term scheme in Fig. 4 would suggest the upper state to be a ®II-state probably 
correlated with the ®D(Cr) + 2S(H) dissociation limit. 


TABLE V 
Summary of constants of the A *2-X ®F system 


By’ = 6.125 cm™ By = 5.268 cm™ 
D,’ = 3.38X10-4 Dj = 2.42107! 
CrH: Bi’ = 5.945 
1 = 3.4010 voo = 11542.88 
a” = 0.180 w, = 1515 
AG = 1581 wixs = 21s 
/’ = 0.15+0.01 y’ = 1.20+0.05 
A” = 0.03+0.01 ’ = 0.12+0.04 
é’ = 0.18+0.06 e’ = 0.66+0.20 
Note: For CrD the constants obtained are: 
We = 1089 Wexe = 13 


The vibrational constants are derived from band heads. 


ACKNOWLEDGMENTS 


Our thanks are due to Drs. A. E. Douglas and G. Herzberg for many 
helpful suggestions during the course of this work. We are also indebted to 
Dr. Lars Akerlind for experimental help with the exposures made in Stockholm. 


REFERENCES 


Bupé, A. and KovAcs, I. 1938. Z. Physik, 111, 633. 

Gaypon, A. G. and PEARSE, R. W. B. 1937. Nature, 140, 110. 

HERZBERG, G. 1957. Molecular spectra and molecular structure (D. Van Nostrand, Com- 
pany, Inc.), p. 341. 

KLEMAN, B. and WERHAGEN, E. 1953. Arkiv Fysik, 6, 399. 

KLEMAN, B. and Litjegvist, B. 1955. Arkiv Fysik, 9, 345. 

KovAcs, I. 1938. Z. Physik, 111, 640. 

Nevin, T. E. 1938. Phil. Trans. Roy. Soc. London, A, 237, 471. 

1940. Proc. Roy. Soc. (London), A, 174, 371. 

1942. Proc. Roy. Irish Acad. 48, 1. 

—— 1945. Proc. Roy. Irish Acad. 50, 123. 

PREMASWARUP, D. 1953. Indian J. Phys. 27, 415. 

Rao, K.S. 1952. Indian J. Phys. 26, 47. 

VAN VLECK, J. H. 1929. Phys. Rev. 33, 467. 




















ABSOLUTE STANDARDIZATION OF RADIOACTIVE NEUTRON 
SOURCES 


II. THE USE OF THE F'(a,n)Na® REACTION! 


kK. W. GEIGER 


ABSTRACT 

Fluorine has only one stable isotope, F'*. If neutrons are produced by the 
F!9(q@,n) Na® reaction the neutron output can be calculated from the yield of the 
resulting radioactive Na*. The growth of Na® (half-life, 2.58 years) has been 
measured in a neutron source consisting originally of 1.6 curies Po#® mixed with 
CaF: powder. Since Na® is a positron emitter, discrimination against y-rays 
from Po?!° and from nuclear reactions could be achieved by detecting the two 
positron annihilation quanta in coincidence. The Na” growth has been followed 
over 20 months and is in agreement with the theoretical growth curve. Com- 
parison with a calibrated Na* source yielded a neutron emission rate of (10.70 
+0.25) X 104 sec™!. This resulted in a neutron emission rate of (3.16+0.10) X10 
sec~! for the Ra-a-Be source of the National Research Council, in good agreement 
with (3.22+0.05) X 10° sec™ obtained by a neutron thermalization method. 


INTRODUCTION 


In a previous publication (Part I, Geiger and Whyte 1959) the principal 
methods for the determination of the neutron emission rate, Q, from a radio- 
active neutron source have been discussed and an experimental determination 
of the output of the N.R.C. standard neutron source has been made. This was 
carried out by using a method where substantially all neutrons were ther- 
malized and captured in an aqueous MnSOQ, bath. Since the proportion of 
captures in manganese can be calculated the neutron emission rate of the 
source can be determined by the absolute activity of the produced radioactive 
manganese-56. 

An alternative for the determination of a neutron yield is to use a reaction 
in the source where the (recoil) particle which is produced simultaneously 
with the neutron can be detected (associated particle method). This has the 
advantage that no knowledge of neutron cross sections is required. Accelerator- 
induced reactions lend themselves easily to the associated particle method; 
however, the comparison of their neutron yield with that of a radioactive 
source is difficult. It is therefore more desirable to use an associated particle 
method within a radioactive neutron source. This principle seems to have been 
used only on photoneutron sources. Gliickauf and Paneth (1938) used the 
Be*(y,m)2He* reaction and analyzed for the amount of helium produced in 
the beryllium. Martin, Bishop, and Halban (1955) counted the photoprotons 
of the H?(y,7)H! reaction in an ionization chamber. The deuterium gas was 
then replaced by heavy water to obtain a reasonable yield. Such standard- 
izations are very difficult to carry out experimentally. 


1Manuscript received January 26, 1959. 
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However, the associated particle method could also be used in the more 
convenient (a,) source when the product nucleus is radioactive. This, as 
far as the author knows, had not been done previously. The following require- 
ments would have to be met by such a source: 

(1) The neutron-producing reaction should give a high yield. 

(2) All neutrons should originate from the same reaction. 

(3) Other possibly concurrent reactions (e.g. a,p) should not produce 
radioactive nuclides. 

(4) Measurement of the growth of the radioactive product should be possible 
without destruction of the source. 

(5) The source should be of small size. 

The 

Bt, E.C. 


2 


> Ne” 


(1) F!9(a,n) Na” (—1.93 Mev); Na” 
2.58 years 


reaction is probably the only one which fulfills the above requirements well. 
The neutron yield from fluorine is only excelled by beryllium and boron (e.g. 
Feld 1953) and fluorine has only one isotope; hence for every neutron formed 
a Na”™ nuclide is produced. To meet requirements (4) and (5) polonium-210 
may be chosen as a-emitter because of the small number of y-quanta emitted 
and its high specific activity. Such a source has been prepared and the neutron 
emission determined absolutely. 


THE Po-a-F NEUTRON SOURCE 
The source consisted of a nickel cylinder, 8 by 18 mm high and filled with 
1 g of CaF. powder; 1.6 curies of Po” were distributed by heat through 
the fluorine.’ Calcium has higher reaction thresholds than fluorine and has not 
been found to react with polonium a-particles (Odeblad 1954; Saha 1938). 
Besides reaction 1 and inelastic scattering of a-particles only the following 
additional reaction is expected to occur 


(2) F'9(a,p)Ne” (1.70 Mev). 


To be useful as a standard the source has to be physically stable with time. 
This has been checked by comparing its neutron emission rate with that of a 
Ra-a-Be source over a period of 20 months, using a BF; proportional counter 
in paraffin. As expected, the neutron emission rate of the Po—F source followed, 
within the uncertainty of the measurements (+0.2%), the well-established 
Po?!® half-life of 138.4 days. 

The y-spectrum of the source has been studied with a scintillation spectro- 
meter, using a 141 in. sodium iodide crystal. The spectrum (Fig. 1) of the 
30-day-old source is in general agreement with the spectrum obtained by 
Breen and Hertz (1955) and shows mainly the 0.80-Mev y-peak of Po*!° and a 
further photopeak at 1.28 Mev. This peak is due to de-excitation of the 


2The source was manufactured by Mound Laboratory, Monsanto Chemical Company, 
Miamisburg, Ohio. 
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first-excited state of Ne”, reached by reaction 2. Although the identical state 
is reached by the 2.58-year decay of Na”, the y-ray cannot be assigned to this 
decay, since after only 30 days very little Na” has been produced. In addition, 
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Fic. 1. y-Spectra from Po-a-F source at various times after preparation. The intensity 
relations between the spectra are arbitrary. Also shown is the spectrum of pure Na”. 


the expected 0.51-Mev photopeak from positron annihilation is not apparent; 
the peak just above 0.51 Mev corresponds to the Compton edge of the 0.80- 
Mev y-ray. The peak at 1.55 Mev is probably identical with a y-ray of the 
same energy observed by Temmer and Heydenburg (1958), who assigned it to 
a transition in Na*. The photopeaks at 0.21 Mev and part of the 0.085-Mev 
peak can be attributed to inelastic scattering of a-particles on F!* which 
produces gamma transitions of 0.197 and 0.109 Mev (Temmer and Heyden- 
burg 1958). Part of the 0.085-Mev peak is from polonium X-rays and from 
fluorescence radiation of a lead shield surrounding the detector. 

With increasing age of the source the growth of Na” becomes noticeable 
by the appearance of the 0.51-Mev photopeak (Fig. 1). For comparison the 
scintillation spectrum of pure Na” is also given on this figure. 


DETERMINATION OF Na® CONTENT IN THE Po-F SOURCE 
The number of Na” nuclei present at the time ¢ after source preparation is: 


aes 2 Vi, 0 


(3) N2= i, [exp(— Ast) —exp(—A,t)] 
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where 


\ is the decay constant, the subscripts 1 and 2 referring to Po®#!° and Na”, 

Ni,0 is the number of Po*° nuclei present at ¢ = 0 (time of source pre- 
paration), 

k is the number of Na* nuclei produced per Po?! disintegration. This is, 
according to reaction 1, equal to the number of neutrons produced. 


The neutron emission rate at the time of source preparation therefore is 
Qo = RriN 1,0. 


It is more convenient to express the Na” content as a disintegration rate. 
The equation for the growth of the Na” activity then becomes: 


(4) AWe = Sen [exp(— Ast) —exp(—Ajt)]. 

Ai—As 
Hence the initial neutron yield of the source can be calculated once the 
absolute decay rate of the Na*® is determined at a time ¢ after source pre- 
paration. It is of course preferable to establish a substantial part of the growth 
curve experimentally. 

Since it is virtually impossible to obtain the Na* content from the scintilla- 
tion spectra of Fig. 1 to the required accuracy, some other method had to 
be chosen to discriminate against the gamma radiation not belonging to Na”. 
This nuclide is essentially a positron emitter, the positrons being absorbed 
within the source. By a coincidence arrangement, use has been made of the 
detection of the two 0.51-Mev annihilation quanta which are emitted simul- 
taneously in opposite directions. Only a negligible fraction of positrons originate 
from another cause: pair production of neutron capture y-rays within the source 
and its surroundings. 

Two scintillation counters with 1} X1 in. Nal crystals were used and com- 
parisons carried out to an originally uncalibrated Na® reference source. The 
scintillation counter bias was set such as to accept y-radiation above either 
250 kev or 380 kev. The results agreed within the statistical accuracy which 
is to be expected since the effect of backscattering is negligible even at 250 kev. 
The counters were placed at a 180° angle for the measurement of the annihi- 
lation radiation and at 90° for background and accidentals. A source-to- 
crystal distance of 13 cm resulted in a convenient counting rate. The positron 
emission from the Po—F source as a function of time is shown in Fig. 2. The 
experimental points follow well the calculated growth curve of Na” as given 
by equation 4. This again strengthens the assumption that the source was 
physically stable and that the contribution of positrons from other causes 
than reaction 1 must be negligible. 

The coincidence method is very sensitive to source position and to y-absorp- 
tion within the source. Therefore, for the absolute Na” determination a capsule, 
identical with that of a source, was made. This was filled with CaF, powder 
and an aliquot of Na® standard solution, dried in the capsule. The Na” 
solution was standardized by the 47 6-counter method as described for Mn*® 
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in Part I. In addition, the self-absorption of 6*+-particles for one of the counting 
samples was determined by the 4% B—y coincidence method (Campion 1958); 
(97.540.5)% of the emitted §+-particles were counted. For the calculation 







Growth of Na**-activity 
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° 100 200 300 400 500 600 700 800 900 1000 
Days 
Fic. 2. Positron emission of Po-a-F source as a function of time, based on comparison 
with a Na” reference source. The growth curve (equation 4) is fitted to the experimental 


points. Also shown for comparison is the decay of Po?! which is identical with the decrease 
in neutron emission from the source. 


of the Na” disintegration rate a branching ratio for electron capture of (10.2 
+0.5)% has been assumed. This is the theoretical value for an allowed 
transition with the error estimated from the more recent experimental de- 
terminations (e.g. Hagedoorn and Konijn 1957, who also give previous 
references). 

The Na” disintegration rate in the neutron source was (1.108+0.020) x 104 
sec~!, 560 days after source preparation. With equation 4 and a Na” half-life 
of (2.58+0.03) years (Merritt, Campion, and Hawkings 1957) the resulting 
emission at the time of source preparation (25 January, 1957) is 


Qo = (10.70+0.25) X 104 neutrons per second. 


COMPARISON WITH THE N.R.C. Ra-a-Be NEUTRON STANDARD AND 
CONCLUSIONS 
Because of its short half-life a Po-a-F source itself is only of limited value 
as a neutron standard. Therefore, the Ra-a-Be standard of the National 
Research Council has been calibrated in terms of the Po—F source. The neutron 
spectra of these two sources are very different; the maximum energy is 12 Mev 
for Ra-a-Be compared with 3 Mev for Po-a-F (Szilvasi, Geiger, and Dixon 
1959). The measurements were carried out by thermalization and capture of 
nearly all source neutrons in the MnSQ, as described in Part I. In this case no 








GEIGER: RADIOACTIVE NEUTRON SOURCES, II 555 


absolute activity determination is necessary. After removal of the source and 
thorough mixing of the bath the Mn** activity was measured with a scintil- 
lation dip counter (13 X1 in. Nal Crystal). Some difficulties were experienced 
with this counter in that the stability was rather poor. This was mainly 
attributed to temperature variations. The rate due to Mn*® had to be fre- 
quently normalized to the rate from a pitchblende reference source. At the 
particular discriminator bias used (0.25 Mev) it was found that a sensitivity 
change of 1% on the pitchblende source corresponded to a change of 1.5% 
in the MnSO, bath and the normalization had to be carried out accordingly. 

During the first 6 months after source preparation three comparisons, 
each consisting of five bath activations, were carried out and the strength 
of the N.R.C. neutron standard was found to be (28.1+0.4) times that of the 
Po-F source for the day of its preparation. Corrections have to be applied 
for fast neutron capture in the bath (2.8+0.4)% and escape from the tank 
(1.4+0.4)% as discussed in Part I. The losses only apply to the Ra-a-Be 
neutrons since there is virtually no fast neutron loss in the bath from the 
Po-a-F source. Absorption of thermal neutrons in the Ra-a-Be source amounts 
to (0.4+0.2)%. For the Po-a-F source, with a total thermal absorption cross 
section of only 0.15 cm?, this is estimated to be less than 0.1%. Hence the 
corrected ratio becomes 


[(Qo(Ra Be)]/[Qo(Po F)] = 29.4+0.5. 


With the neutron emission rate for the Po-a-F source as found in the 
previous section the rate from the N.R.C. standard becomes 29.4 X 10.70 X 104 
= 3.15X10® sec“! as of January 1957, or when corrected for growth of Po*!® 
(de Troyer and Tavernier 1954): 


Q = (3.16+0.10) X 10° neutrons per second as of July 1958. 


This may be compared with the rate as determined by the manganese bath 
method (Part 1) of (3.22+0.05) sec—! as of July 1958. The agreement between 
the two entirely different standardization approaches is satisfactory and makes 
the occurrence of a major systematic error unlikely. Although the accuracy of 
the present method is not quite as high, its main advantage lies in the fact 
that no cross-section data are required, in contrast to the majority of the 
standardization work carried out until the present. 
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SIZE OF IRREGULARITIES IN THE E REGION 
OF THE IONOSPHERE! 


M. SrrrRAMA Rao? 


‘ ABSTRACT 


Three-station fading records taken at Ottawa on frequencies ranging between 
2.0 and 2.6 Mc/s have been analyzed to determine the E region wind velocities 
(V) and frequencies of fading (N) for the period June 1953 to April 1955. An 
empirical relation V = 2.41 NA has been deduced, on a statistical basis, to 
represent the above observations. In this equation, \ represents the wavelength 
of operation. The constant of proportionality (K) in this relation has been 
studied comparatively for three different locations in the northern hemisphere. 
This constant has been related to a characteristic length ZL which specifies the 
horizontal size of the irregularities in the E region of the ionosphere by the 
equation: L = (1/2)KX. Variation of this size with latitude is discussed. 


1. INTRODUCTION 


Systematic measurements have been made, in recent years, on the deter- 
mination of horizontal drifts in the E region of the ionosphere, using the 
spaced receiver method developed by Mitra (1949). The basic assumption 
in this method is that the lower portion of the £ region acts as an irregular 
reflector and forms a diffraction pattern on the ground when illuminated by a 
point source of radio waves on the ground. It can be easily seen that the 
irregularities of the diffraction pattern on the ground would be twice as large 
as those in the diffracting screen itself. As the ionospheric irregularities 
drift horizontally, with a velocity V, the diffraction pattern on the ground 
drifts past the receiver with a velocity 2V and the received signal exhibits 
random fading. It has been shown by Rao and Rao (1958), from a statistical 
study of three-station fading records at Waltair, that for any operating fre- 
quency, the average frequency of fading (V) is proportional to the ionospheric 
drift velocity (V). This is only statistically true, as there may be a wide 
scatter of individual observations about the mean straight line of the N vs. V 
graph. This means that the size of the ionospheric irregularities at Waltair 
varies with time around a statistical average value, as will be shown later in 
this paper. A plot of the individual observations of N and V for Kjeller, given 
by Harang and Pedersen (1956), shows such a wide scatter. Their graph also 
indicates an average increase of the frequency of fading with wind velocity. 
The aim of the present investigation is to verify this result for Ottawa and to 
see how the constant of proportionality (AK) between V and Nd varies with the 
geographical position of the recording station. 
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2. EXPERIMENTAL DETAILS AND RESULTS FOR OTTAWA 


The three-station fading records, used in the present investigation, were 
taken at the Radio Physics Laboratory, Defence Research Telecommuni- 
cations Establishment, Ottawa, using the equipment described by Chapman 
(1953). These records were taken on various frequencies ranging between 
2.0 and 2.6 Mc/s in the daytime hours between 0800 and 1600 hours E.S.T. 
during the period June, 1953, to April, 1955. They have been analyzed to 
determine the horizontal drift velocities in the E region by the method outlined 
by Rao and Rao (1957), using median values of time delays. Average frequency 
of fading has also been determined in each case by the simple method given 
by Rice (1945). The 237 observations, thus obtained, have been divided into 
five frequency groups: 2.0, 2.2, 2.3, 2.4, and 2.6 Mc/s, respectively, according 
to the operating frequency. In each frequency group, the various values of 
drift velocities have been separated into subgroups falling in the ranges 1—10, 
11-20, 21-30, etc., meters/second and the average values of drift velocities 
and fading frequencies in each subgroup estimated. The average values of 
V and N, thus obtained, have been presented in Figs. 1 and 2. Figure 1 contains 
observations taken on 2.0 and 2.6 Mc/s, respectively. A least-square line 
passing through the original and the various experimental points is drawn 
for each of these two frequency groups, as done by Rao and Rao (1958). A 
similar plot of observations, obtained on 2.2, 2.3, and 2.4 Mc/s, respectively, 


(OTTAWA) 


© 2.0 Mc/s 
5 2.6 Mc/s 


FADING FREQUENCY (PEAKS PER MINUTE) 





0 50 100 150 200 
VELOCITY (M/SEC )——___—_> 
Fic. 1. Plot of average N vs. average V for observations at Ottawa, on frequencies 2.0 


and 2.6 Mc/s. An average straight line is shown passing through the origin and the points 
of each frequency group. 
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is shown in Fig. 2. As the total number of observations taken on these three 
frequencies is smaller than those taken on other frequencies, these points are 
randomly distributed about a single least-square line passing through the 
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Fic. 2. Plot of average N vs. average V for observations at Ottawa, on frequencies 2.2, 
2.3, and 2.4 Mc/s. An average straight line passing through the origin and all the points is 
taken to correspond to an average frequency of 2.3 Mc/s. 


origin and all the points. This straight line is taken to correspond to an 
average frequency of 2.3 Mc/s. These N vs. V graphs show a wavelength 
dependence of the gradient of the average straight lines drawn through the 
experimental points, as has already been noticed by Rao and Rao (1958). 
Hence, assuming their general equation 


(1) V = KNX 


the value of A, corresponding to each straight line in Figs. 1 and 2, is estimated 
and shown in Table I. 


TABLE I 


Values of K for each operating frequency 


Operating frequency, 
/ 





Mc/s K 

2.0 2.37 
2.3 2.37 
2.6 2.48 


Average 


2.41 
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The standard error of estimate in the values of K is about 5 to 10%. The values 
of K corresponding to all the three frequencies lie within the limits of this error. 
Using the average K from Table I, the general equation, similar to (1), for 
Ottawa may now be written as 


(2) V = 2.41 Nn. 


3. COMPARISON OF RESULTS FOR DIFFERENT LOCATIONS 


From a comparison of the value of K = 2.41 for Ottawa with that of 
K = 1.86 for Waltair, it is obvious that K varies with the geographic position 
of the recording station on the earth. In order to get a value of K for another 
location, Fig. 16 of the paper by Harang and Pedersen (1956) has been sub- 
jected to the type of analysis outlined above and the resultant plot of the 
average valuestof V and N shown in Fig. 3. 
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Fic. 3. Plot of average N vs. average V for observations at Kjeller, on a frequency 2.0 
Mc/s, taken by Harang and Pedersen. An average straight line is shown passing through 
the origin and the various points. 


In preparing this graph, all the points with values of N greater than 0.75 
cycles/second in the original plot have been omitted for the following reason. 
Harang and Pedersen (1956) have pointed out that higher frequencies of fading 
correspond to smaller sizes of the diffraction pattern on the ground. Following 
Rao and Rao (1957), it may be noticed that these cases involve larger curva- 
tures of the lines of maxima of the diffraction pattern, and that wind velocity 
measurements tend to be inaccurate under such conditions. Further it may 
be noticed that average frequencies of fading higher than 0.6 cycles/second 
have not been used in the graphs for Waltair and Ottawa. 
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A least-square line passing through the origin and the various points in Fig. 3 
has been drawn as in the previous cases. It may be noted that those points 
in Fig. 3 which correspond to a larger number of observations lie close to this 
straight line. These points which deviate most from this straight line corre- 
spond to a smaller number of observations. Using the gradient of this straight 
line, the following equation for Kjeller may be obtained: 


(3) V = 2.29 Nx. 


The operating frequency for the observations at Kjeller was 2.0 Mc/s. 

It has already been pointed out that for the same velocity of drift, smaller 
ionospheric irregularities produce higher frequencies of fading. Pawsey (1935) 
has considered the lower portion of the E region as consisting of ion clouds 
having edges with horizontal dimensions of the order of L/2, where L is the 
least distance on the ground over which the fading is substantially the same. 
This L may be taken as a characteristic length which specifies the size of the 
ionospheric irregularities, following Briggs and Phillips (1950), who have 
defined it as the distance between two points on the ground for which the 
fading correlation is 0.5. Again following Pawsey (1935), the period of fading 
may be approximately given by 2 L/V. Now, from a comparison of this with 
equation (1) we get 
(4) L = (1/2)Ki. 


It may be noted that the observations presented for Waltair, Ottawa, and 
Kjeller have all been made by the same method and during the period of 
minimum solar activity. It would be interesting to compare the results for 
these three locations, which have been summarized in Table II. Calculated 
values of geomagnetic latitudes for these locations have been presented in 


TABLE II 


Comparison of results for Waltair, Ottawa, and Kjeller 


Geographic Geomagnetic Magnetic 
Place latitude latitude latitude L/x A/L 
Waltair 17.8° N. 736° N. 7.6 N. 0.93 1.08 
Ottawa 45.4° N. 56.8° N. 61.2° N, 1.21 0.83 
Kjeller 60.0° N. 60.2° N. 57.8 'N. 1.15 0.87 


this table. These values have been obtained under the assumption of a uniform 
dipole field with the geomagnetic north pole situated geographically at 
78.5° N., 69.1° W. Calculated values of magnetic latitudes for these three 
locations have also been presented in Table II. These values have been 
obtained from the actual magnetic dip angles.* The statistical average sizes 
of ionospheric irregularities (L) for the three locations have been calculated 
in terms of the wavelength \ using equation (4), and presented in Table II. 
It can be seen from this table that the average size of ionospheric irregularities 
probably increases with the actual magnetic latitude. Thus, there seems to be 


*Note that, in this case, magnetic latitude = tan (1/2 tan magnetic dip). 
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some geomagnetic control over the size of the ionospheric irregularities in the 
E region and consequently over the fading frequency. The physical significance 
of this will be discussed in the next section. 

Using a different method, Briggs and Phillips (1950) have obtained.a histo- 
gram of the sizes of ionospheric irregularities observed at Cambridge (geo- 
graphic latitude: 52.2°N.; geomagnetic latitude: 54.9°.N.). From their 
histogram, the most probable size of ionospheric irregularities is found to be 
about 170 meters. Using their operating frequency of 2.4 Mc/s for \ in equation 
(4), it may be found that the most probable L for Cambridge is 1.36 \ which 
is high compared with the values in Table II. In fact it is logical to expect 
this value for Cambridge to be smaller than that for Kjeller in Table II, 
mostly because the observations of Briggs and Phillips (1950) correspond to the 
year 1949 where solar activity was much higher than in 1954, the year of 
minimum solar activity. This discrepancy is probably due to the error in the 
present estimation of ZL through the approximate equation (4). 


4. PROBABLE CAUSE OF VARIATION OF L WITH GEOGRAPHIC 
POSITION ON THE EARTH 

It is well known from meteor studies (Rao 1958; Rao and Armstrong 1958) 
that in the meteoric M region, which usually overlaps the lower portion of the 
E region, prevailing winds are normally turbulent. Ionic irregularities at the 
bottom of the E region of the ionosphere are probably caused as a result of 
such a disturbed or turbulent state. All types of motion in the earth’s atmos- 
phere derive their energy directly or indirectly from the sun (Chapman 1953). 
Hence, it is logical to assume that atmospheric turbulence at the bottom of the 
E region, and hence the size of the irregularities there, are in some way con- 
trolled by the energy of the sun’s radiation. 

Sometimes, turbulence-power w per unit mass has been related to the eddy 
size L, (Booker 1956) by the equation 


(5) L, = #/w 


where v is the velocity of the eddy. The systematic drift velocities in the E 
region appear to have the same range of values at various latitudes 
(Sethuraman 1958). Further, the height of the E region above the ground is 
almost the same at all latitudes, as can be seen from Fig. 4(a). This figure 
shows a plot of the yearly mean values of the virtual height of the E region, 
for 1954, for various ionospheric stations vs. their geographic latitudes. Hence 
it may be reasonably assumed that parameters like the diffusion coefficient, 
which usually vary as a function of height above the ground, are almost the 
same at all latitudes. These parameters and the wind velocity in the E region 
may influence the eddy velocity v. Hence, under the assumption of a constant 
v, we may then consider that 


(6) L. « 1/w. 


Now, it may be assumed that the most probable size of irregularities Z at the 
bottom of the — region at any place is directly related to the most probable 
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eddy size, L., there. Further, the turbulence-power w may be directly related 
to the solar energy absorbed in the E region. As a simple approximation, it 
may be considered that the variation with latitude of the solar energy absorbed 
in the E region at noon is directly proportional to that of the maximum rate 
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Fic. 4. (a) Plot of yearly mean values of virtual height of the E region vs. geographic 
latitude, for the year 1954. 
(b) Plot of yearly mean values of critical frequency of the E region vs. geographic 
latitude, for the year 1954. 


of ion production q(x) at the height h(x) of maximum ionization in the E 
region (x represents the zenith distance of the sun). This may be taken as a 
proportionate measure of the variation of w with latitude. It may be seen from 
Chapman’s derivation (1931) that 


q(x) (_eé 
(7) ? € J = fo 
a mM 


under equilibrium conditions established at noon. In this equation, @ is the 
recombination coefficient, fy is the noon critical frequency of the E region, and 


_ Blocos x 
(8) q(x) = H exp | 


where Jy is the intensity of the sun’s radiation incident upon the E region, 
8 is the number of ions produced by absorption of unit quantity of the radia- 
tion, and His the scale height. It has been found by previous workers(Harnisch- 
macher 1951; Scott 1952; Saha 1953) that the critical frequency of the E 
region at any hour of the day can be represented by a formula of the type 


(9) fe C008 % 
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where the index ” varies with time and latitude. Hence using the same index n 
we may modify equation (7) as 


(10) q(x) = fol” 


because g(x) involves a factor of cos x. Using the above reasoning, it may be 
seen that 


(11) 1/L « fo”. 


In order to verify this relation for the experimental values of \/Z in Table 
II, the following procedure is followed. The yearly mean values of noon 
critical frequencies foE of the E region, for 1954, for various ionospheric 
stations have been plotted against the geographic latitudes of those stations in 
Fig. 4(b). It can be seen from this figure that fo decreases as latitude in- 
creases. Using the mean curve drawn through the various points in this figure, 
the curve of fo! vs. geographic latitude has been drawn in Fig. 5. As the 
Chapman’s theoretical value of 0.25 has been used for m in drawing this 
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Fic. 5. Graph showing the average curve fo! vs. geographic latitude, obtained from Fig. 
4(b). The circles represent values of \/Z from Table II. 


curve, fo! is found to drop off very rapidly with increasing latitude. \/L values 
from Table II have also been plotted in the same figure. The scale for this 
plotting has been adjusted proportionately with the previous scale for the 
fo' curve, in order to get the point for Waltair on this curve. There is not good 
agreement between the curve and the other two points. This may be due 
partly to the use of 2 = 0.25. Hence the variation of m with latitude has been 
considered. As Scott (1952) has pointed out that this variation is a geomagnetic 
rather than a geographic phenomenon, the yearly mean experimental values 
of » for the E region, for 1954, at various ionospheric stations have been 
plotted against their corresponding geomagnetic and magnetic latitudes in 
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Figs. 6(a) and 7(a) respectively. Geomagnetic latitudes in Fig. 6 have been 
obtained using a uniform dipole magnetic field and magnetic latitudes in 
Fig. 7 have been obtained using the actual magnetic dip values. Thanks are 
due to Dr. R. G. Rastogi for furnishing the data, for Figs. 6(@) and 7(@), from 
his unpublished work. It can be easily seen from these figures that the mean 
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Fic. 6. (a) Plot of yearly mean values of the index m for the E region vs. geomagnetic 
latitude, for the year 1954. 
(b) Plot of yearly mean values of critical frequency of the E region vs. geomagnetic 
Jatitude, for the year 1954. 
Fic. 7. (a) Plot of yearly mean values for the index n of the E region vs. magnetic latitude, 
for the year 1954. 
(b) Plot of yearly mean values of critical frequency of the E region vs. magnetic 
latitude, for the year 1954. 
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curve, through the points, is more flattened at higher latitudes in Fig. 7(a) 
than in Fig. 6(a@). However, it may be said that m (£) decreases as geomagnetic 
or magnetic latitude increases. The values of fyE for the various ionospheric 
stations in Fig. 4(6) have been replotted in Figs. 6(6) and 7(b), against geo- 
magnetic and magnetic latitudes respectively. From these figures also, it can 
be seen that the mean curve through the various points is more flattened 
at higher latitudes in Fig. 7(b) than in Fig. 6(b). Now, using Figs. 6 and 7, the 
mean curves of fo!/” vs. geomagnetic and magnetic latitudes have been drawn 
in Figs. 8 and 9, respectively. From a comparison of Figs. 5 and 8 it may be 
said that the effect of using the experimental values of is to cause the value 
of fo!” to decrease more slowly with increasing latitude than when u = 0.25 
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Fic. 8. Graph showing the average curve fo!" vs. geomagnetic latitude, obtained from 
Fig. 6. The circles represent values of \/L from Table IT. 

Fic. 9. Graph showing the average curve fo!” vs. magnetic latitude, obtained from Fig. 7. 
The circles represent values of \/Z from Table IT. 
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is used. Though the curve in Fig. 8 tends to flatten at higher latitudes, there is 
a definite increase of fo!" at higher latitudes shown in Fig. 9, where the actual 
magnetic dip angles are used. This increase is probably due to the ionization 
produced by charged particles from the sun which are more attracted towards 
the magnetic poles. This is also in agreement with the observations of Scott 
(1952), who has pointed out that in the auroral zone ionization by solar 
radiation is modified by other factors like charged particles. The values of 
\ L in Table II have been replotted in Figs. 8 and 9 using the proper geo- 
magnetic and magnetic latitudes respectively. The scales for the Y-axis in 
these figures have been adjusted as in Fig. 5, in order to get the points for 
Waltair on the curve. From either of these two figures, it may be said that 
there is reasonable agreement between the curve and the points, within the 
limits of experimental errors. Thus, in general, we may say that equation (11) 
is verified to the extent that the size of ionospheric irregularities in the # 
region varies with geomagnetic or magnetic latitude in the same proportion as 


fol”. Noon values of fo have been taken as representative values for the whole 


day, whereas values of L have been obtained using observations distributed 
over the whole day. No definite conclusions can be drawn with only three 
points in Figs. 8 and 9. It remains to be seen if observations from other 
latitudes would also fit into the curves of Figs. 8 and 9. It has to be remem- 
bered, however, that the present observations have been confined to the year 
(1954) of minimum sunspot activity. This is probably a good reason why it 
has been possible to obtain a definite linear relation between V and N. It also 
remains to be seen if a similar relation could be obtained with observations 
during the maximum of sunspot activity. It may be expected that during 
periods of higher sunspot activity, the bottom of the E region would be more 
disturbed and turbulent, resulting in smaller values of the average size of 
ionospheric irregularities at each place. Further, as geomagnetic storms are 
more frequent during the year of maximum sunspot activity, there may be 
a wider scatter of points in the N vs. V scatter diagram more towards the 
higher frequencies of fading (Harang and Pedersen 1956). Actually, the scatter 
of points may be so wide that the average size of the irregularity at any place 
may lose the significance given to it in this paper. 

Physically, the effect of a strong geomagnetic storm on N may be envisaged 
as follows. A strong geomagnetic storm may be expected to disturb the eddy 
currents in the E region. As a strong geomagnetic storm is always accompanied 
by a greater amount of solar energy input than normal, into the F region, 
it is logical to believe, from equation (6), that the most probable eddy size 
should decrease. Hence a decrease in the size of ionospheric irregularities in the 
E region may be expected. This, in turn, would increase the frequency of 
fading. Chapman (1953) has found that during geomagnetically disturbed 
conditions, there was an indication of great increase in the drift velocity, 
observed mainly through sporadic E echoes. It is generally known that drifts 
of the sporadic EF layer are different from those of the normal E layer, and 
further wind velocity determinations under disturbed conditions are rela- 
tively unreliable. Hence, no definite conclusions could be drawn about the 
effect of a geomagnetic storm on wind velocity in the E region. Even if this 
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effect is to increase the drift velocity, it would add to the effect of decreased 
ionospheric irregularities in enormously increasing the frequency of fading. 
This is in agreement with Harang and Pedersen’s (1956) important observa- 
tion that the frequency of fading is unusually high during periods of strong 
geomagnetic storms. It would be interesting if more light could be thrown 
upon this problem of relating ionospheric irregularities with atmospheric 
turbulence in the £ region by future investigators. In conclusion, it is suggested 
that the average size of the ionospheric irregularities varies in some smooth 
manner with geomagnetic or magnetic latitude. 
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TRANSIENT DECREASES IN COSMIC RAY INTENSITY DURING 
THE PERIOD OCTOBER 1956 TO JANUARY 1958! 


A. G. FENnTon,? K. G. McCracken,’ D. C. ROsE, AND B. G. WILSON! 


ABSTRACT 


Recent nucleon intensity data obtained from high counting rate recorders 
at Ottawa and Hobart, and subsidiary stations, have been examined for evidence 
for the superposition of transient decreases. It is concluded that, with the statis- 
tical accuracy now available due to the high counting rates, it is possible to 
distinguish two types of transient decreases in the observed variations, super- 
imposed upon the slower 11-year intensity changes. One of these is an almost 
symmetrical event lasting up to 2 weeks and exhibiting a recurrence tendency of 
about 27 days, while the other is the more abrupt Forbush decrease which recovers 
over a period of several days. The evidence indicates that the intensity-controlling 
mechanism responsible for these short-term transient changes is able to influence 
the cosmic ray flux at the earth independently of other events that may be in 
progress at the time. There is also evidence that the physical process controlling 
the Forbush type of decrease operates over a volume large compared with the 
earth because the intensity changes at places as far apart as Ottawa, Canada, 
and Hobart, Tasmania, show changes that are the same within the accuracy of 
the measurements. 


INTRODUCTION 


In a recent paper (Fenton, Fenton, and Rose 1958) it was shown that after 
correction for atmospheric effects the monthly mean intensity recorded by 
cubical meson telescopes and neutron monitors at Ottawa and Resolute de- 
creased appreciably between 1954 and 1957, a period during which the solar 
activity increased from an exceptionally low level in 1954 to an unusually 
high level in 1957. It was also shown that the short-term transient decreases 
observed from time to time in the daily mean intensities are not themselves the 
cause of this longer term decrease, but that they are superimposed upon it. 
It is the purpose of the present report to consider the question of superposition 
of decreases in greater detail in an attempt to determine some of the important 
features of the mechanisms controlling the level of the cosmic ray intensity 
at the earth during the period October 1956 to January 1958. To do this a 
separation of the intensity changes into well-known types has been carried out 
and some comparisons made between corresponding intensity changes at two 
stations, one in the northern and the other in the southern hemispheres, both 
above the latitude knee. 


APPARATUS AND LOCATION OF STATIONS 


During the period considered in this paper, high counting rate neutron 
monitors and meson telescopes were operated at Ottawa and Hobart and at 


1Manuscript received February 6, 1959. 

Contribution from the Division of Pure Physics, National Research Council, Ottawa, 
Canada. 

Issued as N.R.C. No. 5160. 

2Visiting scientist from the Department of Physics, University of Tasmania. 

2Department of Physics, University of Tasmania. 

‘National Research Council Sulphur Mountain Laboratory, Banff, Alberta. 


Can. J. Phys. Vol. 37 (1959) 
569 











570 CANADIAN JOURNAL OF PHYSICS. VOL. 37, 1959 


several other stations in both the northern and southern hemispheres by the 
National Research Council and the University of Tasmania respectively. These 
recorders were developed for use during the International Geophysical Year 
and results from them are considered to be directly comparable, since no 
significant differences in response or cone of acceptance are expected to arise 
as a result of the slight differences in design and construction. The location of 
recording equipment operated at the present time by the two groups is set 
out in Table I. A complete list of the stations is given for reference purposes, 
although not all of these are involved in the present analysis. All equipment 
is operated near sea level except in the case of Sulphur Mountain (2283 meters) 
and the Hobart neutron monitor (725 meters). 


TABLE I 

Station Geographic co-ordinates Geomagnetic latitude 
Sulphur Mountain 51° N., £15.67 W. 58° N. 
Churchill 58.8°N., 94.2° W. 69° N. 
Hobart 42.8°S., 147 .b° E. 52°'S. 
Lae (New Guinea) 6:7°S,, 14775. 16°S. 
Macquarie Island 54.5°S., 158.9° E. 61° S. 
Mawson (Antarctica) 67.6°S,, 62:97. te 3: 
Ottawa 45.4°N., 75.6° W. 57° N. 
Resolute 14.0° Ni. SELF: 83° N. 

NUCLEON INTENSITY VARIATIONS AT OTTAWA AND HOBART 


The Ottawa and Hobart daily mean, pressure corrected, nucleon intensities 
have been compared for the period October 1956 to January 1958 (Fig. 1) 
in order to illustrate the world-wide nature of the intensity variations. Nucleon 




















4a, i; | | | | 
a vid | | | | 
age Vs | | | | 
iy I | | | | 4 
Cn oa. | | | | } | } | 
uy : | | | 
Wea TL A HOBART NUCLEON | | 
; |* yf } heath Ahaha | ant al | f 
Lt, | TR LAME tae AD eo | fy 
ae, da! ow TF Me 1 Viiv YW fal 
no | i | if] i tf 1} i] | 4 Nn 
\ | | | Ie | 4 f iJ | H b if \ ’ fh 
ae 1 AR F oI | iy | +f] fe ak 
QO} hi RPL | A Hl if |e] 
g V WV 1 i . ¢ i Lat 3 a A i | pi iy é 1 
=} . ah Natl fir Th 4, | | +] 
5 Wi | ei TS ety | ' 
| WU a ! [tis = ‘ nf } A fey | 
5 | {} | F [a Z| | Lb }4 | iy | | 
ot 1 1 4b fp doy MR a ie 
| ieee a2 eS } W | 
| iv | 
| | |} | or mawa. NucLeoN | | | 
ae Gi 1 | | | | | 
| | A | 
| | | | | | 
1956 | | 1957 | | | | | | 1958 
oct | Nov DEC, JAN, | FEB | MAR | APR | MAY | JUNE | JULY | AUG | SEPT | ocr | NOW| DEC 
Fic. 1. Daily mean nucleon intensities at Ottawa and Hobart, October 1956 to January 


1958. Intensity scales arbitrarily displaced to avoid overlapping. 
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intensities are considered to be more suitable for the present study than meson 
data, since the variations are much larger in the nucleon than in the meson 
component, and because of the atmospheric influence remaining in the meson 
data after correction for pressure alone. The energy dependence of the in- 
tensity variations will not be considered at this time. In the figure the scales 
have been chosen so that the intensities are comparable on a_ percentage 
basis, and the horizontal lines carrying the time scale provide a constant 
reference level. The Ottawa data have been arbitrarily separated from the 
Hobart data to avoid overlapping. Throughout this paper Greenwich mean 
time will be used, the daily mean intensity referring to the mean hourly 
intensity over a 24-hour period from 00 to 24 hours G.M.T. 

Reference to Fig. 1 shows that there is a close correlation between the two 
sets of data, both in the over-all change of intensity from the beginning to the 
end of the period, and in the short-term fluctuations. The intensity on January 
31, 1958, at both stations was approximately 15% lower than the high point 
reached at the end of October 1956. (A more detailed comparison is presented 
in the next section.) The slow over-all intensity decrease was fairly rapid 
from the beginning of November 1956 to January 1957, then it appears to have 
recovered slightly before a further general decrease from August 1957 onwards. 
This over-all decrease is identified with the 1ll-year solar cycle intensity 
changes reported originally by Forbush (1954) and discussed recently in 
relation to the present cycle of solar activity (Fenton, Fenton, and Rose 1958; 
Neher and Anderson 1958; and Lockwood 1958). Close inspection of Fig. 1 
reveals that there are relatively small differences between the daily mean 
intensities observed at the two stations (the standard deviation for 1 day, 
calculated from the counting rate, is 0.15% and 0.11% at Ottawa and Hobart 
respectively). An examination of the hourly intensities at these and other 
stations has revealed obvious differences in response to transient decreases 
from station to station. These appear to be at least partly due to primary 
anisotropies which are detected at somewhat different times and directions 
in space as the earth rotates. This aspect of the variations has been investi- 
gated, and the results will be reported later. 

The relatively short-term intensity changes shown in Fig. 1 seem to be more 
satisfactorily explained in terms of a succession of transient decreases below 
the prevailing level determined by the 11-year intensity cycle than by a 
series of temporary increases above a lower level. Thus while the maxima of 
the intensity-time plot lie on a smooth curve, indicating that a decrease was 
usually followed by a return to near the predecrease intensity, the minimum 
values show no such systematic behavior. Furthermore, as is clearly evident 
for the months October 1956 to January 1958 (Fig. 1), the maxima tend to be 
broad and relatively free from rapid intensity changes, whereas the minima 
are short-lived and involve appreciable rates of change of intensity. These 
observations suggest that the decreases represent the disturbed condition. 
This view is supported by reference to geomagnetic disturbance data, which 
shows that sudden commencements of magnetic storms often occur near the 
beginning of sudden decreases in the cosmic ray intensity, suggesting that 
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the cosmic ray effect is due to a screening or deceleration process at these 
times rather than to a source or acceleration process occurring at other times. 
Additional arguments have been presented by others (e.g. Van Heerden and 
Thambyahpillai 1955; Neher 1957) to show that decreases rather than in- 
creases are involved, except in the rare cases of solar flare increases, which 
represent injection of additional particles into the solar system by the sun. 
Other evidence supporting the view that the variations during the period 
under discussion consisted only of superimposed transient decreases is provided 
by the correlation of these events with the central meridian passage of active 
solar regions (Venkatesan 1959). 


SUPERPOSITION OF INTENSITY DECREASES 


The Hobart section of Fig. 1 corresponding to the months of November 
1956 to February 1957 has been reproduced in Fig. 2, A. Three classes of 
intensity change are apparent as shown by the straight lines from maxima 
to maxima, the smooth curved lines showing more or less symmetrical dips 
in intensity and the residual smaller variations. This division is somewhat 
arbitrary but it can be repeated in a consistent way at different stations. 
The first class is the slow 11-year cyclic change in intensity; the second is a 
succession of loops of intensity change, which for this series of data appears 
to. be almost periodic and may be associated with solar rotation. The third 
class consists of the residual short-period changes lasting a few days, and the 
large Forbush events. 

The earlier work of Fenton et a/. (1958) and Venkatesan (1959) suggests 
that, of the above classes, two and three may be due to the same basic 
phenomena but are quite different from the 11-year variation because they 
have a different spectral sensitivity. In Fig. 2, B, both the 1l-year variation 
and the variations of duration about a month have been removed and the 
residual variations probably represent small and large Forbush events. 
Figure 3, A and B, show the same analysis covering the same period at Ottawa. 
It is, of course, quite possible to continue this analysis for longer periods. 

The similarity between the curves 2B and 8B is strong evidence for the 
wor!d-wide nature of the modulating mechanism. Attention is directed to the 
close correspondence in all but the smallest variations. On the other hand, 
there are differences and if these are significant they are very important in 
assessing the value of theories about the modulating mechanism causing the 
variation. 

Let us look first for differences in the 11-year cyclic change at the two 
stations. A careful examination makes it apparent that the straight lines 
joining the points of maxima have not exactly the same slope and the intensity 
drop may not have the same ratio at the two stations. A quantitative com- 
parison has been made in two ways. First, the monthly means can be com- 
pared but these might be affected by differences in the magnitude of decreases 
in the Forbush events during the month. Alternatively, we may make a 
comparison of the mean level of the straight lines joining the maxima as 
shown in Figs. 2 and 3. This was done by picking three periods in Fig. 1 (the 
straight lines representing the 11-year variation are not shown) and examining 
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the ratios of the mean levels of the lines. Table II shows these normalized to 
intensities of 100% for the first period chosen. 


SCALED COUNTS PER HOUR 


SCALED COUNTS PER HOUR 





TABLE II 


Intensity ratios 





Period centered on: Hobart Ottawa 
November 18, 1956 100% 100% 
July 3, 1957 92.2% 91.8% 
December 21, 1957 87.0% 87.3% 
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Fic. 2. (A) Daily mean nucleon intensities at Hobart, November 1956 to January 1957 

inclusive, shaded to indicate division into three types of decreases. 
(B) A section of (A) with the slower changes removed. 

Fic. 3. (A) Daily mean nucleon intensities at Ottawa, November 1956 to January 1957 

inclusive, shaded to indicate division into three types of decreases. 
(B) A section of (A) with the slower changes removed. 
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Monthly means treated in the same way are shown in Table III. 


TABLE III 
Intensity ratios 








Hobart Ottawa 








Month 
October 1956 100% 100% 
June 1957 89.4% 89.6% 
January 1958 85.3% 85.8% 





The decrease in intensity is very nearly the same at both stations. It is 
considered that the differences in columns 2 and 3 are not significant in either 
table even though they are somewhat greater than the statistical differences 
that might be expected from a Poisson distribution in counting rate. For 1 
month of data the statistical difference due to counting rate should be about 
03% but the stability of neutron monitors with their amplifiers and dis- 
criminators and the accuracy of barometer correction are also involved. 
Although calibrations are carried out regularly, an over-all confidence of better 
than 0.5% over a year is difficult to achieve. There will be a further discussion 
of accuracies in the following paragraphs. Our conclusion here is that the 
fractional long-term intensity change is the same at the two stations within 
the order of accuracy of the measurements. 

The depth of the Forbush events or the deepest point in the more sym- 
metrical decreases is a measure that might be used in looking for differences 
between the two stations. A technique for comparison of these is shown in 
Fig. 4. The quantity 6//J can be compared for similar events at the two 
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Fic. 4. An idealized Forbush type of decrease showing how 6//Z was measured. 


stations. The ratio (6///) ottawa/ (61/1) nopare Should be a constant if the modu- 
lating effect is truly world-wide in each event. This is particularly so since both 
stations are near or above the geomagnetic knee. Table IV shows a series of 
events all of which are included in Fig. 1. Before considering these let us 
examine the expected accuracy of such a comparison. The quantity 6//TJ is 
estimated as in Fig. 4 from the position of the straight line and the depth of 
the lowest point representing a daily average. Its statistical accuracy, there- 
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fore, can be derived from that of the daily averages. This was stated earlier 
to be 0.11% for Hobart and 0.15% for Ottawa. 


TABLE IV 
‘The magnitude of 11 Forbush decreases at Ottawa and Hobart 
and the ratio between them 











Ottawa, Hobart, 





Date 6I/I 61/1 Ratio 
November 17, 1956 .079 .082 .96 
December 8, 1956 .043 041 1.05 
January 3, 1957 .069 .061 1.13 
January 22, 1957 .123 113 1.08 
April 20, 1957 .088 .075 je: 
July 3, 1957 .080 .081 .99 
August 5, 1957 .061 .070 87 
October 23, 1957 .088 .094 .93 
November 27, 1957 .089 .100 .89 
December 22, 1957 .061 .062 .99 
January 22, 1958 .038 .041 94 


This was based on the total counts for 24 hours and a Poisson distribution. 
There has always been some uncertainty about the validity of using this to 
estimate accuracy for a number of reasons. These include the fact that there 
are errors in barometric readings and possibly in the correction coefficients 
used, and because in the neutron monitor an appreciable number of events 
result in collection and counting of more than one neutron per incoming 
nucleon. 

A simple examination of the latter effect was carried out using data from 
the Ottawa neutron monitor during February and March 1957. The Ottawa 
neutron has three roughly equivalent recording channels. The final figure used 
is the sum of the three labelled A, B, and C. The procedure of data handling 
makes it quite easy to obtain the ratio of the daily averages of say /,/I¢ 
before barometer correction. On the assumption that the same flux goes 
through the two parts of the pile, and no experimental changes occurred 
during the period, the ratio of J,4/Jc¢ should be constant except for random 
fluctuations. A and C were chosen rather than A and B or B and C since A 
and C are farther apart. The standard deviation of J,//¢, that is, the r.m.s. 
value of deviations from the mean was calculated from the actual data. Since 
I, and J¢ are practically the same and J, is only slightly higher, it is easy to 
show that the quantity |o(/,4/Jc)}/(/2 V3) is a fair value for the standard 
deviation of the sum of the three channels. Using 2 months’ data of February 
and March for 1957, selected because fluctuations during those 2 months 
were relatively less than for other parts of the period studied, the derived 
value of o is 0.22%. This is slightly higher than the 0.15% derived from the 
counting rate. Added to this there is a barometric error estimated at 0.1 mb 
and a lack of certainty about the best barometer coefficient. These we estimate 
to represent an additional random error of 0.1%; thus, we may take the 
best estimate of our error in the nucleon intensity as being 0.24%. The average 
value of 6J/J in Table IV is 7.4% and assuming that the error in Hobart is 
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the same as that at Ottawa a simple calculation of the combination of errors 
gives a standard deviation for column 3 in Table IV of 6.4%. The actual r.m.s. 
deviation from the mean of the ratio is 9.3%. This is not significantly different 
from 6.4%, the best estimate we can make of the expected deviation. One can 
conclude, therefore, that the modulating effect on the cosmic ray flux is so 
nearly the same at Ottawa and Hobart that measurements of this type 
cannot show differences. In support of this conclusion it has recently been 
shown (McCracken 1959) from analysis of the data from Hobart and Mawson 
that radiation arriving from directions close to the equatorial plane is affected 
to the same extent as that arriving from large angles to this plane. This is also 
confirmed by an examination of the curves B in Figs. 2 and 3. The small 
variations are as closely in agreement as could be expected in view of the 
statistical accuracies involved and the rather arbitrary method of separating 
the types of variation. On the other hand, a more detailed analysis must show 
differences in the decreases because it is well known that Forbush decreases 
are sometimes followed by an enhanced diurnal variation which follows local 
time. This alone would be expected to result in differences at different stations. 
for a preliminary examination of some of these we should look at a more de- 
tailed analysis of some of the Forbush-type decreases using hourly values. 

In this case we use data from Sulphur Mountain since it is the one in our 
group with the best statistics, having a counting rate of about 100,000 per 
hour. As an example, the hourly intensities observed at Sulphur Mountain 
during the period August 3 to 10, 1957, are plotted in Fig. 5, (a). It is evident 
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Fic. 5. Hourly nucleon intensities, August 3 to 10, 1957. 
(a) Sulphur Mountain. 
(6) Average of Hobart, Ottawa, Mawson, and Sulphur Mountain. 
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that the variations are significantly greater and more regular than expected 
from statistical fluctuations. This event appears to have a different form from 
station to station, and in order to obtain an approximation to the effect that 
would be observed with an omnidirectional (47) detector at the position of 
the earth (i.e. independent of the earth’s rotation) the mean hourly data from 
four stations widely spaced in longitude are presented in Fig. 5, (0). In the 
averaging process each station was given equal weight so that percentage 
changes would be comparable. It is clear that, after averaging, most of the 
fluctuations apparent in the record from a single station have been removed, 
leaving a rather smooth transient decrease on which is superimposed another 
decrease beginning abruptly on August 4 and ending on August 6. It seems 
likely that variations of the type shown in Fig. 5, (a), are to be identified 
with the enhanced daily variations discussed by others (Elliot 1952). 


DISCUSSION 

The evidence presented above supports the view that the cosmic ray inten- 
sity observed at the earth is controlled by separate events and that the main 
controlling influence acts in the same manner over great distances. In other 
words the scale of the modulating effect must be larger than the earth's 
dimensions because at stations separated as far as Ottawa and Hobart the 
effect is the same within the accuracy limits of current measurements. This 
applies both to the ll-year cyclic change and the magnitude of short-term 
decreases of the Forbush or cosmic ray storm type. It is significant that the 
individual short-period events appear to be completely independent and not 
very closely related to the over-all ll-year cyclic change, though the con- 
trolling mechanism of the two may differ only in scale as suggested by Parker 
(1958). It is difficult to classify the short-term decreases lasting from a few 
hours to perhaps 2 weeks or | month on the basis of their magnitude or shape 
since the shape is extremely variable and it is uncertain whether the more 
symmetrical decreases taking several days to reach a minimum are sub- 
stantially different from those which last from a few hours to 2 or 3 days and 
show substantially greater gradients. The former may be just a number of 
smaller decreases so close together that they appear to be one of gentler slope. 
This is in agreement with the analysis of Venkatesan (1959) in comparing 
Forbush decreases with the central meridian passage of active regions on the 
sun. 

An analysis by one of us on the energy dependence and correlation with 
geomagnetic effect (McCracken 1958) suggests that both are produced by a 
common (probably geocentric) process and the difference in appearance of 
the events is due to the differences in the extent or speed of the cloud or 
stream responsible for setting up the process. 

The location of the intensity-reducing agencies relative to the earth is 
important in any theory of the process responsible for cosmic ray intensity 
variations. The rapid onset and brief duration of the Forbush decreases, the 
correlation with active solar regions crossing the central meridian, and the 
association with concurrent geophysical phenomena suggest that these events 
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are produced by a process acting relatively closely to the earth. On the other 
hand the similarity in effect over great distances on the earth’s surface suggest 
that the scale is large compared with the earth’s dimensions. Individual 
isolated examples of the rapid decreases which might give a time scale either 
of the rate of decrease or recovery, practically do not exist, that is, there is no 
unique obvious characteristic shape on which to base a theory. 

The slow changes that make up the 11-year variation may be due to the 
cumulative effect of solar matter streaming into the whole solar system 
from active regions on the sun. This would agree with the observation (Fenton, 
Fenton, and Rose 1958) that, although the long-term intensity is clearly 
controlled by the general level of solar activity, there is not a close correlation 
with changes in this activity from month to month. A possible mechanism 
through which the long-term cosmic ray intensity at the earth may be con- 
trolled by matter streaming outward from the sun has been put forward 
recently by Parker (1958). Whatever the mechanism may be, it is quite 
effective during the present phase of solar activity in reducing the intensity 
below the 1954 level and maintaining it without breaks. Between October 1956 
and January 1958 the over-all intensity level fell by about 15% (Tables II 
and III) which is a larger decrease than observed in any of the transients 
during this period. As reported previously (Fenton, Fenton, and Rose 1958) 
the. over-all decrease in nucleon intensity from April 1954 to December 1957 
was 22%, making the slow 1l-year cycle changes the most substantial of the 
three types of events discussed in this paper. 
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INVESTIGATION OF DEFORMATION PROCESSES IN ARMCO 
IRON BY MEANS OF INTERNAL FRICTION AT MEGACYCLE 
FREQUENCIES! 


W. J. Bratina? 


ABSTRACT 


Internal friction, employing frequencies in the megacycle range, has been used 
in order to study the characteristics of elastic and plastic deformation in Armco 
iron. Tensile stresses only were applied in the present investigation. 

Among the characteristic features observed three were investigated in detail. 
(1) A decrease in ultrasonic attenuation at external static stresses within the 
elastic range, which appears to be related to the ordering effect of the stress on 
magnetic domains, was investigated. (2) An initial increase in attenuation values 
regularly observed in cold-worked specimens upon the application of an 
external load was shown to be associated with the internal stresses. (3) The 
changes in attenuation resulting from plastic deformation were found to be 
insignificant in some instances and large in others. An explanation employing the 
magnetic domain theory and dislocation mechanisms is suggested. 


INTRODUCTION 


Attention has been focused recently on the phenomena associated with the 
elastic deformation of metals and the early stages of plastic deformation in- 
cluding the relaxation processes related to them. To study the fine details 
of deformation characteristics in metals, high-frequency internal friction 
methods appear to be particularly well suited. In one of these methods, a 
beam of ultrasonic energy is sent through the specimen, and the attenuation is 
studied as a function of a static strain or stress. The experimental data re- 
ported in the present paper were obtained by this particular technique applied 
to Armco iron specimens. Static tensile strains employed were small, mostly 
in the range of a few tenths of 1%. 


MATERIAL AND EXPERIMENTAL PROCEDURE 

Armco iron used in all investigations was in the form of cylindrical specimens 
with a diameter of 0.560+0.001 inch. The length of the specimens was 10 
or 12 inches. The end faces were machined and polished afterwards to a high 
degree. 

The annealing was done in all cases in a hydrogen atmosphere. Details of 
the annealing procedure are described later. 

A standard Sperry Reflectoscope with some modifications was employed 
for the ultrasonic measurements, and frequencies of a few megacycles (2.2 
megacycles in most cases) were used. At the time of the experiments, the 
Sperry Ultrasonic Attenuation Comparator, which allows the attenuation 
values to be read directly, was not available. The observations of back reflec- 
tions on the reflectoscope screen were recorded by means of a 35-mm camera 

'1Manuscript received September 30, 1958. 

Contribution from the Department of Engineering and Metallurgy, Ontario Research 


Foundation, Toronto, Ontario. 
2Ontario Research Foundation, Toronto, Ontario. 


Can. J. Phys. Vol. 37 (1959) 


579 











580 CANADIAN JOURNAL OF PHYSICS. VOL. 37, 1959 


with exposure times of 1/10th of a second or less. The quartz crystal which 
acts as a transmitter and as a receiver was attached to one end face of the 
specimen. Dow Corning 200 fluid with a viscosity of 12,500 centistokes at 
25° C, as used also by Hutchinson and Filmer (1956), proved to be a highly 
satisfactory coupling liquid between the crystal and the specimen. The ultra- 
sonic attenuation obtained cannot, therefore, be regarded as an absolute 
value; the losses due to coupling liquid are included in the observed attenuation 
magnitude. 

The external tensile stresses were applied by means of Riehle Universal 
Screw Power machine using specially designed auxiliary grips in order to pre- 
vent excessive indentations. The rate of loading and unloading was kept con- 
stant through the entire elastic region. 

The strains were measured and recorded continuously by an MA-5P Micro- 
former-type recorder in conjunction with a T-1M Extensometer. 

The standard practice followed in the experiments was to strain the speci- 
men up to a predetermined strain value; the specimen was then kept in this 
strained condition for 10 minutes and unloaded afterwards. Measurements of 
ultrasonic attenuation were made at intervals during these operations and 
were continued at zero load for 60 minutes after unloading. 


EXPERIMENTAL DATA 
Internal friction, of which ultrasonic attenuation is a particular case, can 
be defined by the equation, 
tan 6 = log dec./z, 


where the logarithmic decrement is equal to the difference between the 
natural logarithms of the amplitude of successive swings. The ultrasonic 
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_ Fic. 1. Logarithm of the amplitudes of successive reflections plotted versus their number 
for an Armco iron specimen. 
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attenuation values shown in the present work were derived simply from the 
difference between the natural logarithms of the amplitude of successive back 
reflections. Figure 1 corresponds to an actual case; the amplitudes of the 
successive back reflections are plotted on semilogarithmic scale. The value of 
the slope representing attenuation value was taken directly from the graph or, 
as a refinement, it was calculated by the method of least squares. 

The ultrasonic technique has been analyzed by Seki, Granato, and Truell 
(1956) and by Redwood (1957). Our technique corresponds closely to that 
discussed by Redwood, the diameter of the quartz crystal being larger than 
that of the specimen. Therefore, the end face of the specimen is completely 
covered by the crystal. Furthermore, the specimen being in the form of a 
cylinder of circular cross section acts as a solid ‘‘wave-guide’’, which accounts 
for an unusually regular attenuation pattern. Dispersion effects cannot be 
avoided completely, however, and the decay pattern is not ideally exponential. 
The deviations are not large and a relatively good straight line could be drawn 
through the experimental points in all cases (cf. Fig. 1). 

The attenuation values are shown in all graphs plotted against the tensile 
strain. The corresponding stress-strain curves are also included in all cases. 

The data for three specimens—which seem to show characteristic features 
are dealt with in this paper. These actually represent the three prototype 
experiments, and the data shown are supported by a large number of observa- 
tions made on many similar specimens. All the specimens in question were 
annealed in hydrogen for 2 hours at 850°C. Specimen No. 1 was cooled 
slowly in the furnace, Nos. 2 and 3 slightly more rapidly. Microphotographs 
show a structure of relatively regular grains with an average grain diameter 
of 0.08 mm. 

Specimen No. 1 was strained in four steps, the strains 0.04, 0.10, 0.33, and 
1.9% being applied in succession. In the first step (Fig. 2(a)), the metal, on 
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Fic. 2. Ultrasonic attenuation as a function of elastic and plastic strain for specimen No. | 
for the first three stress cycles. Stress-strain curves also shown. 
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the basis of stress-strain data, still behaves elastically. As the strain is con- 
tinuously increased from zero value, a rapid decline in the attenuation magni- 
tude—preceded by a slight initial increase—is experienced. A different path 
is followed on unloading* and a kind of hysteresis results. 

During the second and the third stages of straining up to 0.10 and 0.33 
strain per cent respectively (Figs. 2(b) and 2(c)), the characteristic decline of 
attenuation-strain curve in the elastic region ceases, as soon as plastic deforma- 
tion starts to occur. Surprisingly enough, the phenomena taking place during 
this deformation do not seem to have any significant effect on the attenuation 
values. Neither does the shape of the unloading curve differ from the curve 
obtained on loading in the elastic region. 

Finally, a much higher strain of 1.9% was employed (Fig. 3). A gradual 
rise in the attenuation associated with irregular fluctuations is observed as the 
strain increases progressively. The unloading curve assumes also a different, 
characteristic shape. 

In Fig. 4 the data for specimen No. 2 are shown. In contrast to the previous 
case, this specimen was strained beyond the elastic region in one step. Quite 
different results were obtained as compared with Fig. 2. The initial loading 
causes the attenuation to rise at first and only after a certain stress has been 
reached does it start to decrease rapidly. Further stressing beyond the elastic 
limit causes large changes to occur, and the unloading curve obtained after a 
total strain of only 0.80% assumes a character somewhat similar to that of 
the unloading curve shown in Fig. 3. 


TABLE I 








Specimen Total strain (%) Fig. 
No. Anneal and type of deformation No. Attenuation characteristics 
1 850°C, 2hrin (1) 0.04, elastic 2(a) Decrease on loading; increase on 
Ho, cooled in unloading 
furnace (2) 0.10, plastic 2(b) Decrease on loading, constant dur- 
ing plastic deformation; increase 
on unloading 
(3) 0.33, plastic 2(c) The same as strain cycle 2 
(4) 1.90, _ plastic 3 Decrease on loading, fluctuations 


during plastic deformation; in- 
crease followed by decrease on 


unloading 
2 The same as C1038. plastic 4 Increase followed by decrease on 
specimen 1 loading, increase and fluctuations 
except for during plastic deformation; in- 
faster cool- crease followed by decrease on 
ing unloading 
3 The same as (1) 0.18, plastic 6 Increase followed by decrease on 
specimen 2 loading, large increase during 
plastic deformation; decrease on 
unloading 


Increase followed by decrease on 
loading, little change during plas- 
tic deformation; increase on un- 
loading 


(2) 0.35, plastic 


*The terms ‘loading and unloading curves"’ are used to denote ultrasonic attenuation— 
strain relationship throughout. 
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Fic. 3. Ultrasonic attenuation plotted versus large plastic strain for specimen No. 1. 
The data corresponding to the fourth stress cycle are shown. Stress-strain curve included. 
Fic. 4. Ultrasonic attenuation versus strain for the first stress cycle of specimen No. 2. 


Separate stress-strain curve shown. 
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Further experiments, on specimen No. 3, were done partly in order to test 
the suggested explanations, and partly in order to obtain additional informa- 
tion. The respective experimental data (Figs. 6 and 7) are for these reasons 
shown and discussed separately at the end of the paper. 

No significant changes in the attenuation were observed during the 10- 
minute period with the specimens at maximum load. The same happened 
during the 60-minute period with the specimens at zero load. This applies to 
all the specimens and the strains employed in the present investigation. De- 
tailed information on these stages is not, therefore, presented. 

Plastic deformation was in some instances still progressing after the move- 
ment of the tensile machine was stopped. The load was dropping continuously 
as the elongation increased. This accounts for the break in the end of the 
stress-strain curve noticed on some of the graphs. 

The experimental observations are summarized in the above table. 


ANALYSIS OF RESULTS 

The following experimental observations will require an explanation: 

1. An initial decline in attenuation is observed, in annealed specimens as a 
rule, as the strain is continuously increased within the elastic limit. 

2. In cold-worked (plastically strained) specimens, however, a characteristic 
feature of elastic deformation appears to be an increase in attenuation, 
followed by a decrease afterwards (the unloading curves in Figs. 3 and 4). 

3. A significant rise in the attenuation occurs at the yield stress region in 
some cases. 

4. A change in the initial attenuation (at zero load) is observed in speci- 
mens after a plastic or sometimes even after an elastic deformation. 


(a) Elastic Deformation 

The gradual decline in the attenuation associated with an elastic tensile 
strain (Figs. 2, 3, and 4) appears to be determined mainly by the magnetic 
characteristics of the material. The ferromagnetic under investigation is sub- 
jected to the ultrasonic vibrations of megacycle frequency. There is a loss in 
internal energy associated with mechanical vibrations of any material. The 
internal losses in ferromagnetics are much higher than in non-ferromagnetic 
materials for the reasons of phenomena of magnetic nature. The most signi- 
ficant of these are manifested by magnetomechanical hysteresis, macroeddy 
and microeddy currents. Since our specimens were demagnetized and high- 
frequency ultrasonic waves were employed in these experiments, it is likely 
that microeddy currents comprise the main source of energy loss. These currents 
are generated between the adjacent domains due to the domain wall motion, 
induced in our case by ultrasonic strain. Any change in domain configuration, 
in the number of domains or domain walls, and in their relative mobility will 
be reflected on the ultrasonic attenuation magnitude. 

According to the domain theory of ferromagnetism, a ferromagnetic crystal 
consists of magnetic domains, which should be distributed at random in an 
annealed and unmagnetized material. A static stress, applied in our technique 
by means of a tensile machine, would tend to orient the domains parallel to or 
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perpendicular to the direction of stress, depending on whether the magneto- 
striction of the material is positive or negative. At low stresses the reorientation 
is achieved by the growth of the more favorable oriented domains at the 
expense of those less favorably oriented. (Higher stresses might induce domain 
rotation.) The domain dimensions increase under stress and the number of 
existing domains and the domain walls is reduced accordingly; the domain 
wall motion as induced by ultrasonic strain is reduced also. As a consequence, 
a decrease should be experienced in the losses associated with the domain 
walls and their mobility. The microeddy current losses and the ultrasonic 
attenuation should, therefore, decline upon the application of static elastic 
strain. 

In order to test the proposed mechanism, the effects of an external magnetic 
field were studied. A magnetic field plays a role similar to that of stress and 
the application of a field to a ferromagnetic specimen produces redistributions 
and rotations of domains similar to those produced by the stress. The data 
shown in Fig. 5, curve b, where an external magnetic field was applied— 
its axis coinciding with the longitudinal axis of the specimen—furnishes 
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Fic. 5. The effect of an external magnetic field on the ultrasonic attenuation versus strain 
values: (a) without a magnetic field, (6) with a magnetic field. 


support to the previous interpretation. The attenuation at zero load in the 
presence of an external field is already reduced to a value below the minimum 
achieved in the absence of a field by the highest elastic strain employed (Fig. 
5, curve a). This indicates the magnitude of the magnetic contribution. It is 
reasonable to expect that only insignificant changes in attenuation would occur 
on the application of the load, i.e. by purely non-magnetic phenomena, as con- 
firmed by the data in Fig. 5, curve 0. 
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(b) Plastic Deformation 

It has been shown (Granato and Liicke 1956) that significant changes in 
internal friction at megacycle frequencies occurring in non-ferromagnetic 
materials strained within the elastic limit are associated with dislocation 
characteristics. On the basis of the same mechanism an increase in the ultra- 
sonic attenuation would be expected in specimens undergoing plastic deforma- 
tion for the reasons of increased dislocation density. An eventual decrease 
might follow at larger deformations where, due to high dislocation density, a 
significant interaction between the dislocations can be expected (Hikata et al. 
1956). 

In a ferromagnetic, however, the changes in attenuation in the elastic strain 
region appear to be determined mainly by magnetic phenomena, as shown 
before. In the plastic range, however, both the increase in the number of dis- 
locations and the influence of this increase on the ferromagnetic properties 
should be taken into consideration. The increase in dislocation density will 
cause an increase in the internal friction part associated with dislocation 
mechanisms. The mobility of the domain wall, however, will be reduced due 
to the effect of the increase in dislocation density. The contribution of the 
domain mechanisms to the internal loss is, therefore, expected to decrease, as 
the material is plastically deformed progressively, so that the attenuation 
related to this source should decrease also. 

The observed value of attenuation might be considered as the sum of the 
two above-mentioned effects. The ‘‘magnetic’’ part seems to be predominant 
in elastic straining. The “‘dislocation”’ part appears to be active in the behavior 
of attenuation at the yield stress region, in the sometimes erratic pattern of 
attenuation during plastic deformation, and in the increase in the initial 
attenuation value (at zero load) in heavily deformed ferromagnetics (not con- 
sidered in this paper, but to be published separately) as compared with 
annealed material. 

In the actual experiments, straining beyond the elastic limit sometimes did 
not seem to change the attenuation value at all (specimen 1, Figs. 2(6) and 
2(c)), but in other cases caused appreciable variations (specimen 2, Fig. 4). 
Specimen No. 1, however, was strained in steps, and it is reasonable to expect 
a certain degree of work hardening; this would render the subsequent deforma- 
tion less effective. The difference in the behavior of the different specimens is, 
however, still remarkable. 

The data on specimen 1 (Figs. 2(0) and 2(c)) suggests that the resulting 
plastic deformation was unable to affect the properties of the material as far 
as these are revealed by ultrasonic methods. The character of the attenuation 
changes on unloading seems to be in accordance with this and can be ex- 
plained by the same domain mechanism employed in discussing the curves in 
Fig. 2(a). The initial attenuation value (at zero load) remained nearly un- 
changed, as expected, considering the insignificant changes in the attenuation 
during the plastic deformation. 

A radical increase in strain (as used on specimen 2) was always found to 
produce substantial changes in attenuation, as indicated by Fig. 3. But similar 
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or even larger effects were also obtained at smaller strains (specimen 2, Fig. 
4). Relatively large fluctuations in the attenuation observed during plastic 
deformation very likely correspond to violent changes occurring within the 
metal. More consistent variations might be detected employing high-speed 
photographic technique. As before, the trend in attenuation changes during 
plastic deformation seems to determine the behavior of strained material. 
The unloading curves assume in both cases (Figs. 3 and 4) a characteristic 
shape. As the load is gradually removed, the attenuation starts to rise at first 
and after passing through a maximum, it decreases to a low value at zero 
load. The difference in the values at zero load before and after the stress cycle 
might not be convincing for the case in Fig. 3, but little doubt exists about the 
difference in Fig. 4. This drop in ultrasonic attenuation is very likely related 
to the reduced mobility of domain walls due to a higher dislocation density. 
The contribution of the “dislocation”’ part of attenuation is apparently less 
significant under these conditions than the effect of dislocation density on the 


domain wall mobility. 


(c) Effective Internal Stress 

A characteristic shape of the attenuation-strain curve obtained in sufh- 
ciently deformed specimens (the unloading curves in Figs. 3 and 4) was ob- 
served also in annealed specimens at very low strains (Figs. 4 and 6). Indica- 
tions of a similar trend might be found in Figs. 2(a) and 2(c), at the very 
beginning. In all cases, the characteristic feature is an initial increase in the 
attenuation as the load is increased. 

In aluminum a similar increase was noted at low strains and it was ex- 
plained in terms of dislocation loop mechanism (Hikata et a/. 1956). In Armco 
iron, however, this kind of mechanism does not seem to be probable at the 
low stresses employed. The dislocation loops are believed to be more effectively 
pinned by interstitial atoms (cf. the relatively weak pinning produced by 
substitutional atoms in annealed aluminum). 

A simple explanation in terms of internal stresses and their effect on domain 
configuration is considered here; it appears to account for the observed attenua- 
tion characteristics. It is assumed that the specimen at zero load is not ideally 
stress-free but contains internal stresses. These internal stresses might be 
either the deformation stresses in plastically strained material, or the residual 
stresses in annealed specimens. A strong directionality in the internal stress 
(compression) is expected in the specimens which underwent plastic deforma- 
tion by tensile loading. Therefore, a pronounced order in the domain con- 
figuration is to be expected in these specimens at zero load. As pointed out 
before, a conclusion about the degree of this order in a specimen in annealed 
and deformed state cannot be made solely on the basis of initial attenuation 
(at zero load) data. The observed attenuation is the sum of a ‘‘dislocation’’ and 
a ‘“‘magnetic’’ component. It is, however, possible to obtain an indication of the 
initial configuration by studying the specimen under low external stresses. 
The external stress will counteract the pronounced order, produced by the 
internal stress; the resulting less ordered configuration will manifest itself by 
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higher losses. As the external stress is increased a point will be reached where 
the respective components of the internal and the external stress will be 
equal and opposite. This point should correspond to a maximum in attenua- 
tion. A further increase in the external stress will cause a decrease in losses by 
renewed domain ordering. 


ADDITIONAL EXPERIMENTAL WORK 


The various characteristics found previously on specimens 1 and 2 were also 
obtained on a further specimen, specimen No. 3, which was strained in two 
steps (Figs. 6 and 7). The observed initial attenuation increase (Fig. 6) seems 
to be associated with the residual internal stresses. The sharp drop afterwards 
is caused by the renewed domain ordering discussed above and the similarly 
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Fic. 6. Ultrasonic attenuation versus strain for specimen No. 3. The data correspond 
to the first stress cycle. Stress-strain curve included. 


sharp rise which follows is related to an efficient plastic deformation. The 
character of the unloading curve most certainly indicates internal stresses in 
the deformed material, the attenuation falling continuously as the load is de- 
creased. The value at zero load is lower if compared with the initial attenua- 
tion in the annealed material. 

The second loading (Fig. 7) after a rest of 15 minutes at zero load followed 
in reverse the previous unloading curve very closely. With the external and 
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internal stresses equal, a maximum in the attenuation value is observed. A 
sharp decline which follows appears to be related to the ordering effect of ex- 
ternal stress on the magnetic domains. The plastic deformation taking place 
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Fic. 7. Ultrasonic attenuation versus strain for the second stress cycle of specimen No. 3. 
Stress-strain curve also shown. 


in this stress cycle shows, however, a surprisingly small effect as far as attenua- 
tion is concerned. A similarity with Fig. 2 exists in this respect, and the 
character of the curve obtained on unloading would actually correspond to a 
specimen with low effective internal stresses. No satisfactory explanation for 
this behavior can be offered at present. An increase in the attenuation at zero 
load might be noted also. 


SUMMARY 


An ultrasonic technique employing megacycle frequencies was used in order 
to study the changes in Armco iron on straining. 

The observed changes seem to be accounted for by magnetic domain theory 
and dislocation mechanisms. Three features which appear to be characteristic 
were more closely studied. 

1. In the region where the material behaves elastically, an increase in the 
external stress is accompanied by a decrease in ultrasonic attenuation. The 
effect was shown to be related to the rearrangement of magnetic domains by 
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the externally applied stress. The application of an external magnetic field 
seems to confirm this explanation. 

2. An initial increase in attenuation occurring in the elastic region with some 
annealed and with most strain-hardened specimens appears to be associated 
with internal stresses, i.e. residual stresses left after annealing or stresses intro- 
duced into the material by plastic deformation. These are relieved upon the 
application of an external tensile stress and the result is a magnetic domain 
configuration exhibiting higher internal loss. 

3. The changes in ultrasonic attenuation accompanying plastic deformation 
varied widely among individual specimens—large changes occurred in some 
cases and small changes in others. The difference appears to be related to the 
difference in straining history. The degree of strain hardening as registered by 
ultrasonic attenuation seems to correspond to the behavior manifested by the 
material during the plastic deformation. 
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THE IMPEDANCE, RECTIFICATION, AND ELECTRO- 
LUMINESCENCE OF ANODIC OXIDE FILMS ON ALUMINU}: 


ARCHIBALD W, SMITH 


ABSTRACT 


An experimental investigation was made of the impedance, rectification, and 
electroluminescence of anodic oxide films formed on aluminum, in an electrolyte 
of boric acid and ammonium hydroxide. These properties were studied, where 
appropriate, as a function of the surface condition before anodization, the for- 
mation voltage, the conductivity of the electrolyte, and the frequency. Aconstant- 
voltage formation technique was used. The interpretation of the results suggests 
that the oxide itself is a good dielectric, but that the film contains defects whose 
number and form depend on the conditions of formation. The defects appear to 
control the rectification and electroluminescence and have an appreciable in- 
fluence on the impedance. The frequency dependence of the capacitance was 
much smaller than that found by previous workers although there was an ap- 
preciable contribution arising from geometrical effects in the electrolytic cell. 
The zero-thickness potential for the film was found to be in agreement with a 
theoretical estimate. 


1. INTRODUCTION 


Anodic oxide films on aluminum play an important technological role, 
especially in the construction of capacitors, and have received considerable 
study in the last 60 years. Although many features of their structure and 
properties are now understood, there is still a need for further fundamental 
investigation. The present paper is concerned with the impedance, rectification, 
and electroluminescence of films formed in an electrolyte which does not dissolve 
the oxide. At the present time, the origins of the rectification and electro- 
luminescence have not been established, although a number of theories 
exist. There are many features of the impedance which require further eluci- 
dation, particularly the large frequency dependence of the capacitance found 
by previous workers, and the behavior at low formation voltages. In the work 
described here measurements of these properties have been made as a function 
of the formation voltage for two surface conditions, one rough, the other 
smooth. The impedance was further investigated as a function of the frequency 
and the conductivity of the electrolyte. The interpretation of the results 
suggests that the oxide itself behaves as a good dielectric, but that the con- 
tinuity of the film is disturbed by defects which have an appreciable influence 
on film properties. This is in agreement with recent work on oxide films on 
tantalum and zirconium. 


2. EXPERIMENTAL 

2.1 Preparation of Samples for Anodization. Electrolytes 

The experiments were carried out on aluminum foil of minimum purity 
99.99%. Square samples of nominal area 10 cm? were cut from the foil, which 

1Manuscript received January 29, 1959. a 
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was 7.5X10- cm thick, leaving a tab 0.16 cm wide and 1.5 cm long projecting 
from the center of one side. The tab was used to make electrical connection 
to the sample and to suspend it in the electrolyte. In order to define the area 
anodized, and to eliminate undesirable reactions at the air—metal-solution 
boundary, the tab was coated with wax. The surface of the aluminum was 
prepared for anodization in two ways: (1) washed for 2 minutes in a 5% 
solution of sodium hydroxide at 35° C; (2) etched for 45 seconds in a mixture 
of 45% sulphuric acid (by weight), 43% phosphoric acid, and 12% water at 
100° C, and electropolished (Tegart 1956) in the same mixture for 4 minutes 
at a current density of 100 ma/cm?. The first method gave a superficial 
cleaning only, leaving the surface contours of the aluminum essentially 
unchanged. The second method removed an appreciable amount of aluminum 
and left the surface microscopically smooth. Samples treated by these methods 
will be referred to as ‘washed’ and ‘electropolished’, respectively. Photo- 
micrographs of the washed and electropolished surfaces are shown in Fig. 1. 

The electrolyte used for anodization was a 5% solution of boric acid with 
sufficient ammonium hydroxide added to obtain the desired conductivity. 
Solutions with conductivities between 0.8X10-* and 2.5107? mho/cm 
were used, corresponding to pH values between 5.1 and 5.8. The chemicals were 
of reagent grade, with a maximum stated chloride content of 2 parts in 10°. 
Water purified with an ion-exchange resin was used throughout for making 
up solutions and washing samples. 


2.2 The Electrolytic Cell. Formation Technique 

The formation of the oxide films and the measurements, with the exception 
of the electroluminescent brightness (Section 2.4), were carried out in a cell 
fitted with a platinized-platinum counterelectrode (the cathode during for- 
mation), a calomel reference electrode, and a stirrer. The platinum electrode 
took the form of a rectangular box (13 in. long, 1/2 in. wide, and 1} in. high) 
with the top and bottom open. The aluminum samples were suspended in the 
center of this box, parallel to the long side. The reference electrode was intro- 
duced by means of a liquid bridge filled with the same electrolyte as the cell. 
A thin tube (0.08 in. outer diameter) on the end of the bridge passed through 
the platinum electrode and ended about 0.05 inch from the sample. In order 
to prevent the diffusion of chloride ions from the calomel electrode into the 
cell, asbestos plugs were inserted into the tip of the calomel electrode and into 
the bridge. The temperature of the cell was not controlled but all work was 
done in the range (26+2)° C. The cell was open to the atmosphere. 

The measured capacitance of the cell includes the capacitance of the 
platinum electrode in series with that of the oxide film. The capacitance of 
this electrode was determined to be about 3100 uf at 2 kc/s and 1200 uf at 
20 kc/s in the solution of lowest conductivity, and larger in the other solutions. 
The maximum measured capacitance was 6.2 uf at 2 kc/s and 4.5 uf at 20 
kc/s. Thus the error was only 0.2 or 0.3% in the most unfavorable case. The 
change in series resistance of the platinum electrode between 0.1 and 20 kc/s 
was determined to be less than 0.1 ohm. The absolute value of this resistance 
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was of little significance because the body of the electrolyte added a relatively 
large, unknown series resistance. 

The formation of the oxide films was carried out by applying a constant 
voltage to the cell through a series resistor. This resistor limited the initial 
current density to about 20 ma/cm?, and was reduced as formation proceeded 
in such a way as to keep the current density equal to the above value for as 
long as possible. Below a formation voltage (V;) of about 15 volts, the resis- 
tance of the cell itself resulted in an initial current density of less than 20 
ma/cm?. The time of formation was increased gradually from 15 minutes 
at low voltages (V; < 10) to 60 minutes at high voltages (V; > 200). The final 
current density varied between | and 50 ywa/cm?’, depending on the surface 
treatment and the applied voltage, as shown in Table I. Just before stopping 
formation, the voltage between the platinum and reference electrodes was 
read on a vacuum tube voltmeter. The sum of this voltage and the applied 
voltage gave the total formation voltage with respect to the reference electrode. 


















TABLE I 
Average cell current at the end of film formation 


















Formation voltage, Vr 














Surface —_—— —_———— —— —— — —--——. 
treatment 10 20 50 100 200 300 400 
Electropolished 2.0 1.8 1.3 2.5 4.3 6.0 6.4 
2.2 3.0 2.2 3.9 5.1 7.6 35.0 











Washed 





Note: Current is in units of 10-* amp/cm?*. 











2.3 Impedance Bridge 

The impedance of the cell described above was measured on a series resis- 
tance-capacitance bridge designed to measure capacitance from 0.1 to 10 uf, 
and resistance from 10 to 1000 ohms, over the frequency range 0.1 to 20 kc/s. 
Provision was made for applying d-c. potentials up to 200 volts to the cell 
during measurements. The bridge was equipped with a Wagner ground circuit 
to balance stray capacitance to ground. The maximum a-c. signal applied to 
the cell was 50 millivolts. At frequencies of 5 kc/s and below, the measured 
capacitance was accurate to +0.2% or better. Above 5 kc/s, the errors were 
larger than this, depending on the magnitude of the capacitance, and the 
resistance in series with it. For the measurements reported here, the upper 
limit of error was about +2% at 20 kc/s. The resistance values were accurate 
to +1% or better at all frequencies. 



















2.4. Electroluminescence Measurements 

The electroluminescence of the oxide film was detected with a type 931-A 
photomultiplier. The electrolytic cell was essentially the same as in the 
impedance measurements, except that the reference electrode was omitted, 
and the platinum counterelectrode expanded in size to 33 in. by 3 in. by 1/2 in. 
The samples were correspondingly increased in size to 3 in. square, in order to 
obtain a higher photocurrent. One of the large sides of the counterelectrode 
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was a screen of platinum wire, consisting of 0.010-in. diameter wires placed 
0.25 inch apart. The light output was observed through this screen and a 
plane window in the cell container. The photocathode-to-sample distance was 
about 3 inches. With the particular photomultiplier tube used, the smallest 
detectable photocurrent was judged to be about 1.0X107-!° amp. This figure 
was determined chiefly by the noise in the photomultiplier. 


3. RESULTS 
3.1 Impedance 

The a-c. impedance measurements were carried out in order to investigate 
the relationship of the equivalent series capacitance (C,) and resistance (R,) 
of the oxide film to the formation voltage (V;) and frequency of measurement 
(f). It was also desired to know whether the surface treatment before ano- 
dization and the conductivity of the electrolyte had any significant effect. 
To achieve these aims, measurements were performed on four groups of 
samples, which covered preparation by washing and by electropolishing, and 
anodization in high- and low-conductivity solutions. Each sample in a group was 
anodized at a series of successively increasing V;, and at each V;, C, and R, 
were measured at 0.1, 0.2, 0.5, 1, 2, 5, 10, and 20 kc/s. The series of V; for 
each sample usually extended over the range 10 to 400 volts. The measure- 
ments were made about 10 minutes after ceasing formation, because the 
capacitance dropped about 1% in the first few minutes. 

Some measurements were also made on samples anodized over the range 
2 to 10 volts. The general program of the measurements was the same as 
above. The area of the samples was reduced by one-half to 5 cm? for these 
measurements, the readings at 0.1 and 20 kc/s were omitted, and the con- 
ductivities of the solutions increased somewhat. The measurements were 
again made 10 minutes after ceasing formation. 

The capacitance results will be considered first. Apart from certain exceptions 
noted below, the reproducibility between samples varied from 0.1 to 4%, 
but was usually about 1%. The relationship between 1/C, and V;, (at a fixed 
frequency) was nearly linear, except for washed samples above V; = 200. 
The slopes and intercepts on the 1/C, axis of the straight lines fitting the data 
both increased slightly with frequency. These results are illustrated for electro- 
polished samples (2.5 X 10-* mho/cm electrolyte) in Fig. 2. The situation near 
the 1/C, axis is shown in more detail in the lower insert of this figure. The 
formation voltage was referred to the potential of the saturated calomel elec- 
trode in all the data. The results for washed and electropolished samples were 
much the same below V; = 200. The behavior of washed samples at V; = 300 
and 400 was not very reproducible, but typical measurements are shown by 
the dashed lines in Fig. 2. It will be seen that they deviate considerably from 
the straight line established at lower voltages. Decreasing the conductivity 
of the solution caused the intercepts on the 1/C, axis to increase more rapidly 
with frequency, and reduced the deviations from the straight line for washed 
samples above V,; = 200. Straight lines have been fitted to the data at 0.2 
and 10 ke/s by a least squares procedure, and the resulting constants for the 
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PLATINUM ELECTRODE 


ALUMINUM 


In cm*/uF 


(CAPACITANCE )~! 


(CAPACITANCE) ww cM?/uF 


z 
FORMATION VOLTAGE (REFERRED TO SATURATED CALOMEL ELECTRODE) 
Fic. 2. Reciprocal capacitance plotted against formation voltage. Solid lines: electro- 
polished samples (2.5 X10~* mho/cm electrolyte). Dashed lines: washed samples at Vs = 200, 


300, and 400 (1.6X107-? mho/cm electrolyte). Upper insert: equipotential (solid) and field 
(dashed) lines for one quadrant of horizontal plane through center of electrode system. 


equation 1/C, = a+6V; are shown in Table II. By comparing the various 
cases, it will be seen that extrapolation of high formation voltage data does 
not always yield accurate values of the intercept a. 


TABLE II 


Constants of the equation 1/C, = a+6V; 


Frequency 


10 kc/s 
Surface Range of Conductivity 
treatment Ve of electrolyte b 


Electropolished 10 1.61078 175 
2. 169 

400 0.83 171 
2.: 169 

2-10 2.5 164 
10-200 0.83 41 167 
1.6 37 0.167 


Washed 


Note: Units—e in cm?/yf, 6 in cm*/yf-volt, V¢ in volts, conductivity in mho/cm. 
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The equivalent series resistance results will be discussed next. In general, 
these were less reproducible than the capacitance results because the resistance 
of the electrolyte (R.) appeared to vary slightly between anodizations and 
between samples. The results for electropolished samples in 2.5 X 10-* mho/cm 
electrolyte are plotted against the reciprocal of the frequency (1/f) in Fig. 3. It 
will be seen that at frequencies below 2 kc/s, R, was proportional to 1/f. At 


IN OHMS 


SERIES RESISTANCE OF CELL 





2 4 6 8 10 
( FREQUENCY )~! IN MILLISECONDS 


Fic. 3. Series resistance of the electrolytic cell plotted against reciprocal frequency for 
electropolished samples in 2.5107? mho/cm electrolyte. 


frequencies above 2 kc/s, there was an appreciable curvature, which made it 
impossible to determine R, by extrapolation to 1/f = 0. This curvature in- 
creased as the conductivity of the electrolyte was decreased. For V; less than 
about 10 volts, R, was found to be independent of V;, except perhaps at 0.1 
kc/s. The results for washed samples were very similar except that the slope 
of the linear portion for V; = 400 was about three times as great. 

Three special runs were made to determine the change in impedance when 
the formation voltage was reapplied after having ceased formation for a short 
time. The procedure was as follows: Films were formed to successively higher 
voltages in the same manner as before. At each formation voltage, the im- 
pedance was measured at 0.2, 2, and 20 kc/s, 20 minutes after ceasing for- 
mation. The capacitance and resistance values so obtained have been desig- 
nated C, and Rj, respectively. The formation voltage was then reapplied and 
10 minutes later, the impedance measured at the same frequencies, yielding 
C2 and R». Finally, the voltage was shut off again and 20 minutes later, the 
impedance remeasured, yielding C; and R3. The quantities AC; = (C2—C1)/Co, 
AC; = (C2—C3)/Co, and AR; = (R2—R;)/(R2—R,) are tabulated in Table ITI. 
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The entries in the table are average values of the quantities over the range of 
V; shown in column 2. There was no significant difference between R; and R;, 
and AC did not show any significant trend with the change in thickness. It will 
be seen that the application of the formation voltage increased the capacitance 
by a small percentage, which decreased with increasing frequency. The resis- 
tance was increased by a much larger factor, which also decreased with 


increasing frequency. 


TABLE III 


The change in impedance upon reapplication of the formation voltage 


Frequency in kc/s 
Range bere oie = 





Quantity of V¢ 0.2 2 20 
AC, 5-250 034 .025 .018 
AC; 5-250 .027 .018 O15 
AR 5-25 .27 — 

30-75 .69 eI 
100-250 1.05 .30 





Note: See text for definitions of ACi, AC2, and AR. 


3.2 Current—Voltage Characteristics 

The purpose of the measurements described here was to investigate the 
effect of the surface treatment before anodization on the d-c. current-voltage 
characteristics. To obtain these, a sample was successively anodized for 1 
hour at V; = 100, 200, 300, and 400 in a 1.6X10-* mho/cm electrolyte. At 
each V;, a current-voltage curve was taken by decreasing the applied voltage 
in steps, starting at the formation voltage, passing through zero current, and 
then increasing the voltage with the opposite polarity. Current readings were 
taken about 1 minute after the voltage was set, except at low anodic voltages, 
where about 5 minutes were allowed. The reproducibility between samples 
was rather poor, but certain general trends became evident, as illustrated 
in Fig. 4, for V; = 100 and 400. Before the results were plotted, the zero- 
current (i.e. open-circuit) voltage of the system (usually in the range +1 
volt) was subtracted from the applied voltage. 

It will be seen from Fig. 4 that on a log-log plot, each curve in the anodic 
or reverse direction consisted of a rapidly falling curved section just below 
the formation voltage, followed by a more or less linear section at lower 
voltages (cf. Charlesby 1953; Smith 1957). The experimental points deviated 
least from the straight line for the 300- and 400-volt cases. There was no 
unique relation between the current and the thickness at a given voltage 
along the linear section. However, the current was usually a minimum for 
V; = 100 or 200, and a maximum for V; = 300 or 400. For a given V;, the 
main variation between samples was a displacement of the linear section 
parallel to the current axis. For instance, at V; = 200, the current varied from 
2X 10-'° to 2X 10-8 amp/cm?’, with 10 volts applied. 

The cathodic or forward curves were steeply rising and, at higher voltages, 
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the current consistently decreased with increasing thickness. At lower voltages, 
the curves often crossed over. 
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Fic. 4. Current-voltage characteristics for washed and electropolished samples. 


It is evident that the washed and electropolished samples behaved in much 
the same way, but there were two pronounced differences. First, the anodic 
current for washed samples with V; = 400 was very erratic when the applied 
voltage was close to V;, and was greater by a factor of 5 to 10. Second, above 
5 volts, the cathodic current for washed samples was greater by a factor of 
about 100. 

The determination of the current-voltage curves was not entirely straight- 
forward due to drift of the current. The drift was investigated in a few samples 
by measuring the current at 1-, 2-, 5-, and 10-minute intervals after the voltage 
was set, with the following results: After a decrease in the reverse voltage 
in the range V; > V 2 0.8 V;, the current drifted downward by amounts up 
to 40% in 10 minutes. In the range 0.8 V; 2 V = 0.02 V;, the drift was less 
than 5% in 10 minutes. In the range 0.02 V; = V > 0, a decrease in voltage 
often caused a reversal in current, but the current always drifted back to 
the anodic direction, so that the reversal was temporary. A rough equilibrium 
was usually established within 10 minutes. This effect prevented the measure- 
ment of the zero-current voltage to better than 0.1 or 0.2 volt. In the forward 
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direction, there was little drift up to 3 or 4 volts. Beyond this, the current 
drifted rapidly upward after an increase in voltage, without appearing to 
reach an equilibrium. 


3.3 Electroluminescence 

A study was made of the dependence of the electroluminescence (abbreviated 
EL) on the surface treatment before anodization and on the formation voltage. 
At each formation voltage, both the photocurrent, which is proportional to the 
electroluminescent brightness, and the cell current were measured as a function 
of the applied voltage, in the same way as the cell current in the previous 
section. Observation of the EL in the forward direction is very difficult because 
of the streams of hydrogen gas coming off the surface, and no results for this 
case will be given here. The results in the reverse direction for washed and 
electropolished samples are shown in Fig. 5. The ratio of the photocurrent 
to the cell current (J,/J.) has been plotted against the field in the film for 
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Fic. 5. Ratio of the electroluminescent brightness to the electrolytic cell current plotted 
against the approximate field in the film. The curves are labelled with the formation voltages. 


V, = 100, 200, 300, and 400. The field was computed approximately by as- 
suming that the thickness d of a film formed to a voltage V; was d = 1.31077 
V; cm. For washed samples at high fields, both the photocurrent and cell 
current were rather erratic, and the larger fluctuations in these currents were 
correlated. Under these conditions 7,/J, was nearly constant and therefore 
I, was approximately proportional to /,. This is in agreement with the results 
of van Geel, Pistorius, and Bouma (1957) for EL at the formation voltage with 
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a constant film thickness, although their measurements were made in a some- 
what different way. At low fields J, fell much more rapidly than /,. The 
ratio J,/J, was much smaller for the electropolished samples, and as 7, was 
about the same (cf. Fig. 4), 7, itself was also much smaller. 


4. DISCUSSION AND INTERPRETATION OF RESULTS 
4.1 Impedance 

In this section, an attempt is made to relate the: measurements of the 
equivalent series capacitance and resistance to phenomena within the electro- 
lytic cell. A simple circuit representation of the cell is first proposed on the 
basis of a few simple assumptions. By comparing the predicted and observed 
values of the impedance, it is shown that this simple circuit is inadequate 
due to certain geometrical effects in the cell. Taking these into account, the 
section is concluded with a discussion of the relation of the impedance measure- 
ments to the properties of the oxide film itself. 

The simplest circuit representation of an electrolytic cell containing an 
anodized electrode consists of a capacitor, representing the effect of the 
oxide film, in series with a resistor equal to the resistance of the electrolyte. 
(It is assumed that the impedance of the counterelectrode has been made 
negligible.) If the oxide has a thickness d and a complex dielectric constant 
e* = e(1+7D), its capacitance per unit area will be given by C; = e¢/d, in 
rationalized units. The loss tangent D for the oxide turns out to be sufficiently 
small that C; may be equated directly to the equivalent series capacitance 
of the circuit C, = C; = e/d. The equivalent series resistance of the oxide 
is given by R; = D/wC,; and thus the series resistance of the circuit is 
R, = R.+D/wC;. If the oxide behaves as a good dielectric, « and D will vary 
only slightly with frequency, generally decreasing as the frequency is increased 
from 0.1 to 20 kce/s (cf. von Hippel 1954). The linear relation which holds 
between 1/C, and V;, at least over limited ranges of V;, implies that d = kV, 
where & is a constant. 

This circuit gives a good representation of the frequency and thickness 
dependence of the impedance of cells containing tantalum and niobium (cf. 
Young 1955). However, it is not adequate for the present results for aluminum, 
because extrapolation of the data of Fig. 2 to 1/C, = 0 would give a zero- 
thickness voltage (Vo) which was frequency dependent. This is, of course, 
impossible. The curvature in the R, vs. 1/f plots (Fig. 3) at frequencies below 
2 kc/s is also inconsistent with the above circuit. 

These discrepancies can be explained by taking into account the frequency 
dependence arising from the geometry of the electrolytic cell. Suppose that the 
electrolyte joining the two electrodes is divided into conducting paths of 
equal resistance. For alternating currents, the electrical circuit from each 
path to the anodized electrode is completed through a capacitance due to the 
oxide film. Thus the cell may be effectively represented by a number of 
circuits in parallel, each consisting of a resistance and capacitance in series. 
If any of the capacitances are different, the resistances all being taken equal 
by definition, the impedance of the network will have a different frequency 
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dependence from a single resistance and capacitance in series. Whether or 
not the capacitances are equal is determined by the areas on which the con- 
ducting paths terminate at the anodized electrode, assuming that the film 
has a perfectly uniform structure and thickness. The areas are determined 
by the solution of Laplace’s equation for the electrode geometry of the par- 
ticular cell. The frequency dependence increases with increasing series re- 
sistances. Thus it will be larger with the low-conductivity solutions which 
must be used with aluminum than with the high-conductivity solutions nor- 
mally used with tantalum and niobium. 

Equipotentials and field (or current) lines have been obtained for the 
horizontal plane through the center of the present electrode system by a 
graphical solution of Laplace’s equation. These are shown for one quadrant 
in the upper insert of Fig. 2. The regions between the field lines represent 
conducting paths of equal resistance. It will be seen that paths near the end 
of the aluminum terminate on smaller areas, and hence on smaller capacitances, 
than those near the center. The situation is very similar even if the end of the 
platinum electrode is not closed. Therefore the resultant impedance of the 
complete electrode may be evaluated approximately by considering only the 
quadrant shown in Fig. 2. If the capacitance associated with the mth path 
is C,, and the resistance of each path is R, the equivalent series capacitance 
and resistance of the system are given by: 


l wB G 
C, = G+ y? and R, = Gt Vy" ’ 
where 
1/wC, _ ; — ee 
B = X R+( aC.) and G => » R'+(1/eC,)" : 


There are two limiting cases: (1) At low frequency and/or capacitance 
wR «K1, C,= LC,, and R, = RYC?/(LC,)?. This approximation is 
adequate at 0.1 and 0.2 kc/s under the conditions of the present experiments, 
except perhaps for V; <4. (2) At high frequency and/or capacitance 
wC,R > 1, 1/C, = (1/n?)>-(1/C,), and R, = R/n. This approximation is 
not valid at all under the present conditions. Calculations of 1/C, and R, 
were made at a number of frequencies and formation voltages, using the upper 
insert of Fig. 2, and observed C, at 0.1 kc/s to determine the capacitances 
C,. R was estimated from observed R, at 20 kc/s. The details of these calcu- 
lations will not be given here, but it was found that the observed frequency 
dependence of the intercepts on the V; axis, and of the R, values for V; < 10, 
was approximately reproduced. The 1/C, vs. V; plots at the higher frequencies 
were curved, with the second derivative d?(1/C,)/dV, negative and in- 
creasing in magnitude with frequency. For instance, the slope at 10 kc/s 
was 12% larger than that at 0.1 kc/s over the range V; = 4 to 10, but was 1.2% 
smaller over the range V; = 50 to 100. This curvature explains the extra- 
ordinarily high slopes obtained at 10 kc/s for the low formation voltage cases 


of Table II. 
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The relation of the impedance measurements to actual properties of the 
oxide film can now be discussed. First, the question of the zero-thickness 
voltage Vo will be reconsidered. Experimentally, it is given most accurately 
by extrapolation of the low formation voltage data at 0.2 kc/s to 1/C, = 0. 
An average value of Vo = (—2.4+0.1) volts is obtained from the three cases 


in Table II for which 2 < V; < 10. The frequency dependence arising from 
the cell geometry will cause this value to be about 0.1 volt too low algebraically, 
and thus the best experimental value is Vo = (—2.3+0.1) volts. If certain 
simplifying assumptions are made, Vo» may be calculated theoretically 
(Vermilyea 19540; Young 1957a). When oxide-free aluminum is placed in the 
electrolyte, oxide immediately grows until d = kg, where ¢ is the relative 
oxidation potential associated with the reaction for the production of oxide: 


2Al + 3H.0 = Al.O; + 6H* + 6e. 


If the aluminum is anodically biased (by introducing a counterelectrode and 
an external voltage), more oxide grows on top of the initial layer. Using the 
free-energy data for bulk Al.O3, the value of @ may be estimated as +1.50 
volts, under standard conditions (Latimer 1952). The electrolytes used in this 
work had a hydrogen-ion concentration of about pH = 5.2. When this is taken 
into account, and potentials are referred to the saturated calomel electrode, 
the’value — 2.05 volts is found for Vo, in fair agreement with the experimental 
value of (—2.3+0.1) volts. 

The slope 6 of the 1/C, vs. V; plots will now be briefly discussed. The relation 
b = k/e may be easily derived on the basis of the simple circuit representation. 
The slopes in Table II may be compared with previous measurements of k 
and e. These quantities have been measured by many workers under various 
conditions; k = 1310-8 cm/volt (Dekker and van Geel 1947; Hass 1949) 
and ¢/eo = 8.7 (van Geel and Schelen 1957) may be taken as representative 
values. The value of b derived from these is 0.17 cm/yuf-volt, which is in 
satisfactory agreement with the values in Table II. The increase in } with 
frequency for the data in this table means that ¢ decreases slightly, as ex- 
pected. 

The loss tangent D for thick films (V; > 100) has been evaluated using 
the formula D = wC,R,, with C; equal to C, for 0.1 kc/s and R, determined 
approximately as R,(V; > 100)—R,(V: < 10). Subtraction of the resistance 
R,(V; < 10) for thin films corrects for the resistance of the electrolyte, and 
approximately for the frequency dependence arising from the cell geometry. 
The values obtained are shown in Table IV for V; = 200 and 400. The 
values for V; = 100 and 300 volts were much the same as for V; = 200. 
At high frequencies R, is the small difference of large numbers so the un- 
certainty in D is quite large. Therefore, D is given at low frequencies only 
in Table IV. For each surface treatment, D tends to be lower in the solution 
of lower conductivity. For washed samples at V; = 400, in the 2.51073 
mho/cm electrolyte, D is about three times greater than average, and this 
fact appears to be related to the higher currents at the end of formation for 
the same conditions (Table I). 
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TABLE IV 


Frequency in kc/s 











Surface Formation Conductivity ———_ a ee 
treatment voltage, Vr of electrolyte 0.1 0.2 0.5 1.0 
Electropolished 200 0.83 X 107% 6.8 7.9 6.3 6.1 
2.5 9.6 9.1 8.0 6.9 
400 0.83 10.0 12.6 8.7 a 
2.5 12.6 12.0 10.6 10.1 
Washed 200 ; 9.5 9.7 8.2 8.9 
1.6 9.9 9.7 8.1 4.2 
400 0.83 12.5 10.8 11.3 13.0 
1.6 32 31 30 30 


Note: Entries give the loss tangent in units of 10-3. Conductivity is in units of mho/cm. 


The increase in capacitance (or e) which results from the application of the 
d-c. formation voltage (Table III) appears to be related to effects in the body 
of the oxide film, because percentagewise it is independent of thickness. The 
corresponding percentage increase in resistance increases with thickness. This 
implies an increase in D which tends to be independent of thickness. These 
changes are probably due to the introduction of extra carriers through the 
action of the applied field. The decrease in the change of ¢ and D with in- 
creasing frequency suggests that these carriers have a relaxation time of 107° 
or 10-4 seconds, and therefore are ions. The extra ions could be injected into 
the oxide itself or into defects in the film opened up by the d-c. field. 

The impedance of electrolytic cells containing anodized aluminum electrodes 
has previously been measured as a function of formation voltage and frequency 
by van Geel and Scholte (1951), Scholte and van Geel (1953), and Lilienfeld 
and Miller (1953). Although valid comparisons are difficult to make because 
of the differing conditions, their results are qualitatively similar to those 
reported here. One point worth discussing in some detail is the frequency 
dependence of the capacitance. The situation is illustrated in Fig. 6, in which 
the reciprocal capacitance, normalized at 1 kc/s, is plotted against frequency 
on a log scale. It will be seen that the frequency dependence of the present 
results is considerably smaller, even if the extra dependence caused by the 
geometry of the cell and by the application of a d-c. field is included. Neither 
pair of authors considered explicitly the possibility of frequency-dependent 
effects due to the geometry of their cells, although the conductivity of their 
electrolytes was low enough to make such effects important. However, such 
effects alone cannot account for the difference in the frequency dependence, 
and it appears that actual differences in the structure of the film are involved. 
An important factor here may be the number and form of defects in films 
produced under different conditions (see Section 4.2). 

Both pairs of authors found non-zero intercepts on the 1/C, axis and in- 
terpreted these as evidence of interfacial effects between the oxide and the 
electrolyte. However, in neither case was a reference electrode employed, nor 
was any attempt made to calculate the theoretical voltage at which the 
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Fic. 6. Reciprocal capacitance, normalized at 1 kc/s, plotted against frequency, on a log 
scale. Curve A: results of this paper for washed samples with V;_ = 100, and 1.61073 mho/cm 
electrolyte. Curve B: same as curve A, but corrected approximately for frequency dependence 
of the cell. Curve C: results of Scholte and van Geel (1953, Table V) for V; = 100. Curve D: 
results of Lilienfeld and Miller (1953, Fig. 5) for V¢ = 180. Curve E: typical frequency de- 
pendence of a niobium film (Young 1955, Table I). 


thickness would be zero. Lilienfeld and Miller used an aluminum cathode and 
extrapolated to V; = 0 using data over the range V; = 100 to 400. It is 
unlikely that meaningful intercepts can be obtained under these conditions. 
Unfortunately, van Geel and Scholte analyzed their results by a method 
which is sensitive to systematic errors in the high-frequency measurements 
(such as would be caused by geometrical effects in the cell). For these reasons, 
the evidence for interfacial effects appears to be of questionable validity. 


4.2 The Role of Defects in the Oxide Film 

Oxide films formed on tantalum and zirconium appear to have defects whose 
number and size depend on the roughness of the surface and on the thickness 
of the film (Vermilyea 1954a; Young 19575 and 1958). The defects have an 
observable effect on the capacitance and loss tangent, and on the rectification. 
Young (19576) has suggested that they arise from compressional stresses build- 
ing up due to preferential oxidation of sharp projections of metal, and that they 
might take the form of cracks filled with electrolyte or gas. However, many 
of the defects would probably take less violent forms, ranging in size down 
to minute pores containing only a few water molecules. Impurities in the 
metal might cause similar defects by interfering with the normal formation 
of oxide. 

The present experimental results for aluminum are consistent with the 
concept of defects in the film. This is brought out most clearly by the difference 
in the capacitance (Fig. 2), the loss tangent (Table IV), and the current 
at the end of formation (Table I) for washed and electropolished samples. It 
is evident from this data that the defects become more serious as the for- 
mation voltage and the conductivity of the electrolyte are increased. If it is 
assumed that most of the d-c. current passes through defects, rather than 
the oxide, two aspects of the current-voltage characteristics of Fig. 4 can 





: 


Sa neaLH UALS NETETTON 


Se RS TEP 





\ 
t 
i 
i 
i 


SMITH: ANODIC OXIDE FILMS 605 


be more easily understood, namely the increase in the reverse current with 
increasing thickness, and the relatively greater forward current for washed 
samples. Young (1958) has concluded that the observed current-voltage 
characteristics for the tantalum are chiefly determined by defects. The present 
results for aluminum appear to be more consistent with this view than with a 
previous theory (Smith 1957) in which rectification arose from the difference 
in the behavior of holes and electrons at the oxide/electrolyte interface. 
However, the mechanism by which the defects rectify still remains to be 
determined. 

The influence of the defects on film properties will be more extensive if it 
is possible for electronic breakdown to occur at a defect. It is well known that 
penetration of water vapor into pores lowers the breakdown strength of 
insulators (e.g. Inuishi and Powers 1957). The extra conduction due to the 
water probably leads to a local heating which causes breakdown in the sur- 
rounding insulator. A similar type of breakdown would be expected to occur 
at pores in anodic oxide films, but in this case the breakdown would be partially 
blocked off by the formation of new oxide and/or a gas bubble. If new oxide 
were formed, it would probably be thicker than normal because of the local 
temperature rise. By this mechanism, minute pores having little or no effect 
at low voltages could give rise to more serious defects and changes in film 
properties when the formation voltage was raised. Other, though less probable, 
types of breakdown are also possible. For instance, the intrinsic “high tempera- 
ture”’ type of breakdown (Frohlich and Simpson 1950) could occur in regions 
of oxide containing a large number of electron traps. 

Two further experimental observations can be more easily understood on 
the basis of electronic breakdown. First, the erratic nature of the cell and 
photocurrents for washed samples (Section 3.3) could be due to the random 
occurrence of breakdown. Second, the conductivity of the electrolyte could 
influence the film properties through its effect on breakdown. This may be seen 
as follows: The overvoltage necessary to drive a given current through a pore 
will be a function of the conductivity of the electrolyte. In general, it would 
be expected that the overvoltage would be higher for a solution of lower 
conductivity. This was found to be true experimentally at a platinum electrode 
in the electrolytes used to anodize the aluminum in the present work. If this 
is true for the oxide electrode also, it would be expected that a higher applied 
voltage would be required to cause breakdown with the solution of lower 
conductivity. This is consistent with the observed effect of the conductivity 
on the impedance. 


4.3 Electroluminescence 

There are two principal explanations for the steady EL in the anodic 
direction: (1) Collisions of electrons with impurity atoms which act as 
luminescent activators. This mechanism has been discussed previously (Smith 
1957; van Geel et a/. 1957). (2) Breakdown effects in the oxide. The emission 
of visible radiation would be expected from a region of breakdown. The well- 
known scintillations which become visible at high formation voltages are a 
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striking example of this. When these occur, the damage to the film is extensive 
enough to be seen under the microscope (Vermilyea 1954). 

The relative importance of the first and second mechanisms will depend 
on conditions in the particular film under study. For the present results 
(Fig. 5) it is possible that the EL for the electropolished samples is determined 
chiefly by the first mechanism. But the difference between the electropolished 
and washed samples in the magnitude and field dependence of J,/J, can be 
most easily understood in terms of an increasing number of defects, and hence 
an increasing amount of breakdown. The erratic nature of the EL for washed 
samples, under conditions previously noted, also follows more naturally on 
the breakdown mechanism. On the collision mechanism, the decrease in 
brightness upon polishing can only be understood by making the more special 
assumption that the activating impurity is concentrated at the surface. 


5. CONCLUSIONS 

The chief conclusions may be summarized as follows: 

(1) Under the conditions of formation used here, the oxide behaves as a 
good dielectric. 

(2) The observed voltage corresponding to zero thickness of the film is in 
reasonable agreement with the theoretical value. 

(3) The films contain defects whose number and size increase with increasing 
surface roughness before anodization, and with increasing thickness. The 
presence of the defects probably leads to a randomly occurring breakdown 
of the oxide in the regions surrounding them. 

(4) The characteristics of the rectification and electroluminescence are 
probably determined by the defects. 
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PHOTONEUTRON CROSS SECTION OF Ni**! 


J. P. Roatsvic, R. N. H. HAsLtaAm, AND D. J. MCKENZIE 


ABSTRACT 


The reaction Ni®§(y,2)Ni®? has been investigated between 12- and 24-Mev 
maximum photon energy. The absolute yield value at 22 Mev was found to be 
much lower than reported previously. Also, a shorter half-life of Ni? was ob- 
tained. The cross section was determined using two different methods. 


1. RELATIVE ACTIVATION CURVE 

The activation curve and the cross-section curve for the reaction 
Ni®8(y,2) Ni5? have been determined from threshold (12.0 Mev) up to 24.0 
Mev. Nickel disks with diameter 2.5cm and thickness 0.15 mm were irradiated 
in the center of the X-ray beam from the University of Saskatchewan 24-Mev 
betatron at a fixed distance from the target. The disks were irradiated at 
different maximum betatron energies, the irradiation procedure being the same 
for all disks. 

In order to obtain a relative yield curve, no correction for the finite irradia- 
tion time was necessary. However, the relative dose given to each disk was 
measured by a monitor in connection with an R-C circuit whose time constant 
was equal to the mean life of Ni5’. Also the relative number of roentgens 
per dose was measured for different energies to find the relative number of 
roentgens given to each disk. This ratio changes slightly with the energy. 

The secondary activity from Ni? was followed for several days in a pro- 
portional flow counter. The disks were mounted on aluminum of sufficient 
thickness to obtain saturated backscattering. Since all disks were of the 
same thickness, no correction for self-absorption and backscattering was 
necessary in determining a relative activation curve. 

The activity in counts/minute, No, at the end of the irradiation was obtained 
from the decay curve of each sample. After the short-period activities (due 
largely to Ni®(y,p)Co") had disappeared, in about 24 hours, the least squares 
method was used to obtain the best-fitting straight line relating In N and t. 
The No was divided by the dose given to the sample and by the weight of the 
sample to give a point on the relative activation curve. 

The somewhat high yield value at 24-Mev maximum photon energy (see 
Fig. 2) could be explained by the onset of the reaction Ni®§(y,2”)Ni*® with 
threshold 23.2 Mev and a half-life of 6.4 days. However, the slopes of the 
decay curves of samples exposed at this energy did not yield any higher 
half-life than the others; nor did the points show any departure from a straight 
line even after 120 hours. 

Also the half-life of Ni®? was obtained by this method. It was found to be 
34.141.4 hours, which is slightly lower than the values reported previously 

1Manuscript received February 16, 1959. 

Contribution from the Physics Department, University of Saskatchewan, Saskatoon, 
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by Maienschein and Meem (1949), Friedlander et a/. (1950), and Rudstam 
(1956). Their values were 35.7+1 hours, 36.4+1 hours, and 37.6 hours 
respectively. 


2. ABSOLUTE YIELD VALUE AT 22 MEV 


The absolute value of the activation curve was determined at 22-Mev 
maximum bremsstrahlung. Nickel disks, with the same diameter as before 
but with thicknesses ranging from 10 to several hundred milligrams per 
square centimeter were irradiated 54.7 cm from the target, the betatron 
being operated at a maximum energy of 22 Mev. Recent investigations have 
shown that the energy stability is very good, so that deviations from 22 Mev 
were assumed to be less than 20 kev. 

A slightly different irradiation procedure was followed: The disks were 
irradiated for a time, tg, and the counting rates were corrected for finite 
irradiation time by multiplying by the factor tg(1 —e~*‘x)~!. The relative dose 
given to the sample was recorded by a monitor connected to a condenser and 
was given in arbitrary units (see Johns et a/. (1950)). The disks were counted 
in a flow counter as before and the counting rates at zero time obtained. 
These were corrected for self-absorption in the samples, backscattering, ab- 
sorption of the backscattered 8-particles, and geometry of the counter by 
using a method outlined in another paper (Roalsvig and Haslam 1959). 

The corrected number of counts were divided by the dose and by the number 
of moles* of Ni** in the disks to give saturated activity per mole per unit dose. 

The peaking effect of the X-ray beam, that is the distribution of intensity 
as a function of angular distance from the beam axis at the position of the 
disks, was measured by two different methods: (a) Small copper squares, 
lcmX1 cm X0.5 mm, were irradiated at the same time at different distances 
from the center of the beam. They were counted in the flow counter, and the 
relative counting rates at a fixed time were obtained. (b) A copper disk, 5 cm in 
diameter and 1 mm thick, was irradiated with its center at the center of the 
beam. It was then placed on a X-ray film for 1 hour. The film was developed 
and measured by a microphotometer. The relative intensity distributions 
F(r) found by these two methods were in agreement. F(r) is shown in Fig. 1, 
by the curve marked ‘position I’’. The ‘effective area”’ of the disks, defined as 
on F(r)2mr dr, was calculated, and the ratio of the true area to the 
“effective area’’ was found to be 1.104. 

The number of roentgens per arbitrary dose unit was obtained by using a 
“Farmer Sub-Standard Dosemeter’’, manufactured by the Baldwin Instru- 
ment Company, referred to as a Baldwin—Farmer r-meter. The instrument 
has been calibrated for Co™ radiation, the calibration factor being 1.08. This 
same factor was applied to obtain the dose in roentgens of the high-energy 
bremsstrahlung used in the present experiment. The thimble chamber of the 
Baldwin—Farmer r-meter was placed in the center of a lucite block, 8X88 
(cm), and given a dose corresponding to a certain number of arbitrary units. 


*By definition: | mole = 1 gram atom. 
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This was done several times at different distances from the betatron target. 
A 1/D* dependence was obtained for the roentgens per unit dose, where D 
is the distance from the target to the center of the thimble chamber. The 
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_Fic. 1. Spread of the relative intensity of the X-ray beam from its center at two different 
distances from the target. Position I: 54.7 cm; position II: 182.0 cm. 


number of roentgens per unit dose at the position where the nickel disks were 
irradiated was obtained from this relationship. These values were corrected 
for differences in pressure and temperature from the values (76 cm Hg and 
22° C) at which the instrument was calibrated. 

Since the Baldwin—Farmer r-meter, however, is placed inside a lucite 
block, the nickel disks should theoretically also be placed in the center of the 
lucite block (see Johns et a/. 1950; Katz and Cameron 1951). This introduces a 


correction factor, 


222 Mev 
{ og(y,n) P(E, 22 Mev) "dE 


threshold jccatinipetenansnenbaaietintiis 


222 Mev 
| on(y.n) P(E, 22 Mev) dE 
© threshold 
where og(y,) is the cross section at the energy E for the reaction Ni5*(y,m) Ni®’. 
Since only ratios are involved, relative values can be used. P(E, 22 Mev) is 
the number of photons per cm? per 100 r with energy between E+1/2dE and 
E—1/2dE, when the maximum photon energy is 22 Mev. These values were 
taken from Katz and Cameron (1951). 4(£) is the X-ray attenuation coefficient 
in the lucite of X rays with energy E. Lucite was taken as C;02Hs (i.e. 8% H, 
60% C, and 32% O). The values of u(£) were taken from Grodstein (1957); 
x is the thickness of 4 cm of lucite (4.13 g/cm*). The integrations have to be 
carried through from threshold to 22 Mev. The correction factor was calcu- 
lated to be 1.074. 

The absolute value obtained in this way is 3.17 X 107 n/100 r mole. However, 
this value has to be corrected due to electron capture. It was found by Fried- 
lander et al. (1950) and recently confirmed by Konijn et a/. (1956, 1958) that 
the ratio of the probability of electron capture to positron emission in Ni*? is 
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very nearly unity. The absolute activation value at 22 Mev is therefore 
0.63108 n/100 r mole. The absolute activation curve normalized to this 
value is given in Fig. 2. 


neutrons /iOOr/ mole nieve 
Oo O 
ND w Q 8 9 


Activations, 


oO 


13 15 17 19 2! 23 25 
Maximum photon energy (MeV) 


Fic. 2. Activation values for the reaction Ni55(7,2)Ni®? as a function of maximum photon 
energy. 


3. COMPARISON OF ABSOLUTE YIELD VALUES 

The absolute yield value was first measured by Price and Kerst (1950) for 
the natural element of nickel to be 1.04108 n/100 r mole at 22 Mev. Katz 
et al. (1951) obtained the value 1.12X10° n/100 r mole for the reaction 
Ni®§(y,2) Ni? alone, also at 22 Mev. This value is much higher than the 
present one. Skarsgard (1955) measured the ratio of the yield of Cu®(y,2)Cu® 
to the yield of Ni®8(y,7) Ni”? at 23 Mev to be 3.4. 

In order to compare this ratio with the present work, the absolute yield of 
the reaction Cu™(y,7z)Cu® at 22 Mev was redetermined, using the same 
method as for nickel. The correction factor for the absorption in the lucite 


block is essentially the same, being 1.075 for this reaction. Also, the copper 


disks were irradiated at two distances from the target: (1) 54.7 cm as for 
nickel, and (11) 182.0 cm. The peaking effect was also measured in the second 
position and is given in Fig. 1 (curve marked ‘“‘position II’’). The correction 
factor was 1.007 in this position. After these corrections were applied, the 
two saturated activities showed a 1/D? attenuation to within 0.49%, where D 
is the distance from the target to the disk. The absolute yield value obtained 
was 2.25 10° n/mole 100 r, which is somewhat lower than the value 2.51 10° 
n/mole 100 r used previously in this laboratory (Katz and Cameron 1951). 
The activation curve of the reaction Cu®(y,2)Cu® was also remeasured from 
20 to 23 Mev. Excellent agreement was found with the older curve of Johns 
et al. (1950) when the latter was normalized to the present value at 22 Mev. 
The half-life of Cu® which was obtained from our experiment was (9.72 
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+0.03) minutes. This value agrees well with the value of Berman and Brown 
(1954), (9.73+0.02) minutes, but is lower than the value 9.9 minutes as 
given by the General Electric Chart of the Nuclides 1956. 

The ratio of the yield of Cu®(y,u)Cu® to the yield of Ni®§(y,2)Ni®? was 
calculated to be 3.5 at 22 Mev, and 3.4 at 23 Mev, in excellent agreement with 
Skarsgard. Titterton (1958a) has also measured this ratio at 30 Mev to be 3.3. 
Since the cross-section curve of Cu®(y,2)Cu® has its maximum at a lower 
energy than that of Ni5*(y,7)Ni*’, a lower value of the yield ratio is expected 
at higher energy. Titterton’s result is thus in better agreement with the 
present value of the ratio than with the old one of Katz et a/. (2.2 at 22 Mev). 

The difference between the present activation value and that previously 
reported from this laboratory may be due to several factors. The calibration 
of the betatron energy scale is more reliable than it was at the time of the 
earlier investigation (1951). Longer irradiations and the use of improved 
counting methods have made better counting statistics possible. The decay 
curves were also followed for longer times to make sure that short-period 
activities had died out. It seems probable that some contribution from the 
(y,P) reaction in Ni® might have been included in the earlier work with a 
consequent increase in the observed yield. Perhaps the most important factor 
is an improved knowledge of the corrections to be applied for self-absorption 























and backscattering. 







t. CROSS SECTION 









Two different methods were used in calculating the cross-section curve. 






These were: 
(a) the ‘‘photon difference method” of Katz and Cameron (1951), and 
(6b) the ‘modified spectrum method” of Penfold and Leiss (1958). 

The two curves are given on Fig. 3 and the characteristic values in Table I. 
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Fic. 3. Cross-section curves for the reaction Ni®8(y,n)Ni*? as a function of photon energy 
Curve marked K and C: photon difference method; curve marked P and L: modified spectrum 
method. 
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The general shapes of the two cross-section curves are the same. They both 
have their maxima at 20.5 Mev, and their half widths are nearly the same. 
However, the absolute value of om,x by using the ‘‘photon difference method”’ 
is about 40% higher than the value obtained by the ‘modified spectrum 
method”’. 


TABLE I 


Yield and cross-section data for the reaction Ni(y,m) Ni®? 








Yield at Ratio of yields, 





24 Mev 
22Mev, {Cu®(yn)Cu"l/ Em, om Tym Jo %nmdb, 
References n/100rmole = { Ni58(y,2) Ni57} Mev mb Mev Mev-barn 
Price and Kerst 
(1950) 1.04 1054 
Katz et al. (1951) 1.12108 2.2at 22 Mev 18.5 60 4.6 0.33 
Katz and Cameron 
(1951) 18.5 54 5.6 0.34 
Katz and Cameron 
renormalized (0.63 X 108) 18.5 30 5.6 0.18 
Skarsgard (1955) 3.4at 23 Mev 
Titterton (19582) 3.3 at 30 Mev 
Titterton (1958)) 19.0 32 4.4 0.22 (to 32 
Mev) 
Present work 0.63 X 108 3.5at 22 Mev 20.5° 29 Tie 0.20 


20.5° 21 8.0 0.16 


*For natural element. 
’Katz and Cameron's method. 
©Penfold and Leiss’ method. 


A comparison between the present work and that of Katz and Cameron 
(1951), and of Titterton (19580), is presented in Table I. This indicates that 
the cross-section curve in this work is much flatter and has a lower maximum 
value. Also, this maximum value is found at a higher energy than reported 
previously. However, good agreement is found for the integrated cross sections. 
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STUDY OF PHASE PROPAGATION IN SUPERCONDUCTING 
ALUMINUM BY ULTRASONIC ATTENUATION MEASUREMENTS! 


K. L. CHopra AND T. S. HuTCHISON 


ABSTRACT 


The phase propagation in superconducting aluminum has_ been studied 
by measuring the time rate of change of ultrasonic attenuation. The time taken 
for the destruction of the superconducting phase in a cylindrical specimen, by 
means of a magnetic field, H, greater than the critical field, H., is approximately 
proportional to {H/(H—H,)} in agreement with eddy-current theory. In the 
converse case, where the superconducting phase is restored by switching off 
the magnetic field H (>H,), the total time taken is nearly independent of the 
temperature (or H,) as well as H. The superconducting phase grows at a non- 
uniform volume rate which is considerably less than the uniform rate of collapse. 


INTRODUCTION 


It is now well known (B6Ommel 1954; MacKinnon 1955, 1957; Mason and 
Bémmel 1956; Morse, Tamarkin, and Bohm 1956; Chopra and Hutchison 
1958) that when a metal specimen undergoes a transition from normal to 
superconducting phase the high-frequency ultrasonic attenuation decreases 
rapidly. This is shown for polycrystalline aluminum in Fig. 1, which is repro- 
duced from our previous paper (Chopra and Hutchison 1958) where the 
temperature measurements below 1° K were in error and have been corrected. 
Thus the change in ultrasonic attenuation provides another tool to detect 
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Fic. 1. Ultrasonic attenuation vs. temperature in polycrystalline aluminum. 
‘Manuscript received February 4, 1959. 
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the superconducting transition. Rates of growth of a phase following a mag- 
netically induced transition are readily followed. It is clear that this method 
gives a measure of the volume rate of growth of the new phase in a localized 
region determined by the path of the sound beam in the specimen. 

Faber (1952, 1953, 1954) has studied thoroughly the nucleation process 
which initiates the phase transition and the subsequent growth of the phase 
in the case of tin. We shall be concerned here with the latter aspect of the 
phase transition, namely the propagation of the phase. According to the 
theory of Pippard (1950) and Lifshitz (1950), the rate at which the super- 
conducting phase collapses radially in cylindrical specimens is controlled by 
an electromagnetic damping associated with the eddy currents. The latter 
are generated due to the displacement of magnetic lines of force through the 
normal phase during the movement of the interphase boundary. If the 
cylindrical core of the superconducting phase is contracting radially and 
symmetrically in a long rod of radius, 7, under the influence of an external 
longitudinal field H (>H,), the time it takes to collapse completely is pro- 
portional to or?}H/(H—H,)} where o is the effective electrical conductivity 
of the normal phase. Faber (1953) has checked this result for various tin 
specimens. 

In the converse case, when the superconducting phase is allowed to expand 
radially, the process is expected to be asymmetrical and the above-mentioned 
eddy-current theory not valid. Faber has explained this on the basis of the 
general principle that the closer the boundary is to the outside surface of a 
specimen, the less difficulty it experiences in displacing flux through the 
intervening normal phase. Thus the boundary moving inwards will be quite 
stable while the outward-moving boundary will be composed of a tortuous 
arrangement of radially expanding superconducting filaments. 

This paper presents results on the average volume rate of growth of the 
superconducting or normal phase in cylindrical specimens of polycrystalline 
aluminum. Results on a single crystal of different geometry are also men- 


tioned. 


EXPERIMENTAL DETAILS 

Measurements of ultrasonic attenuation are made by the conventional pulse 
technique. A pulse of 20 Mc/sec longitudinal waves is axially applied to one 
end of the cylindrical specimen by a }-in. quartz crystal which is also used 
to receive the reflected pulses. After amplification, the pulse and its echoes 
are displayed on an oscilloscope and photographed with a 16-mm Cine Kodak 
Special Movie Camera. The rate of change of the attenuation is measured 
by the rate of change of the logarithm of the amplitude of one of the decaying 
echo pulses. A pair of Helmholtz coils (diameter, 12 inches) around the liquid 
air Dewar provides a variable magnetic field. Temperatures down to 0.88° K 
are obtained in a booster cryostat (Atkins et a/. 1955). 

The polycrystalline cylindrical specimens are prepared from aluminum of 
99.99,% purity. They are recrystallized from the cold-rolled condition by 
annealing at 500° C for 5 hours. This process gives a grain size of ~3 mm. 
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The single crystal specimen is obtained from aluminum of similar purity and 
is of cubical shape with side 20 mm. 


RESULTS 

Two simple cases have been studied. In the first case, the attenuation is 
measured as a function of time when the superconducting phase of the cylin- 
drical specimen collapses radially under the influence of an applied magnetic 
field H (>H,). In the second case the time rate of change of attenuation is 
measured when the superconductivity of the cylindrical specimen is restored 
by removing the applied magnetic field H (>4H,). 

Figure 2 shows a plot of attenuation coefficient against time for a poly- 
crystalline cylindrical specimen of diameter 20 mm and length 40 mm. The 
initial time taken for the attenuation to show any change is the time taken 
for the collapsing process to reach the region traversed by the 3-in. cylindrical 
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Fic. 2. Ultrasonic attenuation coefficient as a function of time for radially collapsing super- 
conducting phase. 


Fic. 3. Total time taken for collapsing of superconducting phase as a function of 
{H/(H—H,)}. 
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sound beam. The linearity of the graph indicates that the gross rate of change 
of superconducting volume in the }3-in. cylinder is constant and that the 
process is symmetrical. Similar volume rates are obtained with the single 
crystal which is probably due to the poorer geometry. 

The total time taken for the collapse of the superconducting phase is 
measured as a function of the temperature and applied magnetic field. Figure 3 
shows a plot of time against {H/(H—H,)}. The straight-line relations verify 
the validity of the eddy-current theory. Marked departures from this relation 
occur for H close to H,. The dependence of the slopes of the straight lines 
on temperature indicates that the electrical conductivity o is a function of 
temperature. 

It is clear that the average volume rate of propagation of the supercon- 
ducting phase in the above case is a function of both temperature (or H,) 
and the magnetic field H. At 7 = 0.88° K and H = 97 oersted, the rate is 
5 cm?/sec. 

In the converse case, a plot of attenuation coefficient against time for the 
same specimen is shown in Fig. 4. The non-linear curve indicates that, as 
expected, the growth of the superconducting phase is not symmetrical. More- 
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Fic. 4. Ultrasonic attenuation coefficient as a function of time for the growth of super- 
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over, the total time taken for the restoration of the superconducting phase is 
nearly independent of temperature (or H,) as well as applied magnetic field 
H. The average rate of propagation is about 0.5 cm*/sec. The single crystal 
gives similar rates—this probably being due to the poorer geometry in the 
single crystal than in the polycrystalline cylinder. 

It should be mentioned that in thicker specimens (for example, one with 
diameter 30 mm and length 40 mm) the plot of attenuation coefficient against 
time shows steps as shown in Fig. 5. This seems to indicate the movement 
of domains at different rates. 


CONCLUSIONS 


Measurements of ultrasonic attenuation have been successfully utilized to 
study the gross rate of propagation of the superconducting or normal phase. 
The radial collapse of the superconducting phase has been found to be sym- 
metrical and to obey eddy-current theory. 

In the converse case, when superconductivity is allowed to expand radially, 
the process is found to be asymmetrical. In fact, the rate being independent 
of temperature as well as magnetic field, the transition process appears to 
be quite devious. 
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ON TOROIDAL FUNCTIONS! 
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ABSTRACT 


The general solutions of Laplace’s equation in toroidal co-ordinates have been 
derived, both for the three-dimensional and the axial symmetrical cases. Also, 
the general theory and properties of toroidal functions have been briefly outlined. 

The numerical tables for the zonal and the tesseral toroidal functions were 
calculated by the I.B.M. 650 Computer at the University of Toronto, and are 
included in the present paper. 









INTRODUCTION 


Toroidal functions, which arise when Laplace’s equation is separated in 
toroidal co-ordinates and are suitable for problems in which the boundary 
conditions are given over the surfaces of anchor rings, were first investigated 
by Neumann (1864), Hicks (1881, 1884), Basset (1893), Niven (1892), and 
others. Although Relton (1931, 1933) and Fock (1932) emphasized the desira- 
bility of tabulating the toroidal functions which are essentially the Legendre 
functions of half odd integral order, and also suggested that the tabulation of 
the functions should be through the hypergeometric functions, in terms of 
which they may be expressed, the author has only been able to trace one 
numerical table of the zonal toroidal functions, namely that published by 
Fouquet (1937), apart from the numerical table calculated by Carter and the 
author (1958). 

Although toroidal functions were studied many years ago, the lack of 
numerical tables greatly limited their application to physical problems. There 
are, however, a number of physical problems where toroidal functions may 
be found useful, such as the potential of a ring, the electric potential of a tore 
and its capacity, the velocity potential for a tore moving parallel to its axis, 
as well as the energy of the motion, and the self-inductance of a tore. 

In the present paper, a brief discussion of the general theory of toroidal 
functions is presented and a complete set of numerical tables for zonal toroidal 
functions and tesseral toroidal functions has been computed and included. 
In addition, expressions for the approximate values of p,(u) and q,(u) are also 


























derived. 








TOROIDAL CO-ORDINATES 

Toroidal co-ordinates are a natural co-ordinate system to describe a ring. 
The relationship between the toroidal co-ordinates (u, v, w) and the Cartesian 
co-ordinates (x, y, 2) is given by 
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cosh u—cos v 


| see a sinh u cos w 


a sinh uw sin w 


™ » = cosh u—cos 0” 


z= ae 


cosh u—cosv’ 





where 0 <u < ~,0 <u < 2m, and 0 < w < 2n. The surface u = up is an 
anchor ring with axial circle in the x—y plane, centered at the origin, of radius 
acoth wu, having a circular cross section of radius acsch up. The surface 
v = vo (for vo < m) is that part of the sphere of radius @ csc vo, with center at 
x = y = 0, 2 = a cot %, which is above the x—-y plane, and the rest of the 
same sphere below the x—y plane, is the surface v = 247—vp (see Fig. 1). The 
limiting tore u = © is an infinitely thin ring of radius a. It is evident that 
the value of a can be so chosen that any given anchor ring will be included in 
the system. 





Fic. 1. Toroidal co-ordinates. 


For the element of length ds we have in the new co-ordinates u, v, w, 
(2) (ds)? = (h, du)?+(h, dv)?+(h, dw)? 


which is easily found to reduce to 
a” ‘ $4 2 se 2) 
aaa | (du)°+(dv)°+sinh” u(dw)*}. 


(3) (ds)? = 


Expressed in terms of the toroidal co-ordinates, Laplace’s equation becomes 


(4) 0 ( sinh u _ ie) 42 ( sinhu | $8) 
“ @u \cosh u—cosv Ou dv \cosh u—cosv dv 


a Se, 
(cosh u—cos v) sinh u dw 


= 0. 
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TESSERAL TOROIDAL FUNCTIONS 


In order to reduce equation (4) to a form for which integrals can be ob- 
tained, let ¢ = ¥/(cosh u—cos v); Laplace’s equation then becomes 
aya is a 
sat! a ~——% aie u 42 y= 


dv" sinh? u dw 









(5) 













Let us endeavor to satisfy the equation by means of a function y = U-V-W, 
where U is a function of u only, V of v only, and W of w only; substituting in 
the differential equation and dividing throughout by U- V-W, we have 


6) fee sae a Lalas i? i 2 
U \ du V dv" ! sinh? u dw ” 

In this equation the first term involves u only, the second v only, and the 
third w only; it is consequently obvious that the only case in which it can be 
satisfied is when each term is equal to a constant, the sum of the three con- 
stants being equal to zero. We thus obtain three ordinary differential equations 


(7) oO a! ‘Ww =0, 














(8) as s+n°V = 0, 










and 





fu. ( 1 <) = 
ero ’ du P 47 sinh? w _ 


The independent solution of the equation (7) is 


(9) 










(10) W = a, cos mwt+b,, sin mw 





where a,, and },, are arbitrary constants. 
Similarly, the solution of the equation (8) is 






(11) V = c, cos nv+d, sin nv 





where c, and d, are arbitrary constants, 
In the theory of toroidal functions, m is any positive integer, and m is any 
positive integer which is not greater than , the value zero being included. 
The differential equation (9) is satisfied by 







"(4) = 1 (—1)"C(n+}3) J cos m6 dé 
Pn ~ e T(n—m+3) Jo (cosh u+sinh u cos gy 
) 
da ‘Sta (—1)"P(n+}3) cosh mé dé 
dn ~  T(n—m+}3) 0 (c osh u+sinh u cosh 6)"* 





where p,; (uw) and g;'(u) are known as ‘‘tesseral or associated toroidal functions”’ 
of first and second kinds of degree m and order m. 

Using equation (12) the following tables (I and II) of values for both kinds 
up to the fourth order have been worked out by the I.B.M. 650 Computer at 
the University of Toronto. The numbers in parentheses in the tables indicate 
the power of 10 by which tabulated values are to be multiplied, i.e., p2(1.0) = 
0.95216. 
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The tesseral toroidal functions py'(w) and gf'(u) are essentially the associated 
Legendre’s functions of first and second kinds of degree (n—}) and order m 
respectively; the relations between p;'(u) and g;'(u), and the corresponding 
associated Legendre’s functions are 

bu (u) = Pry (cosh x) 


(13) and 
Gn (u) = Qn-y (cosh x). 


The new notations are suggested by Carter and Loh (1958). 
The required solutions of Laplace’s equation are 


c 


/ (cosh u—cos v) pr’ (u) Sin 2V Sin Mw 


(14) and 
/ (cosh u—cos v) gr (u) Sin 20 Sin MW. 


The spaces interior to and exterior to a tore are not simply connected; and 
thus the toroidal functions are not always directly applicable to the solution 
of potential problems for those spaces in which circulation is involved, as the 
potentials in such cases are not always uniquely determined by their values 
on the surface of the tore. It may be shown easily that, for the space exterior 
to the given tore (uw = m,)—that is, 0 < u < u,;—the p-functions must be used 
as potential functions. Similarly, for the space interior to the given tore—that 
is, “4, <u < o—the g-functions are chosen as potential functions. Heine 
(1881) has mistakenly stated that the p-functions are to be adapted as poten- 
tial functions to the interior space and the q-functions to the exterior space. 

However, the general solution of Laplace’s equation in toroidal co-ordinates 
(uw, v, w) may be given by the series, 

(15) @ = V(cosh u—cos v) . z (dm cos mw+b,, sin mw) 
r 


m 
° m m 
X (c, cos nv+d, sin nv) X {Amn Pn (U)+Bmn Gn (u)}, 
where Gm, Om, Cn» Any Amn, Bmn are arbitrary constants to be determined by 


the boundary conditions. 


ZONAL TOROIDAL FUNCTIONS 
For the particular problems having cylindrical symmetry, the variable w 
becomes irrelevant. Expressed in terms of co-ordinates “ and v, Laplace’s 


equation becomes 


0 sinh u at) 0 ( sinh u 4) 

6 -{-- 50 <? +5 ae Py) = 0. 
(16) Ou 3 u—Cosv Ou T5p cosh u—cos v dv 
Again we substitute, 
(17) @ = ¥ V (cosh u—cos 2). 
Equation (16) reduces to 

ay, ay dy 1 

3 ) oe > 5 h = = =o) 

(18) we ay? t cot u ta y= 0 
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This is satisfied by a function of the form 


(19) y = F(u)-G(v) 
if 






€6. 36 =0 
dw 






(20) 






and 





TE coh u E+ (tn) F = 0, 

The first equation is satisfied by G = cos(mv+a,) and the second by p,(x) 
and g,(u); pa(u) and q,(#) are known as the zonal toroidal functions of first 
and second kinds respectively. 

The general solutions of Laplace’s equation in toroidal co-ordinates having 
axial symmetry are of the form 






(21) ¢ = V (cosh u—cos v) . {An Pn(u)+By dn(u)} cos (nv+an) 





where Ay, By, a, are arbitrary constants. The constants are to be chosen so as 


to fit the given boundary conditions. 
The zonal toroidal functions of the two kinds are here denoted by 


oe 10 
Pr(u) = | ‘ ati ’ 


mw Jo (cosh u—sinh u cos 0 








qn(u) = | OF ay 
J» (cosh u+sinh u cosh 6) 





For the particular cases, where » = 0 and n = 1, it has been shown by 
Hicks (1881) that 


po(u) = 2\/k' F(k) = 26°” F(R), 
92 










(23 & . u/2 yp 
ss pilu) Ve E(k) = 2 e“”” E(k), 
and 
(24) go(u) “a 2\/k’ F(k’) =e (u/2) F(k’), 
2 2 “(Be >/L? 
(25) gi(u) = a [F(k) —E(k)] = 2c” [F(k’) —E(k’)] 







where 





F and E are the complete elliptic integrals of the first and second kinds, and 
k?, the square of the modulus. 

Starting with the values of po(u), pi(u), go(u), qi(u), the zonal toroidal 
functions of higher order can be derived from the recurrence formulae (Hicks 


é¢,.? 


1881), argument ‘‘z’’ understood: 












(26) (2n+1) Paii—4n cosh u pyp+(2n—1) par = 0 








626 CANADIAN JOURNAL OF PHYSICS. VOL, 37, 1959 


and 
(27) (2n+1) dn41—4m cosh u ga +(2n—1) qr = 0. 


Using equations (23), (24), (25), (26), and (27), the numerical tables (III 
and IV) for the first four functions of each kind have been worked out by 
Carter and Loh (1958) in connection with their investigation of the electric 
field between a rod and a concentric ring. These numerical tables have been 
checked by the I.B.M. 650 Computer at the University of Toronto. There 
are differences between these figures and those contained in Fouquet’s table 


TABLE III 
Zonal toroidal functions of the first kind 


u ph(u) PY(u) b3(u) p3(u) p9(u) 

0 1.0000 1.0000 1.0000 1.0000 1.0000 
0.1 (—1) 9.9938 1.0019 1.0094 1.0220 1.0397 
0.2 9.9751 1.0075 1.0378 1.0893 1.1635 
0.3 9.9440 1.0169 1.0858 1.2060 1.3856 
0.4 9.9009 1.0300 1.1547 1.3793 1.7314 
0.5 9.8460 1.0469 1.2459 1.6196 2.2410 
0.6 9.7795 1.0676 1.3615 1.9419 2.9739 
0.7 9.7021 1.0921 1.5043 2.3659 4.0161 
0.8 9.6141 1.1205 1.6776 2.9176 5.4910 
0.9 9.5161 1.1526 1.8852 3.6311 7.5740 
1.0 9.4086 1.1887 2.1320 4.5506 (1) 1.0515 
| 9.2924 1.2287 2.4236 5.7330 1.4667 
1.2 9.1679 1.2726 2.7667 7.2518 2.0533 
is 9.0359 1.3206 3.1692 9.2015 2.8826 
1.4 8.8971 1.3727 3.6401 (1) 1.1704 4.0554 
1.5 8.7521 1.4289 4.1902 1.4914 5.7151 
1.6 8.6015 1.4895 4.8321 1.9033 8.0648 
1.7 8.4460 1.5544 5.5802 2.4319 (2) 1.1393 
1.8 8.2863 1.6238 6.4515 3.1103 1.6108 
1.9 8.1229 1.6978 7.4659 3.9808 2.2790 
2.0 7.9565 1.7765 8.6462 5.0980 3.2262 
2 7.7876 1.8602 (1) 1.0019 6.5320 4.5691 
£2 7.6167 1.9489 1.1616 8.3727 6.4734 
2:3 7.4444 2.0429 1.3472 (2) 1.0736 9.1741 
2.4 791 2.1423 1.5631 1.3769 (3) 1.3005 
2.5 7.0973 2.2474 1.8139 1.7663 1.8438 
2.6 6.9234 2.3584 2.1055 2.2662 2.6147 
27 6.7497 2.4756 2.4444 2.9080 3.7082 
2.8 6.5766 2.5992 2.8382 3.7320 5.2597 
2.9 6.4045 2.7296 3.2958 4.7900 7.4609 
3.0 6.2337 2.8669 3.8276 6.1484 (4) 1.0584 
Sut 6.0643 3.0116 4.4455 7.8925 1.5016 
32 5.8967 3.1639 5.1635 (3) 1.0132 2.1304 
3:3 5.7311 3.3243 5.9978 1.3007 3.0226 
3.4 5.5677 3.4931 6.9671 1.6699 4.2886 
3.5 5.4066 3.6708 8.0934 2.1439 6.0851 
3.6 5.2480 3.8578 9.4021 2.7525 8.6344 
3.7 5.0920 4.0545 (2) 1.0923 3.5340 (5) 1.2252 
3.8 4.9388 4.2614 1.2689 4.5375 1.7385 
3.9 4.7884 4.4790 1.4742 5.8259 2.4669 
4.0 4 4 1.7126 7.4802 3.5006 


.6408 .7079 
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TABLE IV 


Zonal toroidal functions of the second kind 


a9(u) g3(u) f a9(u) 


a8 (u) 

@ @ co @o 
.3799 2.3896 1.7421 1.3675 117 
6805 1.7124 1.1022 (—1)7.7142 6171 
. 2669 1.3295 (—1)7.6406 4.8023 3.1482 
9693 1.0691 5.5122 3.1202 8452 
7350 (—1) 8.7684 4.0668 2.0763 . 1088 


NNwOw 


5404 7.2828 3.0433 1.4027 7680 
3733 5.1010 2.2995 (—2)9.5739 1755 
2262 5.1425 7495 6.5825 2.5956 
0947 4.3543 .3378 4.5501 .6225 
9754 3.6995 .0270 3.1581 -0185 


=m tw bw 


. 8661 3.1515 -9066 2.1987 4135 
-7654 2.6902 6.1013 1.5345 .0489 
.6719 2.3001 .7166 1.0729 2.5611 
5846 .9691 6515 (—3)7.5135 -6225 
5030 .6875 2.8300 5.2679 -0292 


ee a 


4263 4473 1954 3.6970 ( 6.5344 
3542 . 2422 .7044 2.5965 .1522 
. 2861 .0667 3239 1.8247 2.6402 
. 2219 1635 .0289 . 2830 6796 
-1611 .8749 .9993 (—4) 9.0250 0690 


OD mt eet es 


3.7694 5.2211 6.3504 5. 8058 
8205 4.8394 4.4697 .3342 
.0055 3.7654 3.1466 .7608 
3053 2.9303 2.2156 -7589 
7035 2.2807 1.5603 . 1208 


. 1035 
0490 
9726 
9.4819 
9.0160 


com 


2. 
2 
2 
2. 
2. 


1862 .7753 1.0989 (—6)7.1424 
7413 3820. = (—5) 7.7406 4.5521 
3588 0759 5.4528 2.9015 
0297 3771 3.8414 1.8495 
.7466 ).5220 2.7063 1.1790 


5737 
. 1535 
7542 
3748 
0142 


tN tw 


—wnrwhds w 


won do 


5.6714 5031 5.0789 1.9067 5162 
5. 3454 . 2935 3.9548 1.3434 7917 
5.0355 .1132 3.0796 (—6) 9.4658 3.0549 
7408 5808 2.3982 6.6697 9477 
4605 2456 8675 4.6997 .2418 


wwe 


ww 


5.1940 .0966 4543 3.3116 (—8)7.9174 
.9405 5.1078 . 1326 2.3335 5.0481 
4.6994 5.2568 (—5) 8.8201 1.6443 3.2186 
.4701 4.5244 6.8689 1.1587 2.0522 
. 2520 3.8941 5.3493 (—7) 8.16417 1.3085 


in which there are undoubtedly a number of errors. In addition, the numerical 
table for the first derivative of the zonal toroidal functions p,(u), gi(u) has 
been computed and included in the paper. The relations between p),(u) and 
q,(u) and the zonal toroidal functions p,(u) and q,(u) are 

dp,(u) _ 2n+1 


du Zeinh g MOO), 


dq, (u 2n+1 
tank) a ee (Qn+1—COsh U gn). 


(28) 
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Using the above expressions derived, the numerical tables for the first deriva- 
tive of toroidal functions have been calculated and are shown in Tables V and 
VI. 


TABLE V 


The first derivative of zonal toroidal functions of the first kind 





u bo (u) P,(u) po(u) p,(u) p(u) 

0 0 0 0 0 0 

0.1 (—2) —1.2484 (—2)3.7512 — (1) 1.8824 (—1) 4.4194 (—1) 8.0245 
0.2 — 2.4886 7.5024 3.8081 9.1081 (0)* 1.6967 
0.3 —3.7117 (—1) 1.1261 5.8220 (0)* 1.4353 2.7857 
0.4 — 4.9097 1.5027 7.9708 2.0474 4.1963 
0.5 — 6.0753 1.8803 (0)* 1.0304 2.7849 6.0943 
0.6 —7.2012 2.2594 1.2874 3.6933 (1) 1.1509 
0.7 —8.2813 2.6405 1.5741 4.8287 1.2342 
0.8 —9.3101 3.0240 1.8971 6.2619 1.7444 
0.9 (—1) —1.0283 3.4106 2.2636 8.0819 2.4628 
1.0 —1.1196 3.8009 2.6819 (1) 1.0402 3.4769 
1.1 — 1.2046 4.1957 3.1616 1.3368 4.9105 
1.2 — 1.2832 4.5959 3.7136 1.7165 6.9392 
1.3 — 1.3551 5.0024 4.3502 2.2030 9.8117 
1.4 —1.4205 5.4161 5.0860 2.8269 (2) 1.3881 
1.5 —1.4791 5.8383 5.9373 3.6272 1.9647 
1.6 —1.5313 6.2699 6.9236 +. 6543 2.7821 
1.7 1.5770 6.7123 8.0669 5.9724 3.9411 
1.8 ~1.6165 7.1668 9.3931 7.6643 5.5845 
1.9 — 1.6499 7.6346 (1) 1.0933 9.8362 7.9152 
2.0 - 1.6776 8.1171 1.2719 (2) 1.2625 (3) 1.1221 
2 — 1.6998 8.6158 1.4793 1.6204 1.5910 
2.2 — 1.7168 9.1322 1.7201 2.0800 2.2563 
2.3 1.7289 9.6678 1.9998 2.6702 3.2001 
2.4 Bir Elian = 2.3247 3.4278 4.5390 
2.5 — 1.7395 1.0803 2.7021 4.4006 6.4388 
2.6 — 1.7387 1.1406 3.1404 5.6497 9. 1343 
2.7 — 1.7342 1.2035 3.6497 7.2535 (4) 1.2959 
2.8 1.7263 1.2691 4.2413 9.3128 1. 8386 
2.9 1.7153 1.3377 4.9286 (3) 1.1957 2.6087 
3.0 ~1.7014 1.4095 5.7270 1.5352 3.7014 
3.1 —1,6851 1.4845 6.6546 1.9711 5.2519 
3.2 1.6664 1.5631 7.7323 2.5308 7.4521 
3.3 1.6456 1.6455 8.9844 3.2495 (5) 1.0574 
3.4 1.6230 1.7319 (2) 1.0439 4.1724 1 1066 
$.5 1.5987 1.8224 1.2129 5.3573 2.1291 
3.6 1.5730 1.9174 1.4093 6. 8787 4.0213 
3.7 1.5461 2.0171 1.6374 8.8324 4.2872 
3.8 1.5181 2.1218 1.9024 (4) 1.1341 6.0838 
3.9 1.4891 2.2316 2.2103 1.4562 8.6331 
10 1.4594 2.3470 2.5681 1.8697 (6) 1.2251 


For the general properties and a discussion of toroidal functions the reader 
should consult the articles by Hicks and Basset. 
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TABLE VI 
The first derivative of zonal toroidal functions of the second kind 


u a (u) q, (u) ay (u) a, (u) a, (u) 


© © © co 
0044 —9 8756 —9.5660 9.1781 — 8.7448 
-O711 —4.8024 —4.3815 3.9147 —3.4496 
3.4243 —3.0821 —2.6146 .1514 — 1.7371 
2.6062 — 2.2076 —1.7279 .3020 (—1) —9.5895 
2.1178 — 1.6755 —1.2040 . 2784 —5.5481 


7932 —1.3171 -6586 — 5.4208 —3.3005 
-5613 — 1.0595 5.3568 —3.6179 — 1.9992 
.3868 -6615 -7336 — 2.4466 — 1.2261 
. 2501 —7.1633 -5611 —1.6701 (—2) —7.5863 
1395 —5.9756 2.6993 —1.1480 — 4.7249 


-0476 5.0174 2.0577 (—2) —7.9326 — 2.9574 
.6971 . 2342 .5755 — 5.5036 — 1.8581 
9.0235 3.5873 . 2104 —3.8301 —1.1709 
4321 3.0486 3243 — 2.6721 (—3) —7.3953 
9062 2.5972 . 1980 — 1.8678 —4.6795 


4331 2.2169 .5660 — 1.3076 — 2.9654 
0037 . 8953 .3098 (—3) —9.1659 —1.8814 
}.6107 .6224 3.3407 — 6.4313 —1.1947 
. 2487 . 3902 2.5918 —4.5162 (—4) —7.5927 
9134 . 1922 .0122 —3.1735 —4.8283 


5.6015 0231 5632 —2.2311 —3.0719 
.3100 2 7849 .2149 — 1.5693 — 1.9552 
.0369 5464 4454 —1,1041 —1.2449 
. 7802 4847 .3459 (—4) —7.7708 (—5) —7.9282 
5386 5741 5.7140 —5.4703 — 5.0503 


.3106 - 4.7926 4464 ~ 3.8516 
0952 -1214 3.4603 2.7123 
8915 3.5447 2.6932 —1.9103 
6986 3.0492 2.0964 — 1.3455 
5158 - 2.6232 6320 (—5) —9.4780 


3425 2. 2560 — 1.2706 — 6.6770 
1780 -1,.9419 (- —9 8926 ~4.7040 
0219 — 1.6710 —7.7026 —3.3142 
8736 4379 5.9977 — 2.3350 ( 
7327 2374 6703 - 6453 


5987 0649 3.6368 1593 
4718 (—3 9.1646 2.8320 —8. 1684 
3509 8872 2.2054 -5.7558 
2360 1.7880 7175 4.0558 (— 
1267 - 8422 3375 — 2.8579 


APPROXIMATE EXPRESSIONS FOR ?,(u) 
The zonal toroidal functions of the first kind may be written as 


ag de ; 
w J» (cosh u—sinh u cos 6)" a 


(29) pb, (u) = 


Two special cases will now be considered. 
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(1) For small values of u (u < 1.0) 
If the above integral is expanded as a power series of cos 6, term-by-term 


integration leads to 
1 rs : ‘ 
(30) Pr(u) = 2 (cosh u)~"~? yet (n+1)(n+3) tanh’ 5+ ne :) : 


The substitutions 


cosh u = (1+3 w+...) 


and 
(31) tanhuw =ut+... 


reduce equation (30) to 
4n°—1 4 
(32) Pru) = (4 = ee ’) 
16 


neglecting the terms of order higher than the second. 


; LG. 

po(u) = (:-4.4) 

pi(u) = (43 .’) 
16 


o 15 ‘) 
p2(u) = +8. 


The approximate values for p,(u) having small arguments given by equation 
(32) have been calculated and plotted in Fig. 2. This approximate expression 
derived for p,(u) for small values of u is surprisingly accurate, even beyond 
the expected range of validity. 

(it) For large values of u (u > 1.5) 
By replacing cosh uw and sinh u with their exponential forms in equation 


(29), it is seen that 
(33 fi I { d6 
aad Pru) = wJdo [i(e'+e “)—4(e'—e “cos 0}""* 


The substitution 


reduces the above integral to 


le 2 (n+4)u {~ do’ 
( 5¢ = ~ € 3 iu SCD 
34) pa(u) T Jo (cos'6’+e ™ sin’d’) 


nt f- 


By comparison with the approximate expression for the complete elliptic 
integral, 


e dx et 
om j : r= 2172 = log (4/k), if k? is small. 
Jy (cos x+k* sin’x)'”” g (4/k) 
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We then find 
po(u) = =e ™ (loge 4+u) = 


Pile) : eM, 
2-4-6 eure (2a —2) (n—})u 
‘3 


Pu(u) 3 ~ (2u—1) ° 


The approximate values have been computed and plotted in Fig. 2. 
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APPROXIMATE EXPRESSIONS FOR gqn(u) 
The zonal toroidal functions of the second kind, g,(u), are defined as 
“u) = SR 
qn(u) 0 (cosh u+sinh u cosh ey 
Also it may easily be shown that q,() can be expressed in terms of p,(u) in 
the form 


ae an(w) = Andale) fsa te tayp 


where A, is an arbitrary constant and is found to be unity in the present case. 
Two cases are now to be considered separately. 


(1) For large values of u 
Using the approximate expression, equation (36), derived for p,(u) in the 
preceding section, equation (37) yields 
(38) (u) eames 3:5 ae . (2n—1) e” })u . een du : 
; dn 22-4... (2n—2) © 4. sinh ue” 
The integral in equation (38) can be evaluated by expressing sinh u in terms 
of its exponential equivalent; therefore, 


3-5... (2n—1) w-px { du 
3S = ae wea. eas 2n—1)u 7 ums 
a WH) = Oa. (On—2)° Sup (Pe) 
3°5 ae 2 —!] n- u ia —2nu 2a, 
= ce a » J e™(1—e ™)" du. 


For purposes of the analysis, let us expand (1 —e~™")~ by use of the binomial 
theorem, and integrate term by term. We have then an approximate ex- 


pression for g,(u) given by: 


Ce ae: or oe ee { 

ees a ae eee 
Ca eo ee eee \ 

“sh... 20 9 “eee ee 


For large values of u (« > 2.0), the result with one term is quite satisfactory ; 


that is 
qo(u) = me *™, 


qi(u) =e = ’ 


q2(u) = " e 


and 


1-3-5... (2n—1) —(ntaru 
“=——-== ¢ : 


(41) W(u) = * TG... On 
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Fic. 3. Zonal toroidal functions of second kind. 


(tt) Series Expression 
A simple but good approximate expression for g,(u) with small values of u 
is somewhat difficult to achieve. However, the author has derived a series 
expression for g,(“), which is valid for both large and small values of wu. 
By replacing cosh uw and sinh uw with their exponential forms, equation (22) 
yields 








r 
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a d0 ’ 
(42) da(u) = J {4e"(1+cosh 6) +3e “(1 —cosh 6)}"*# 


5s penne 0 
= 5 — . 2 
0 (cosh? 40—e¢ * sinh” 49)"*# 


= es | cosh” "* 30{1+(n+})e™ tanh® 30+ .. .} do 
0 








» (1+sinh? 30)"*" 3 





ee ” cosh 36 dé 2n+1 _2% foc 19 sinh’ 36 de 
0 


nae 


If we introduce 
sinh 3 6 = x, 


then, 
4 cosh 3 6 d@ = dx. 


Equation (42) can be reduced to the following expression: 


ee a pe) "x ds \ 
Qa(u) = WS, aaete (FEE) [at f- 
The substitution x = tan @ leads to 
Ss —(n+})u (2n—1) (2n— Pl: . i —u 2n+1 = 
(43) qn(u) = me 2n(2n—2) . er 2 2n+2 


tu (2n+1)(2n43)-1-3 | 
Fe (Qn+2)(Qn+4)-2-4 - 


Equation (43) is the general expression for q,(u), which is valid for all values 
of 1. 
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THE LYMAN BANDS OF MOLECULAR HYDROGEN! 


G. HERZBERG AND L. L. Howe 


ABSTRACT 


The Lyman bands of Hz have been investigated under high resolution with 
a view to improving the rotational and vibrational constants of He in its ground 
state. Precise B, and AG values have been obtained for all vibrational levels of 
the ground state. One or two of the highest rotational levels of the last vibrational 
level (v = 14) lie above the dissociation limit. Both the Bi’ and AG” curves have 
a point of inflection at about v’ = 3. This makes it difficult to represent the 
whole course of each of these curves by a single formula and therefore makes 
the resulting equilibrium constants somewhat uncertain. This uncertainty is 
not very great for the rotational constants for which we find 


B, = 60.864 a, = 3.0764 cm™ 


but is considerable for the vibrational constants we and w,x,. for which three-, 
four-, five-, and six-term formulae give results diverging by +1cm™~. The 
rotational and vibrational constants for the upper state (B'!Zi) of the Lyman 
bands are also determined. An appreciable correction to the position of the 
upper state is found. 


1, INTRODUCTION 


The most precise values at present available for the rotational and vibra- 
tional constants of the hydrogen molecule in its ground electronic state are 
those derived from the quadrupole spectrum in the photographic infrared 
(Herzberg 1950a) and from Raman data (Stoicheff 1957). However, these 
constants represent precisely only the lowest vibrational and rotational levels, 
up to v = 3 and J = 5. For the higher vibrational and rotational levels one 
is dependent on the emission spectrum of Hg, in the vacuum ultraviolet. The 
frequency dispersion in the vacuum ultraviolet is low, and therefore, in the 
early work of Hyman (1930) and Jeppesen (1933), vibrational and rotational 
differences were obtained only with rather modest wave number accuracy. As 
a result there were differences of up to 4cm = between the predicted and 
observed wave numbers of the quadrupole infrared bands. 

Recent improvements in vacuum ultraviolet spectroscopy have made it 
possible to measure lines in this region with much greater accuracy. Because 
of the importance of the Hz molecule in general and of a precise knowledge 
of all its energy levels in particular, it was considered to be of interest to 
remeasure a sufficient part of the far ultraviolet spectrum in order to obtain 
the constants of all vibrational levels of Hz in the ground state. The behavior 
of the levels near the dissociation limit promised to be of particular interest 
(compare the early work of Beutler 1934). 

There are two well-known emission band systems of He in the vacuum 
ultraviolet; the Lyman bands (B—X) representing a 'Z{-'Z} transition and 
the Werner bands (C—X) representing a 'II,—'Z} transition. Since the former 


1Manuscript received February 23, 1959. 
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lie at longer wavelengths and therefore allow one to obtain a higher wave 
number accuracy they were chosen for the present work. Moreover, the higher 
vibrational levels of the ground state are more favored in this transition than 
in the Werner bands on account of the Franck—Condon principle. 


2. EXPERIMENTAL 

The Lyman bands of H: can be produced in any discharge through hydrogen. 
In the present work most spectra were obtained by means of an electrodeless 
discharge tube operated by a high-frequency oscillator (a Raytheon micro- 
therm unit). In order to bring out lines of high rotational quantum number 
a few spectra were taken with a high-current tube with electrodes of con- 
ventional type. The electrodeless tube was also used in some experiments 
with liquid nitrogen cooling in order to reduce the Doppler width of the lines 
and to simplify the spectrum by the exclusion of lines with high J. 

The spectra were taken in the fourth and fifth orders of our 3-meter vacuum 
spectrograph, that is with a dispersion of 0.66 or 0.5 A/mm. The orders were 
separated by means of a foregrating as described by Douglas and Herzberg 
(1957). Standards were supplied by a Fe hollow cathode discharge tube. The 
light from this source was sent through the Hz discharge tube onto the slit 
of the spectrograph. The Fe lines were photographed in second order. No 
attempt to eliminate the error introduced by the method of overlapping orders 
was made, that is the absolute accuracy of the wavelengths may only be 
+0.01 A corresponding to +0.5 cm™' but the relative accuracy (and therefore 
the accuracy of rotational and vibrational differences) was +0.05 cm~' for 
unblended lines. Lines that were more than 0.5 cm~! apart were resolved, that 
is the attained resolving power was well above 100,000. The measurements 
were made with a photoelectric comparator built in this laboratory by Dr. 
B. P. Stoicheff according to the design of Tomkins and Fred (1951). 

In Fig. 1 a section of the long wavelength end of the spectrum is reproduced 
showing some of the bands with the highest v’’ values. 


3. ANALYSIS OF THE SPECTRUM 

The Lyman bands represent a typical many-line spectrum (see Fig. 1); 
individual bands or branches are not at all obvious. However, with the bands 
analyzed by Hyman and Jeppesen, it was relatively easy to proceed and 
find other bands. Some difficulty was encountered in assigning the lines of 
the bands of highest v’”’ since, because of the high D’’ values, the form of the 
branches is rather unusual and since in addition the branches are very short. 
However, the use of combination relations has in each case made a reliable 
analysis possible. 

In Table I are listed the wave numbers of the lines of all the bands that 
have been analyzed. Lines used more than once are indicated by asterisks. 
There are a considerable number of lines which remain unidentified. These 
are listed in Table II for the region 1685-1380 A. It is very difficult to see 
what their identification may be unless they are ‘extra lines’ of some of 
the weaker bands with perturbations. 
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0-3 band 
vo = 78420.99 cm™! 


J P(J) RJ) 
0 78459 . 84 
1 78319.93 436.16 
2 157 .67 349.94 
3 77935. 62 202.32 
4 655.82 77994.53 
5 320.81 728.34 
6 
8 
9 
10 
ll 
0-8 band 
vo = 63372.37 em! 
J P(J) RJ) 
0 
1 63300 .32* 63416.49 
2 195.99 388.63 
3 060.76 327.55 
4 62896. 43 235 .35* 
5 704.70 112.42 
6 
1-7 band 
vo = 67233.92 em~! 
J PJ) RJ) 
0 67270.80 
1 67155. 85 266.16 
2 037.36* 220.13 
3 66879 .80 133.60 
4 685.20 007.89 
5 455.94 66844. 81 
6 194.78 646. 66 
7 65904.64* 416.10 
8 588.72 156.21 
9 250.44 65870.02 
10 64890. 37 
2-4 band 
vo = 77552.82 em"! 
J P(J) RV) 
0 77588 .08 
1 77457 .48 563.14* 
2 302.70 477.65 
3 090.45 333.02* 
4 76822 .57 131.76 
5 502.03 76874.63 
6 131.92 565.26 
7 75715. 46 
5 256.63 


*Blended line. 
tDoubtful. 
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TABLE I 


Vacuum wave numbers of the lines of the Lyman bands of H2 























0-4 band 0-6 band 0-7 band 
vo = 74952.93cm7! vo = 71711.42 em™! vo = 68697.69 cm™! vo = 65915.5: cem™! 
RJ) P(J) RV) P(J) RJ) PJ) R(J) 
74992 .03 71750 .37 68736 .57 65953 . 59 * 
74857.76 973.88 71621.79 737.96*  68613.71 729.85  65837.40 953.59* 
706.60 898.97 481.84 674.29 485.32 677.80 720.83 913.13 
501.38 768.14 293.77 560.40 314.14 §80.90* 567.13 833.90 
243.96 582.66 058.70 397.40 101.51* 440.50 378.21 716.92 
73936.79 344.51 70779 .47 187.06 67850.78 258.44 156.27 563.88 
582.91 055.89 459.80 70931.58 563.71* 036.99*  64903.86 
185.21 73719.74 099.09 633.66 243.13 67777.74 623.79 
72747 .19 336.30 69704.27 296. 16 
272.22 72916.92 271.89 69922.21 
71764.28 454.32 
225.96 
1-3 band 1-4 band 1-5 band 1-6 band 
vo = 79739.40cm™! vo = 76271.39 cm! vo = 73029.83 cm=! vo = 70016.10 em™! 
PJ) RV) P(J) R(J) PJ) RV) P(J) R(J) 
79776.47 76308. 15 73066. 80 70052.99 
79638.50 749.13 76175.98 286.36 72940.14. 050.49 69932 .01 042.48 
474.13 657.09 022.92 205.68 798.37 801.77 69984.63 
248.19 §502.17¢ 75813.85 067.78 606.27 72860.21 626.72 880.53 
78962. 84 285.79 550.97* 75873 .52 365.68 408.76 731.16 
620.79 236.56 625.36 079.10 150.47 539.65 
227.84 68854. 47 306.48 
74466.05 74976.15 71380.09 524.11 035.63 
1-8 band 1-9 band 1-10 band 2-3 band 
vo = 64690.7s cm™! vo = 62397.82 em™! vo = 60371.56 cm! vo = 81020.83 cm™ 
P(J) R(J) P(J) RJ) P(J) RJ) P(J) R(J) 
64727 .89 62434. 14 60408. 40 
64618.70 729.02 62332.28 442.67  60312.92 423.30  80919.91 81025.36 
512.46 695.33 239.04 422.54* 233.84 416.57 80928. 63 
373.30 627.16 119.34 373.22 135.27 389.14 524.79 
203.20 525.85  61975.07 298.41* 019.41 342.11 
004.28 393.25 808.77 197.94  59888.56 277.45 
63779 .54 231.52 622.98 745.54 197.57 
531.82 043.40 420.89 61930.33 593.42 105.10 
266.71¢ 63831. 89 435.51 
2-5 band 2-6 band 2-7 band 2-8 band 
vo = 74311.26cm™! vo = 71297.53cm™! vo = 68515.35 cm™! vo = 65972.2; cm™! 
P(J) RV) PJ) R(J) P(J) RiJ) PJ) RJ) 
74346 .49 71332. 85 66007 .47 
74221.53 327.12 71213.56 319.06  68437.24* 68542.70  65900.21 005.74 
078.10" 253.02 081.62 256.49 317.02 491.10* 792.24 65967.17 
73882 .73 125.56 70903 32 146.17 156.38 649.85 892.71 
637.34 73946.19 680.63 70989.51 67956 .99 474.93 783. 86 
344.59 717.16 416.14 788.76 721.58 269.92 642.58 
007.56 440.99 112.61 546.20 452.85 037 .68 471.11 
72629 38 120.51 69773.44 268.10 64781.24* 272.37 
504.08 
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TABLE I (Continued) 


Vacuum wave numbers of the lines of the Lyman bands of He 


2-9 band 
vo = 63679.25 cm™! 


P(J) RW) 


! 
| 
| 
| 
| 


63714. 45 
719.20 
693 . 66 
638.74 
555.90 
447.10 
314.92 
161.32 


63613.78 
518.82 
395.94 
247.00* 
074.50 

62881. 22 
670.13 
445.08 


CONaukhonre oo | w 


3-5 ban 
vo = 75558.03 cm™! 


P(J) R(J) 


| 
| 
| 
| 


75591. 80 
569.40 
491.12 
358.19 
172.13 

74935.06 
650.26 
320.77 


75468. 26 
323. 
125 

74875 
577 
233. 

73847 


0 
1 
2 
3 
4 
5 
6 
7 
8 
9 


3-10 band 
vo = 62899.76cm™! 


P(J) RV) 


2-11 band 
vo = 59916.33 cm™! 


3-3 band 
vo = 82267.60cem™! 


2-10 band 
vo = 61652.99 cm=! 


PJ) RV) P(J) RJ) P(J) RV) 
61688. 28* 
699. 82 
688. 28* 
654.66 
600.08 
526.82 
437.12 


334.15 


59951. 42 
970.98 
975.28 
965.15 
942.13 
908.54 


59865. 51 82267.90 
800.39 
722.29 
633.41 


536.04 


61594.39 
513.62 
411.92 
291.23 
154.20* 
003.63 

60842.79 
675 


81767 .27 


3-7 band 
vo = 69762.12 em 


vo = 67218.93 cm™! 


P(J) 


vo = 72544.30 em™! 


PJ) RV) RJ) 


69795.95 
785.14 
730.10 
631.86 
491. 87 

68954. 18 312.22 
678.28" 095.71 
371.92 68844.78 
036.99* 

67679. 26* 


P(J) RW) 
67253 .20 
248.08 
205.10 


72578 .05 
561.40 
494.58 
378.75 


67146. 88 
037 .36* 
66892.20 125.29 
713.08 009.50 
502.46 66860.43 
263.90 
65994 . 66 


69683. 97 
562.35 
398.74 
195.10 


72460.04 
326.93 
145.71 

71918.73 
648.41* 
338.49 


70991. 45 466.83 


3-11 band 
vo = 61163.10em™! 


PJ) 


4-3 band 
vo = 83480.5s em7! 


3-12 band 
vo = 59748.12em™! 


RJ) P(J) RJ) P(J) RJ) 





62933. 49 
942.35 
926.59* 
887.35 
826.05 
744.82 
645.91 
533.92 
411.49 
281.07 


62841.18 
758. 88 
654.25 
529.42 
386.83 
229.45 
060.96 

61884.73 
705.63 


4-5 band 
vo = 76771.0; cm™ 
P(J) RJ) 
76803.19 
76681. 11 778.28 
534.70 
334.15 
080.61 
75777 17 
427.21 
034.22 


*Blended line. 
tDoubtful. 


61197.18 
213.41* 
213.41* 
197.70 
168.04 
126.71 
076.57 
021.17 

60959. 67 


59781. 65 
807.34 


21° 83379. 66 
210.94 


82976. 29 


61112 59706. 64 83476. 86 
045.70 
60964.73 
871.46 
768.68 
659.62 
548.21 
438.78 
336.50 
247.30 
179.44 


59395. 


4-8 band 
ve = 68431.96cem~! 


4-7 band 
vo = 70975.10em™! 


RJ) 

71007 .65 

70994 . 22 
935.30 
831.94 
685.94 
498.94 


4-6 band 
vo = 73757.2sem™! 


P(J) RJ) 


P(J) RV) P(J) 
68464 40 
457.12 
410.50 
325. 46 
203.68 
047.32 
67856. 51 


73789. 64 
770.36 
699.66 
578.70 
409.18 
193.34 

72931 .78* 


68359. 94 
249.26 
101.51* 

67918.29 
702.71 
457.42 
185.98 


70897 .06 
773.97 
607.85 
400.39 
154.27 

69872. 46 


73673. 16 
538.46 
354.68 
123.83 

72848. 81 
532.36 
178.21 

71789. 89 


vo = 78799.59 em™! 


P(J) R(J) 


78833 . 48 
805.41 
715. 86 
566.00* 


78548. 18 
332.90 
061.06* 

77734. 66 
357.78 

76933. 83 
465.72 

75958. 74 


092.80 
77773.77 
406.50 


76527 53 


3-9 band 

vo = 64926.02 em=! 
PJ) RJ) 
64959. 81 
961.59 
932.02 
871.35 
781.24* 
665.09 
521.97* 


357.10 


64860. 47 
764.08 
638.30 
485.16 
307.08 
106.99 

63888 .28 
654.33 


4-4 band 
vo = 80012.5; em™! 
PJ) RV) 
80045. 25 
014.55 
79921. 12 


79917 37 
759. 88 
542.05 
266.09 

78934. 86 


279.57 


120.22 78576.07 


4-9 band 
vo = 66139.00em™! 
Pid) RJ) 

66170. 87 
170.87* 
137.05 
071.51 

65975. 81* 
851.81 
702.38 
530.06* 


66073. 55 
65975. 81* 
847.37* 
690.29 
507.14 
300.85 
074.88 
































640 CANADIAN JOURNAL OF PHYSICS. VOL. 37, 1959 
TABLE I (Continued) 
Vacuum wave numbers of the lines of the Lyman bands of He 
4-10 band 4-11 band 4-12 band 4-13 band 5-3 band 
vo = 64112.74cm™ vo = 62376.03 cm™! vo = 60961.10cm™! vo = 59911.92 cem™! vo = 84660.60cm™! 
J P(J) RJ) P(J) RJ) P(J) RJ) P(J) RW) P(J) RJ) 
0 64145.09 62408. 42 60993. 38 59944.07 84670.99 
1 64054. 17 151.38 62325 .33 422.54* 60919.21 61016.44 59880 .73 977.98 84559.75 653. 23 
2 63970. 54 131.82 257.41 418.62 869.32 030.59* 852.37 60013.80 
3 863.34 087 .52 173.79 397 .92 812.92 037.02* 828.97 053.14 152.70 
1 734.57 019.98 076.69 362.10 752.48 037.02* 813.24 098.64 
5 586.88* 63931. 54 61968.80 313.50 691.13 036.08 809.42 154.15 
6 423.41 853.58 252.26 632.79 823.11 224.73 
7 734. 86 582.38 862.13 317.92 
8 617.42 540.24 935.72 
5-4 band 5-5 band 5-6 band 7 band 5-8 band 
vo = 81192.59 cm™ vo = 77951.03 cm™! vo = 74937.30cem7! vo = 72155.12 em=! vo = 69611.9s cm™! 
J P(J) RJ) P(J) R(J) P(J) RJ) P(J) RJ) P(J) R(J) 
0 77982 .35 74968. 39 72186. 59* 69643. 23 
l $1097.32 8119090 77861 .32 954.87 74853 .21 946.78 72077.00 170.56 69539 .93 633.51* 
2 80916.90 713.92 870.15 717.30 872.46 71952.81 108.03 427.99 583.20 
3 718.38 510.46 726 .27* 530.99 746.85 784.12 71999.93 277.61 493.42 
4 254.00 528.28 296.54 571.53 573.02 091.03 365.95 
5 76945.25 277.39 016.69 348 . 82 322.36 68870.71* 202.79 
6 589.33 73694.40 081.80 034.97 619.65 006.09* 
7 189.52 76630.97 73773 .24 70714.27 341.54 68781.18 
8 75749. 36t 039.55 529.30 
5-9 band 5-10 band 5-11 band 5-12 band 5-13 band 
vo = 67319.02 cm“! vo = 65292.76 cm™! vo = 63556.10 em™! vo = 62141.12 em™ vo = 61091.92 cm™! 
J P(J) R(J) P(J) RJ) P(J) RV) P(J) R(J) P(J) RW) 
0 67350. 12 65323 .84 63586. 88* 62172.33 61122.77 
| 67253.20* 346.92 65234. 12 327.62 63505. 27 598.81 62099. 21 192.77 61060.84 154.20* 
2 149.28 304.50 435.88* 591.59 048.11 203.32 030.59 186.46 
3 023.72 039.65 255.39 350.13 565.91 61989. 24 205.05 005.29 221.16 
4 64907.34 182.15 249.41 524.14 925.15 200.09 260.94 
5 66675. 14 754.84 087.05 136.84 468.90 859.06 191.50 60977 .49 
6 462.93 585.60 64972.95 015.93 794.98 182.48 982.68 
7 230.45 403.09 842.90 62890. 55 738.10 178.75 
8 209.52 765.40 693.55 183.83 
9 669.63 
5-14 band 6-5 band 6-6 band 6-7 band 6-8 band 
vo = 60469.93 cm™! vo = 79098.6; cm™! vo = 76084.85 cm! vo = 73302.70 em! vo = 70759.56 cm™! 
J P(J) RV) P(J) RV) PJ) RW) P(J) Ri) P(J) RV) 
0 60501.01 79128.74 73332 .70 70790 .97 
1 60452 .52 546.10 79008 .94 098.94 76090.53 73224.40 314.48 70687 .30 777.94 
2 451.32 606.30 78860 .32 75863 . 87 099.16 248.56 723.84 
3 468.70 684.49 654.62 78862.49 675.49 75882.87 72928 .23 136.00 421.91 629.65 
4 511,40 786.31 394.08 713.74 72978.54 
5 081.56 401.88 458.51 779.40 
6 77720.17 165.64 
7 314.05 71838. 59 262.59 


*Blended line. 
tDoubtful. 
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TABLE I (Continued) 
Vacuum wave numbers of the lines of the Lyman bands of Hz 






6-9 band 6-10 band 6-11 band 6-12 band 6-13 band 
vo = 68466.60 cm™! vo = 66440.34cm™! vo = 64703.6s cm™ vo = 63288.79 cm™! vo = 62239.52 cm™! 

























J P(J) RJ) P(J) R(J) P(J) RV) P(J) RV) P(J) RV) 
0 68496. 63 64733 .51 63318.53 62269.34 
1 68401.00*  491.10*  66381.52 66471.78 64652. 80 742.85 63247.00* 336.80  62208.32 298.41* 
2 300.93 450.40 445.04 582.50 732.00 194.42 343.94 177.69 327.12 
3 167.78 375.75* 494.17 702.04 133.30 341.23 149.44 357.31 
4 003. 68 268.71 389.99 654.94 065.83 330.77 126.68 391.78 
5 67811.37 131.75 211.43 272.98 593.35 62995. 38 315.60 113.66* 433.03 
6 593.94 146.52 521.97* 926.59*  300.32* 113. 66* 
7 354.73 67777.74 014 86 438.81 862.35 286.89 

098.01 810.88 










778.93 


com 






7-5 band 7-6 band 7-7 band 7-8 band 






vo = 61617.56 cm™! vo = 80213.97 em"! vo = 77200.24em™ vo = 74418.0¢6 em7! vo = 71874.92 em™! 

J PW) RJ) P(J) R(J) P(J) RJ) P(J) RJ) P(J) RJ) 

0 61647.08 80243 .15 77229 .00 74446 92 71903. 82 
1 61600.08* 690.06 80124. 38 211.00 77116. 12 202.71 74339. 87 426.52 71802.90 889. 53* 
2 597.46 747.14 79974.71 118.49 76977 .79 121.67 213.36 357.25 688 . 63 832.55 
3 612.74 820.58 766.66 79967.05 786.93* 76987 .25 040.01 240.38 533.71 733.79 
4 651.98 916.98 502.17* 545.94 801.15 73822.81* 078.10* 340.41 

5 78986.04 257.31 565.94 562.68* 73871. 87 111.32 

6 75924. 49 264.79 625.73*  70849.67 

7 550.97* 72931.78* 342.62 559.06 

















7-9 band 7-10 band 7-11 band 7-12 band 7-13 band 
vo = 69581.96 cm™! vo = 67555.70em7! vo = 65819.04 cm~! vo = 64404.06 em™! vo = 63354.8s em™! 

















J P(J) RJ) P(J) R(J) PJ) RJ) P(J) R(J) PJ) RW) 

0 69610. 88 65847 .37* 64432.76 63383 . 62 
1 69516. 45 603.12 67497.09 67583.86 65768. 21 854.82 64362.17 448.74 63323 . 80 410.40 
2 415.19 559.05 409.97 553. 86 696.65 840.58 308.65 452.60 291.98 435.88 
3 279.75 480.09 295.73 496. 33° 606.13 806.44 245.31 445.61 261.40 461.71 
4 112.38 367.99 * 412.13 754.19 174.56 430.10 235.35* 490.86 
5 68915. 81 224.92 66995. 50 304.65 686.54 099.76 408.88 218.25 527.21 
6 815.55 212.94 576.18 
7 448.08 620.51 031.67 230.22 





7-14 band 8-6 band 8-7 band 8-8 band 8-10 band 






vo = 62732.92 cm™! vo = 78284.2¢cm™! vo = 75502.0s em"! vo = 72958.94 cm™! vo = 68639.72 em™! 
J P(J) RJ) P(J) RJ) P(J) RJ) P(J) R(J) P(J) RJ) 
0 62761.61 78312.24 75529. 89 72986 .90 68667 . 61 
1 62715.60 802.21 78200.19* 284.08 75423.78 507.66 72886. 86 970.83 68580.90* 665.07 
2 711.61 855.68 061.06*  200.19* 296.58 435.83 771.70 911.22 492.92 632.49 
3 77868. 42 062.66 121.42 315.65 615.07 809.32 376.84 571.36 
4 624.36 77872.94 149.88 419.07 667.97 235.33 
5 332.88" 632.48 186.59* 070.81 364.03* 
6 76996 .70 339.49 71921 .60 67887 .55 
7 679 .26* 









*Blended line. 
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TABLE I (Continued) 
Vacuum wave numbers of the lines of the Lyman bands of H2 
8-11 band 8-12 band 8-13 band 8-14 band 9-6 band 
vo = 66903.06 em™! vo = 65488.0s cm=! vo = 64438.90 cm™! vo = 63816.94 cm™! vo = 79337.02 cem™! 
J P(J) R(J) P(J) RV) P(J) RJ) P(J) R(J) P(J) R(J) 
0 66930 .81 65515.74 63844 .69 79363 .91 
1 66852 .12 936.15 65446.04 530.06* 63799. 52 883.41 79252 .98 333.53 
2 779.76 919.32 391.95*  530.06* 794.82 934.29 112.64 
3 687.34 881.62 520.85 64536 .04 806.06 64000.36 78917.53* 104.37 
4 577.25 501.87 64313 .92 839.36 087 .52* 670.65 
5 452.56 480.42 374.26 
6 317.71 
9-7 band 9-8 band 9-9 band 9-10 band 9-11 band 
vo = 76554.84 em™! vo = 74011.70cm™! vo = 71718.74 cm™! vo = 69692.4s cm™! vo = 67955.82 cm™! 
J P(J) RJ) P(J) RJ) P(J) P(J) RJ) P(J) R(J) 
0 76581 .90 74038.72 71745. 41 69719.22 67982 .52 
1 76476. 56 557.06 7393951 020.09 71653. 13 733.79* 69633.51* 714.38 67905.02* 985.57 
2 347.98 482.09 822.81* 73957.12 549.85 683.59 545.04 678.42 831.51 965.40 
3 170.78 357 .36 664.28 850. 86* 410.50 596.98 426.20 613.04 736.80 923.43 
4 75947 .33 464.97 703.19 237.10 623.43 861.76 
5 679.70 228.14 516.70 032.76 494.37 782.87 
6 372.19 70800.69 353. 46 
7 028.10 544.37 
s 264.72 
9-12 band 9-13 band 9-14 band 10-6 band 10-7 band 
vo = 66540.8, cm™! vo = 65491.66 cm™! vo = 64869.70 cm™! vo = 80359.24em™! vo = 77577.06em™! 
J P(J) RJ) P(J) RJ) P(J) RJ) P(J) RJ) P(J) RJ) 
0 66567 .30 65518. 26 64896 . 36 80385. 82 77603 .33 
1 66498.85 579.36*  65460.43 541.00 64852. 32 932.68 80275.22 354.02 77498.79 577.62 
2 443.33 426.57 560.61 846.53 980.40 134.26 369.68 500.98 
3 375.85 562.01 391.95* 578.59 855.47 65042.00 79938 .33 191.28 376.72 
4 299.16 359.97 598.26 885.46 123.69 76966 .02 195.79 
5 216.68 505.08 334.93 623.48 956.20 244.87 699.02 76952.00 
6 132.48 382.24 
7 051.88 
10-8 band 10-9 band 10-10 band 10-11 band 10-12 band 
vo = 75033.92 em™! vo = 72740.96cm7! vo = 70714.70cm™! vo = 68978.04cm™! vo = 67563.06 em™! 
J P(J) RJ) PJ) RV) P(J) RJ) PJ) RJ) PJ) RJ) 
0 75060.17 72767 .18 69004.20 67589.05 
1 74961. 66 040.63 72675.43 754.28 68927 .19 006.09*  67521.09 599.97 
2 845.08 74976.15 571.52 702.98 853.47 68984.48 465.02 596.38 
3 684.95 870.30 431.02 615.85 70446.94 757.35 942.80 396.45 582.06 
4 712.74 256.48 642.55 318.34 
5 052.22 131.85 513.22° 236.08 


*Blended line. 
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TABLE I (Concluded) 


Vacuum wave numbers of the lines of the Lyman bands of He 




















10-13 band 





11-9 band 
vo = 73733.03 em! 

















304.99 


10-14 band 11-7 band 11-8 band 

vo = 66513.88 em™! vo = 65891.92 em™! vo = 78569.13 cm! vo = 76025.99 em=! 
J P(J) RV) P(J) RV) P(J) RV) P(J) RV) 
0 66539 . 83 65918.00 78594. 16 
1 66482 .61 561.53  65874.59  953.59*%  78491.62  566.00*  75954.00 
2 448.18  579.36* 868.14 999.76 360.59 485.51 835.72 75960.64 
3 412.50 597.96 875.93 66061.41 179.80 354.14 
4 379.10 904.64* 77950 . 57 173.54 
5 354.29 676.42 77946.27* 224.84 
6 328.29 359.99 672.42 
7 76992. 10 








11-10 band 
vo = 71706.77 cm“ 























11-14 band 
vo = 66883.99 cm! 


11-13 band 
vo = 67505.95 cm™! 


11-12 band 
vo = 68555.13 em! 








RiJ) 





P(J) 
73758.00 
73667.30 742.56 
562.68* 687.26 
419.26 593.59 
240.66 
029.47 302.18 










12-7 band 
vo = 79531.50cem7! 




























































tDoubtful. 


J P(J) P(J) RiJ) P(J) Ri) P(J) RJ) P(J) R(J) 
0 67530. 83 66908. 84 79556. 56 
1 71648.41* 723.32 68513.11* 68588 .38 67474.78 549.95 66866. 56 941.79 79453 . 57 529.13 
2 557.16 682.14 580.90* 439.16 563.71* 858.96 984.04 323.09 
3 435.36 609.64 384.92 559.27 400.98 575.25 864.34 67038.34 142.70 293.46 
4 285.12 302.64 363.65 889.08 111.80 
5 109.15 213.35 925.38 
12-8 band 12-9 band 12-10 band 12-11 band 12-12 band 
vo = 76988.3¢ cm™! vo = 74695.40 cm™! vo = 72669.14cm™! vo = 70932.4s cem=! vo = 69517.50cm™ 
J P(J) RJ) P(J) RV) PJ) R(J) P(J) RJ) P(J) R(J) 
0 77013.30 74720.15 72694.19 70957. 44* 
1 76916.27 76992.10*  74629.90 72610 .67 686.33 70881.73 957.44* 69476.13 69551.23 
2 797.97 929.86 525.61 656.59 519.68 651.33 806.52 
3 636.11 786 .93* 382.02 532.59 398.31 548.95 708.66 859.33 498.92 
4 437.90 209.34 254.49 
5 164.10 73968. 65 048.40 430.27 
12-13 band 12-14 band 13-7 band 13-8 band 13-9 band 
vo = 68468.32 cm™! vo = 67846.36em™! vo = 80464.9, cm™! vo = 77921.80cm™! vo = 75628.84cm7! 
J P(J) R(J) P(J) R(J) P(J) R(J) P(J) RJ) P(J) RJ) 
0 77946. 27* 
1 68437.24* 68513.22* 67905.02* 80457 .32 77849. 84 920.13 75563 .23 
2 401.00* 847.21 75572.55 
3 364.03* 514.16*  67827.33 977.97 80070.90 563.58* 726.27° 310.55 
4 332.87 354.95 563.58" 
5 270.01 104.78 351.62 
6 76817 .47¢ 
13-10 band 13-11 band 13-12 band 13-14 band 
vo = 73602.5s cem™! vo = 71865.92 em! vo = 70450.94cem™! vo = 68779.80 cm™! 
J P(J) RV) P(J) RV) P(J) R(J) P(J) RV) 
0 73625 .73* 71889. 53* 68803. 24 
1 73544 .04 614.32 71815.19 885.56 68762. 59 832.96 
2 451.36 568.30 737.96* 855.33 753.33 870.71° 
3 326.35 489.62 636.78 800.10 70276. 40 755.57 918.83 
4 171.22 379.65 513.43 775.45 
5 72989 . 19 242.86 371.19 093.39 
*Blended line. 
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HERZBERG AND HOWE: LYMAN BANDS 









Most of our assignments were checked by combination differences. The 
agreement of combination differences formed from unblended lines is within 
+0.1 cm. Table III illustrates this agreement for the level v’’ = 11. Similar 
tables were set up for all other vibrational levels of the upper and lower state 
but are not reproduced here in order to save space. 











TABLE III 
AsF’(J) v valees derived from bands with v’ 









429. i 





341. 





69 

















151.03 
3-11 (151 .48*) 248. 68* 341 .95* 429.02 508 . 42 578.50 
4-11 151.01 248 .75* 341.93 429.12 508.52 (578.64) 
5-11 151 .00* 248 .68 (342.18) 429.07 508.21 (578.35) 
6-11 151.01 248 .68 342.01 429 .06 508 . 42 578.49 
7-11 150.72* 248 .69 
8-11 151.05 248.81 342.07 429 .06* 
9-11 151.01 248.77 341.97 429 06 508 . 30 
10-11 150.73 248 .74* 341.93 (429.58) 
12-11 150 .92* 248. 78* 429 .06* 
13-11 78 341.90 428.91 




















95 05 





73 





150. 94 508 . 37 578.50 


Mean 





341. 429 


*F rom blended ane. 
( ) Not included in the mean. 












(a) The Lower State, X 'X} 

In Fig. 2 the values of {A2F’’(J)}/{4(J+4)} are plotted against (J+ 4)? for 
all vibrational levels of the ground electronic state. As is well known, according 
to the relation 






att. | 3 ame eeu .) 1)? 1) 
the intercepts of the lines in Fig. 2 give the B, values, the slopes the D, values, 
and the curvatures the H, values. For the actual determination of these 
quantities the plots were made for each v” separately after addition of 


a(J+3)?—b(J+3)! 


with suitable values of a and 6 to make the resulting line straight and of 
sufficiently small slope that both the remaining slope and the intercept could 
be determined accurately. The value of 6 is then 3H, while the slope of the 
straight line added toa is {2D,—(17/2)H,} and the intercept is {B,—(3/2)D, 
+(27/16)H,}. 

A check on the B, and D, values is obtained in the process of determining 
the vibrational differences AG’’(v+4). This was done by forming the differences 
of corresponding lines of bands with the same upper state and differing by 
Av’’ = 1 in the lower state according to the formula 


(2) Rye(I)—Rerg(J) = Pe(J)—Perai(J) = 46" (v" +4) 


— (By — Bers) I(T +1) + (De — Derry) (I +1)? 
— (He Her) J(I+1P + 
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Sei 4(J+4) 


40 








0 10 20 30 40 £50 60 70 (J+4)" 


Fic. 2. Values of A?F’’(J)/4(J+4) for all observed vibrational levels of the ground state 
of Ho. 


If one adds to the numerical value of the left-hand side of this equation the 
numerical value of the quantity 


(Bor — Byrn) II FY) = (Dor — Dor) PI HIP + Mr — Ber) (I +18 


with the previously determined B,, D,, and H, values and plots against 
J(J+1) one should obtain a horizontal line if the B,, D,, and H, values were 
correct. The slope and curvature of the line actually obtained gives small 
corrections to these. constants. Actually for hydrogen the terms in D, and H, 
are relatively much larger than for most other molecules and it is difficult to 
obtain them with sufficient accuracy from the plots of equation (2). In most 
cases, therefore, the D, and H, were taken from the A:F’’ values and only the 
B, values were slightly adjusted. However, for v’’ = 13 and 14 the number of 
AF” values available is insufficient and the constants were determined entirely 
from the AG plots. The curvature of these last AG plots is so large that it 
cannot be satisfactorily represented by (2). Therefore the B,3, Diz, Bys, Dis 
have much reduced accuracy. 
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In Table IV all the BY’, D’’, and H7’ values thus obtained are collected. 
The constants for v = 0, 1, 2 are those of Stoicheff (1957), except that the 
values for v = 1 have been recalculated from his data using an approximate 
H, value from Fig. 6 below. The level v = 3 is the only one for which our 
measurements can be compared with Stoicheff’s values. Stoicheff found 
B; = 50.6282 and D; = 0.0412;cm™ while we find B; = 50.623. and D; 
= 0.04065. The agreement is quite satisfactory. In Table IV the average of 
these two sets is listed for v = 3. For comparison with the new values the 
earlier B, data of Jeppesen are listed in column 3. Experimental D, values 
had previously been determined only for v = 0, 2, 3 while H, had only been 
determined for v = 0 (Stoicheff 1957). In spite of the lower wave number 
accuracy the present data yield D, and H, of an accuracy comparable to 
the Raman work since they extend to higher J values. 


TABLE IV 


Rotational constants of the ground state of H2 





B,”, cm™ 
ee ee eee D.”, H.", 

v’ This paper Jeppesen cm7 X 10-5 cm 

0 59 .3392¢ 59.373 0.045942 5.29 

1 56 .3685° 56.391 0.0437 5° 

2 53 .47542 53.508 0.0427;4 

3 50.62574 50.676 0.0409,4 3.75 

4 47.801; 47 .850 0.0397, 3.50 

5 44.958, 44.954 0.0385, 3.4 

6 42.096, 42.075 0.03775 3.33 

t 39.1455 39.116 0.03677 3.2; 

8 36. 088, 36.068 0.0365; 3.20 

9 32.833» 32.694 0.03635 3.05 
10 29.311. 0.03665 2.50 
11 25.399, 0.0378, 1.8; 
12 20.926; 0.04135 (1.853) ¢ 
13 15.582s° 0.0492; 3 


14 8.770,° 0.0846; - 








*Stoicheff's values. 

*Obtained from Stoicheff's data using an interpolated H”. 
‘Obtained from AG” plots. 

4A verage of Stoicheff’s and our values. 

*Assumed value. 


The AG” values obtained from the intercepts of the plots of equation (2) 
are listed in Table V. The first three values are from Stoicheff. The values 
of Jeppesen are given in column 3. The fourth column gives 4,G’’(v+1) and 
the fifth column Go(v) obtained by simply adding up the values in the first 
column. 

In Figs. 3a, 4, 5a, 6, 7a, 8 are plotted the values of B,, AB,, D,, H,, 
AG(v+4), and A’?G(v+1). The B, and AG, curves have very similar behavior 
showing an initial slight positive curvature and toward the end a rapid con- 
vergence to zero. The D, curve is of a quite different type reaching very 
high values at the last v’’ values. H, seems to change sign near v = 12. 

The B, and AG(v+3) values are much more accurate than can be repre- 
sented by plots like Figs. 3a and 7a. Even after subtraction of linear terms 
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TABLE V 
Vibrational differences and vibrational levels in the ground state of He 


AG" (v+3) L2G” Go(v) 








v+} This paper Jeppesen This paper This paper v 
0 0 

} 4161. 13,* 4160.96 
235. 162 4161.14 1 

1} 3925 .97;* 3922.18 
230.735 8087.11 2 

2} 3695 .240* 3693 .82 
227 . 230 11782.35 3 

33 3468.01 3467.82 
226.45 15250 .36 4 

4} 3241.56 3239.51 
227 .83 18491 .92 5 

53 3013.73 3013.96 
231.54; 21505. 65 6 

63 2782.18 2782.01 
239.04; 24287 .83 7 

73 2543.14 2543 . 36 
250.18 26830 .97 8 

8} 2292 .96 2292.55 
266.70 29123 .93 9 

9} 2026 . 26 2026 . 28 
289 .60 31150.19 10 

103 1736.66 1737 .53 
321.68 32886. 85 1] 

113 1414.98 1416.92 
365. 80 34301 . 83 12 

123 1049.18 
427 .22 35351 .01 13 
133 621.96 


35972 .97 14 


*Stoicheff's values. 


as in Figs. 3b and 7b a very smooth curve remains which in each case shows 
a point of inflection at about v’’ = 3. The rise of the first few points above the 
straight line representing the next few points is quite striking (see the initial 
drop of the AB and A°G curves) and apparently must be related to a pecu- 
liarity of the potential energy curve. A similar point of inflection has been 
observed for the AG values of the ground state of the O2 molecule (Curry 
and Herzberg 1934). While Hz and O:2 are the only cases of homonuclear 
molecules for which this anomaly has been observed in the ground state, it 
is quite common for polar molecules, for example HF, HCl, and others. It 
may be due to the contribution of ionic states to the ground state, since for 
ionic states a AG curve with positive curvature throughout is expected. 

On account of the initial inflection it is difficult to represent the whole 
B‘'; and AG” curves by polynomial expressions except by using a very large 
number of terms. This is shown by Figs. 3c and 7c where the deviations from 
the five-term formulae 


(3) By = 60.818.—2.9621,(v+4) +0.00458(v+ 4)? 
+0.004548(v+4)*—0.0005105(v+4)4, 

(4) AG’ (v+3) = 4278.70.— 235.4967 (v+ 3) +0.5886,(v+4)? 
+0.28036(v+4)*—0.037903(v+4)4 
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8; - [60.50 -27.68 (v+ $)] 





Fic. 3. By curve for the X !Z} ground state of Ho. 
(a) Plotted directly. 
(6) Plotted after subtraction of a linear term. 
(c) Residuals of equation (3). 
(d) Residuals of equation (9). 


12 1a v" 


Fic. 4. AB, curve for the X 'Z} ground state of Ho. 
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° 2 4 6 8 10 


Fic. 5. D, curve for the X '2} ground state of Hs 
(a) Plotted directly. 
(b) Residuals of equation (15). 


r 1 


aeeemestelieeinenesce ee eS 


a oe 8 10 12 14 v" 


Fic. 6. Hy, curve for the X 'Zf ground state of H2. 
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AG"(v+$) 


AG, - [4265.0 -227.3 (v+#)] 











12 14 v" 


AG(v+4) curve for the X 'Z} ground state of H. 
(a) Plotted directly. 

(b) Plotted after subtraction of a linear term. 
(c) Residuals of equation (4). 

(d) Residuals of equation (14). 


are shown. These formulae were so chosen that the observed B, and AG(v+ 4) 
for v = 0, 3, 6, 9, 12 are exactly represented. It is seen that there are large 
systematic deviations which are much larger than the possible error of the 
individual B, and AG(v+}) values and which cannot be reduced by a different 
choice of coefficients in the equations. It is clear from Figs. 3c and 7c that the 
first term in (3) and (4) does not give the correct B, and w, respectively. Rather, 
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al seiainiaaaaalh a 
8 10 12 14 ing 


Fic. 8. A’G(v+1) curve for the X '!2} ground state of Hz. 
g 


in order to obtain these constants the curves in Figs. 3¢ and 7c have to be 


extrapolated to —}. Such an extrapolation is inherently beset with difficulties 
when the B, and AG values do not follow a simple formula and, therefore, 
the B, and w, values derived in the following discussion have appreciably 
larger limits of error than the individual B, and AG values. 

Stoicheff (1957) obtained from his Bo, B,, Bs, Bz; values, by the method 
of least squares, the three-term formula 


(5) BY = 60.840;—3.0177,(v+4) +0.02855(0+4)? 


which, when the slightly revised B,; and B; values of Table IV are used, changes 
to 


(6) BY = 60.842; —3.0243.6(v+ 4) +0.0302,(v+3)?. 

If the same B, values are represented by a four-term formula one obtains 
(7) BY, = 60.864; —3.08103(v+ 4) +0.0644,(v+ 3)? —0.00570(v+ 3)°, 
while a five-term formula for the first five B, values is 


(8) BY = 60.868;—3.0958;(v+4) +0.0788;(v+4)?—0.0110;(v+4)* 
+0.0006;(v+ 4)‘. 


Finally a least squares four-term formula was developed from the first eight 
B, values (By... Bz): 


(9) BY, = 60.8635—3.07633(v+ 4) +0.0601;(v+ 4)? — 0.0048, (v+ 4)°. 
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As shown in Fig. 3d, the residuals for the first five B, values are less than 
0.002 cm—. This together with the agreement with B, from (7) and (8) suggests 
that a fairly satisfactory extrapolation to —} has been obtained and that 
we can put B, = 60.864, a, = 3.0764, y. = 0.06017, 6, = 0.0048, cm—'. These 
constants are significantly different from Stoicheff’s values which were based 
on a three-term formula. It is believed that the new B, value represents an 
improvement. 

With regard to the vibrational levels the situation is not as satisfactory. 
Using the first three, four, five, and six AG values one obtains the following 
three-, four-, five-, and six-term formulae (the first already given by Stoicheff) 


(10) AG’(v+4) = 4280.299 — 239.457 o(v+3) +2.17257(0+3)?, 

(11) AG’ (v+4) = 4280.666 —240.47253(v+ 3) +2.90500(v+4)? 
—0.1536;(v+4)3, 

4280.173 —238.81983(v+ 3) + 1.2898, (v+4)? 
+0.44733(v+4)3—0.07512,(v+3)4, 


(12) AG” (v+4) 


ll 


4279.590 — 236.7367 (v+4) — 1.05352(v+4)? 
+ 1.58200(v+ 3)? —0.321792(v+4)*+0.019733(0+4)°. 


il 


(13)  AG’(v+3) 


There is no clear convergence. If one uses the first eight AG values for a four- 
term least square solution one finds 


(14) AG’ (v4) = 4282.041;—243.57640(0+ 4) +4.54174(0+ 4)? 
—0.3966;(v+3)*. 


The residuals from this formula are plotted in Fig. 7d. The deviation of 
AG(14) and AG(64) from the formula is far larger than any conceivable error 
of the observed values. One must assume that the true AG curve is strongly 
curved with respect to eq. (14) (see Fig. 7c) and therefore w,, w,.x,... cannot 
be obtained from this formula. The w, values corresponding to the formulae 
(4), (10), (11), (12), (13), (14) are 4396.55, 4400.39, 4401.39, 4399.80, 4398.30, 
4404.59 cm~' respectively. Both the first and the last of these values are 
clearly not the true w, values since the residuals of (4) and (14) plotted in 
Figs. 7c and 7d indicate definitely a large deviation for v = —}. It is very 
difficult to determine which of the remaining four w, values (if any) is the 
correct one since there is no convergence. For lack of better information it 
seems best to retain Stoicheff’s original value 4400.39 cm~ based on the first 
three AG values since there is no good reason why the initial part of the curve 
should depend on the higher terms necessary to account for the higher AG 
values. One must, however, realize that while the w, value in itself is quite 
accurate, it depends strongly on the method of evaluation (i.e. the number of 
terms used). This large uncertainty of w, and the smaller one of B, is unfor- 
tunate since it is the B,, w,.,... values that can be most directly interpreted 
in terms of internuclear distances and vibrational frequencies and which are 
required if one wants to predict the constants of the isotopic molecules HD, 


D., HT, To, ete. 
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The D, values can be represented by a polynomial more easily than the 
B, and AG values. Figure 5b shows the residuals using the four-term expression 


(15) Di = 0.04657;—0.001616(v” +4) — 1.43 X 10-5(v"” +4)? 
+7.6X10-%(0" +3)*. 


As shown by Figs. 3a and 7a, both the B, and the AG curves are likely to 
go through zero between v = 14 and v = 15 in agreement with the fact that 
no bands with v = 15 were found in spite of considerable search. The breaking- 
off of the vibrational structure in the ground state of He does not lead to a 
very accurate value for the dissociation energy Dj since the last observed 
interval is still fairly large. However, it is consistent with the precise value 
D} = 3611646 cm~ obtained by Beutler (1935) from the onset of the far 
ultraviolet absorption continuum. 

The energy of the v’’ = 14 level above v’’ = 0 (see Table V) is 35973.0 
cm. This is the J = 0 level. The J = 1 level lies 2B—4D—8H = 17.2 cm“ 
above this. Adding the observed A:F values successively to the energies of 
the J = Oand J = 1 levels, one obtains the following values for the rotational 
levels of the last vibrational level 


J=0 35973.0, J=1 35990.2, J= 36022.8, 
J=3 36067.6, J=4 36117.8, J=5 36153.4. 


The last observed rotational level is definitely above the dissociation limit, 
while the second last is possibly above it. The existence of sharp rotational 
levels above the dissociation limit is easily understood since the effective 
potential energy curves for J # 0 have maxima which have to be passed 
before the molecule can dissociate (see e.g. Herzberg 19506). Unfortunately 
it is not readily possible to predict the exact height of these maxima unless 
the potential curve for J = 0 were very precisely known for large r values. 
From the observation that lines with J = 5 are still quite sharp while lines 
with J = 6 are absent, presumably on account of extreme width, one must 
conclude that for barriers of more than 60cm~' the H atoms readily pass 
through. It will be interesting to compare the transmission of D atoms with 
those of H atoms through these barriers. 


(6) The Upper State, B'3} 

The molecular constants of the upper state were determined in substantially 
the same way as for the ground state. The combination differences A.F”’ as 
well as the AG’ values agree very well with the values derived by Richardson 
and Davidson (1929) for the lower v values from an analysis of some visible 
bands of H. which have the B'Z state as lower state. We have preferred 
not to mix our data with theirs but have evaluated Bi, Di, AG’ independently. 
On account of a number of perturbations the rotational constant Hj could 
be determined only for a few levels. The average of these values Hj = 1.1 
X10-° cm~ was assumed to apply to all vibrational levels of the B 'Z state. 

Our final values of BY, and D‘, are given in Table VI and the values of AG’ 
and voot+Go(v) in Table VII. The B, values of Jeppesen (1933) which are 
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partially based on those of Richardson and Davidson (1929) are included in 
Table VI for comparison. Similarly the vo0+G’(v) values given by Dieke 
(1958) and also based on Jeppesen’s data are included in Table VII. The 
systematic almost constant deviations in the latter case (up to 9 cm~') are 
apparently due to lack of adequate standards in Jeppesen’s work, which was 








TABLE VI 


Rotational constants of the upper state (B !Z7) 





B,’, cm 















$$ D,’, 
v This paper Jeppesen cm! 
0 19.453; 19.449 0.0154; 
1 18.469, 18.475 0.0133, 
2 17.6535 17.646 0.01215 
3 16.927; 16.927 0.0111, 
4 16.265; 16.260 0.0102, 
5 15.650; 15.621 0.0094; 
6 15.0666 15.000 0.0085; 
7 14.4935 14.406 0.0076; 
8 14.027: 0.0060, 
9 13.4685 0.0058; 
10 13.165 0.0030, 
1] 12.554, 0.0042, 
12 12.610; 
& ] . 00542 









TABLE VII 


Vibrational constants of the upper state (B 'Z*) 






voot+G’ (v), cm! 
AG’ (v+}4), siseea hg ee Ree Sen 























v+} cm"! This paper Dieke v 
90203 .35 90195 .46 0 

j 1318.41 
91521 .76 91513.80 1 

1} 1281.43 
92803 .19 92795 .28 2 

23 1246.77 
94049 .96 94041 .99 3 

34 1212.98 
95262 .94 95255 .00 + 

4} 1180.02 
96442 .96 96434 .95 5 

5} 1147.58 
97590. 54 97582 .45 6 

6} 1115.36 
98705 .90 98698 .0 7 

74 1084.02 
99789 .92 99781 .2 8 

8} 1052.76 
100842. 68 100834 .0 9 

9} 1022.22 
101864 .90 101858 .0 10 

10} 992.07 
102856 . 97 102852 .3 11 

114 962.37 
103819. 34 103812 .0 12 

124 933.44 








104752.78 104751 .3 13 
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particularly troublesome since he used a grazing incidence instrument. Jeppe- 
sen quotes an error of the relative \ values of 0.05 A corresponding to 3 cm! 
in the region in which his work overlaps ours. Our wave numbers are believed 
to be free from such large systematic errors (see Section 2) and, therefore, the 
values of column 3 are considered to supersede those of Jeppesen.? 

In Figs. 9a, 10, and lla the Bi, Di, and AG’ values are plotted. It is seen 
that the shapes of the curves are very different from those of the ground 
state. There is a pronounced positive curvature for both B, and AG’, much 
greater than the initial positive curvature in the ground state. The B12} 
state is generally considered to be an ionic state and the behavior of the 

, and AG’ curves can be understood on this basis. 





. 0 a Mo 
B, curve for the BZ upper state of the Lyman bands of H2. 


(a) Plotted directly. 
(6) Residuals of equation (18). 


*It must be noted that Jeppesen’s band origins vo are larger by approximately 4(B’’ —B’) 
than Dieke’s and ours since he uses the rotational term F(J) = B(J +4)? instead of BJ(J+1). 
His vo unlike ours is not the wave number of the transition J’ = 0 ~ J” = 0. 
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0.000 





Fic. 10. D, curve for the B'2{ state of Hz. 


OG'(v+}) 


4 © 8 12 14 v 


Fic. 11. AG(v+4) curve for the B '!2f state of H:. 
(a) Plotted directly. 
(6b) Residuals of equation (19). 
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As in the ground state, a large number of terms would be required to repre- 
sent all the B’ and AG’ values within the accuracy of the measurements. The 
following two formulae represent the first four B, and AG’ values exactly 
but the higher values increasingly poorly. 


(16) Bi, = 20.035—1.2312(v+4) +0.14518(v+4)?—0.013450(v+4)* 

(17) AG’(v+3) = 1338.22 —40.690(v+ 3) +2.2475(v+ 3)? —0.24167(v+ 3)° 
while a least squares representation of the first 10 values yields 

(18) Bi, = 19.952—1.0619(v+4) +0.06494(v+ 4)?—0.002628(v+ 3)?, 

(19)  AG’(v+43) = 1336.68 —37.503(v+ 3) +0.69516(v-+ 3)? —0.024640(u+ 3)’. 
There are considerable differences between the coefficients in these two sets 
of formulae. It would appear that the first set yields more significant values 
for the equilibrium constants B,, ay, we, WX, . . . and these are therefore chosen 
as the final values in Table VIII. On the other hand the formulae (18) and 
(19) give a much better representation for the higher v’ values. As is shown 
in Figs. 96 and 116. Nevertheless there are systematic deviations which can be 
removed only by introducing still higher terms. Such a calculation did not 
seem worth while. 


TABLE VIII 


Molecular constants of Hz in its ground state (X1Z}) and first-excited state (BZ) 





X13; Bs 


Bo 
Be 


Ge 


59.3392 

60.864 
3.07635 
0.0601, 


19.453; cm™ 


1.2312 
0.14515 


from equations (8) and (14) 


20.035, | 
J 


ye 
Te 0.74158 1.29253 
va" 0.74116 
D, 0.04657 
Be —0.00161, 
Vvoo 0 

ro 0 

We 4400.39 
WeXe 120.815 21.52, 

Wee 0. 7242 0.87005 } 
WZ —0.0604, 


Bed } from equations (12) and (18) 


90203 .35 

91699. 15 
1358.9, | 
from equations (10) and (17) 


The D4, values (Tables VI and Fig. 10) can be represented by 


D‘, = 0.01625 —0.00216;(v+ 3) +2.28) X 10-4(v+ 4)? 
— 1.185 10-5(v+3)8. 


(20) 


Figures 9b and 10 show that the Bi and Dj, valves for v’ = 8, 10, and 12 
are strongly perturbed. The higher vibrational levels of the B 'Z} state starting 
with v’ = 8 are in the same energy range as the lower vibrational levels of 
C'Il,. Indeed, if one plots the (unperturbed) total term values of the B 'Z, 
and C'll, against J(J+1) one finds intersections between B'Z, (v = 8) and 
Cll, (v = 0) at J ~ 7, between BZ, (v = 10) and C'll, (v = 1) at J ~ 6, 
and between B'!Z, (v = 12) and C'll, (v = 2) at J ~ 4, and no intersections 
for B'S, v = 7, 9, 11, 13. These intersections of the (unperturbed) term 
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curves are clearly the reason for the anomalously high B and low D values of 
the levels with v = 8, 10, 12. Unfortunately no observations are available 
for high J values beyond the perturbations, and therefore it is not possible 
to obtain the true Bi, values for these levels. 

The B, and AG’ curves (Figs. 9b and 110) also seem to indicate perturba- 
tions for v’ = 0 or 1, or both, but no perturbing levels are located in this 
energy region. One must therefore conclude that the apparent irregularities 
are similar in nature to those found for the ground state at the lowest v 
values. 


6. CONCLUSION 


In Table VIII are collected the molecular constants derived from the 
analysis of the Lyman bands for the equilibrium positions of the upper and 
lower states. While these constants suffer from the difficulties mentioned 
earlier, it is believed that they represent improvements over those at present 
in the literature. It must be emphasized again that these constants represent 
well only the lower vibrational and rotational levels. For the higher levels 
the observed B,, AG, and D values should be used. 

In Table VIII two internuclear distances are given, one (r,) based directly 
on the observed B, and the other (r,*) after correction for various minor effects 
described more fully by Stoicheff (1957). With the new r?’ value of Hz the 
agreement of the 7, values of Hz, HD, Dz is no longer quite as good as according 
to Stoicheff. Improvements similar to the present ones for Hz in the con- 
stants of De, Tz, and other isotopes appear desirable in order to find the cause 
of the apparent small discrepancies of the r, values. Such investigations are 
being planned in this laboratory. 

It would be interesting if theoretical work could be carried out which 
would give a more detailed explanation of the behavior of the B, and AG 
curves for low v values. Such work would probably also enable one to derive 
more reliable equilibrium values, B, and w,, and might remove the discrepancies 


in the r, values. 
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ANALYSIS OF TWO INTERFERING 1/2+ LEVELS IN THE 
DISINTEGRATION OF C" BY PROTONS! 


A. J. FERGUSON AND H. E. Gove 


ABSTRACT 


Two interfering resonances occurring at bombarding energies of 1.32 Mev 
and 1.51 Mev in the reactions C“(pn)N" and C(py)N* have been analyzed 
in terms of Wigner’s two-level formula. Resonant energies and reduced widths 
have been found which give good fits to the shapes of the yield curves for both 
processes in the neighborhood of the resonances. The absolute cross sections for 
the C'(pn)N" and N4(np)C* reactions are obtained from the analysis without 
reference to the measured cross sections. The calculated cross sections are found 
to be about 40% lower than the measured ones. 


INTRODUCTION 


The analysis of two interfering resonances of the same spin and parity 
presents some features not found in a single resonance or in the interference 
of several levels of different spins and parities. Such a feature is the property 
discovered by Newton (Newton 1952; Thomas 1955) that the normal compound 
nucleus formation may not be, and probably generally will not be, obtained 
in these circumstances. The first study of interfering levels of the same spin 
and parity appears to be that of Fowler and Thomas (1953). Other studies 
have been made by Krotkov (1955) and Vogt (19585). A case bearing some 
resemblance to the two-level one occurs in the reaction C?(py)N™ (Woodbury, 
Tollestrup, and Day 1954) at a bombarding energy of 1.7 Mev where inter- 
ference between a resonance and a non-resonant background is present. Two 
preliminary reports of the present work have been given (Ferguson et al. 
1955; Ferguson and Gove 1956). A complete revision of the calculation has 
been made since the last of these to remove the rather substantial approxima- 
tions contained in the earlier work. 

Two interfering 1/2+ resonances have been observed by Bartholomew 
et al. (1954, 1955) in the reactions C4(pn)N" and C'*(py)N'. These are parti- 
cularly advantageous for the present study because only a single channel spin 
is involved, so that the cross section is given by a single intensity. One proton 
channel and two neutron channels are possible in the energy interval. One 
of the neutron channels involves D waves and is assumed to be negligible. 
This assumption is supported by the neutron elastic scattering studies of 
Johnson, Petree, and Adair (1951). The gamma-ray widths are some orders of 
magnitude smaller than the particle widths, and their contributions to the 
total widths are neglected. Only the ground state gamma ray has been studied, 
which is assumed to be a pure electric dipole transition. Figure 1 gives the 
yield curves obtained by Bartholomew et al. (1955). 


1Manuscript received February 4, 1959. 
Contribution from the Physics Division, Atomic Energy of Canada Limited, Chalk River, 
Ontario. 
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Fic. 1. Yields of ground state gamma rays and of neutrons from C' bombarded by protons 
as functions of the bombarding energy in the laboratory system. The yields are in arbitrary 
units. The gamma-ray yield was measured at 90°. The neutron yield was measured at 90° 
up to 2.2 Mev and at 0° above this. The resonances concerned in the present study have 
energies 1.32 Mev and 1.51 Mev. 





The two-level formula of Wigner (Wigner 1946; Lane and Thomas 1958; 
Vogt 1958a) has been used. Fitting was done by the least-squares method, 
the computations being done on FERUT, the Ferranti electronic computer 
at the University of Toronto. The Wigner formula can be approximated by 
two resonance-like terms combined in amplitude, with an arbitrary phase 
between them. The early work was done on this basis prior to the decision to 
fit directly in terms of the reduced widths of the resonance theory, and this 


will be discussed briefly. 
THEORY 
It will be convenient to use the comprehensive development recently given 
by Lane and Thomas (1958). We will adhere to their notation throughout. 
If only a single spin, parity, and channel spin contribute, then the reaction 
cross section can be written 


2 
7 


(1) Cc’ = Rk? 


Jr 
Vie 








The formulae appropriate to the present work are given in Chapter IX, 
Section 1, of Lane and Thomas. We write for the scattering matrix 
(2) U = aWaQ 
Q(1+%'"(1—RL’) RB’ ’wja. 


Here Q, %, L’, and w are diagonal matrices whose elements are 
Qeer = beer Exp t (we— ) 
Beer = Seer pe/(Fe+Ge) 
Lie = bce(|Se-—-B.+i P,] 
Weer = 24 bee 
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and where w, is the Coulomb phase shift, 


le 
7 tan *y/I, 
l=1 


@, is the hard-sphere phase shift tan“! F./G,; F, and G, are the regular and 
irregular solutions for channel c in the external region; P, is the penetrability, 
p,/(F2?+G?); S, is the energy shift term, P,(F.Ft+G.G~); and B, is the 
arbitrary boundary value parameter. The elements of w are the Wronskians 
of the incoming and outgoing solutions for the various channels; p, = k,r, 
where k, is the wave number and ,, the radial co-ordinate for channel c. The 
prime indicates differentiation with respect to p. 

R is the Wigner R-matrix given by 

(yr. X Ya) 
"= » E,-E 

where y, is a vector whose elements, y,,, are the reduced width amplitudes 
for the various channels. The £, are the characteristic energies of the nuclear 
states; (y, vx) is the direct product of the vectors y, and yy. 

It will be assumed that the entire yield in the energy range considered 
arises from the two 1/2+ resonances. In this case the division of R, 


R = R°+R’, 


given by Lane and Thomas is not useful, so that we put R° = 0 with some 
simplification of the formulae. The results required are given by the two-level, 
many-channel form. To find U we then have 


(3) (1-RL")"R = 2) (1X14) Am 
(Lane and Thomas 1958, Chapter [X, equation (1.13)). Substitution in 
equation (2) gives for a non-diagonal element of U 
(4) Ucer = 2iexp i(wetwe — be— be) (PePer)” a YreYuc'A mu: 
A, is an element of the level matrix A given by 
A = (e—E-§)"; 


e is a diagonal matrix whose elements are the £). E is the energy, E, times the 
unit matrix. — has elements 


0 
iy = > LecVnc¥ pe 
c 


= (Se—Bet+tPc)Yrc¥ne- 


c 
— is separated into its real and imaginary parts 


g = -A+hir 
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where A is the shift matrix and F the width matrix. A and F have elements 


Ay = — be (S.— Be) neues 
Ra 8 27. P tite 
c 


Since S, varies only slowly with energy, we have assumed the arbitrary para- 
meters B, chosen so that 
S,.—B, = 0.* 


The diagonal elements of f are the ordinary widths of the levels familiar from 
the one-level formula. 
For the present two-level case the elements of A are 


An = (E,— E—£&x)/D, 
Ax» = (E,—E—£én)/D, 
Aye = Ag = £12, D. 


Here D = (E,\—E-—£&n)(E,—E—£22) —&?, is the determinant of the matrix 
(e—E-—£). U.., may now be evaluated by substituting these quantities in 
equation (4) and finally the cross section, ¢,.,, found from equation (1). The 
explicit formula which results is given by Lane and Thomas, Chapter XII, 
equation (4.4). 


DETAILS OF THE COMPUTATION 

The two levels involved in the present study occur at proton-bombarding 
energies of 1.31 Mev and 1.50 Mev. The lower resonance has a total width 
of about 40 kev, and the proton and neutron widths have comparable magni- 
tudes. The upper one has a total width of approximately .5 Mev which is 
composed almost entirely of proton width. Figure 1 shows the yields of the 
(pyo) and the (pm) processes from the work of Bartholomew et al. (1955). The 
gamma-ray curve suggests that most of the gamma yield throughout the 
range results from the broad resonance at 1.5 Mev and this yield curve gives 
the best indications available of the width of this resonance. The fit for the 
gamma rays was made over the region 0.91 Mev to 1.61 Mev. A narrow reso- 
nance with spin and parity 1/2— occurs in the range at 1.165 Mev. The points 
in a 60 kev long region at this resonance were omitted from the fit. The fit to 
the (pm) yield was initially made over the same range as for the gamma rays 
but it was soon found that the neutron yield below 1.28 Mev predicted by the 
theory was many times smaller than the observed yield here, suggesting a 
contribution from another process. Thus in subsequent work the fit to the 
neutron yield was confined to the range 1.28 Mev to 1.61 Mev. 


*This approximation need not have been made. At the time the calculations were performed, 
the present development was not available to the authors and the one used required this 
choice. If B ischosen to obtain this result at 1.1 Mev, then the shift term for the broad resonance 
will be 50 kev at .85 Mev and 1.50 Mev, the ends of the regions considered. The shift term 
for the narrow resonance is 1 kev at these points and is negligible where this resonance contri- 
butes strongly. 
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The method of fitting was the least-squares one. However, since the compu- 
tation of the derivatives required for the usual analytical procedure appeared 
somewhat complicated, a trial and error method was adopted wherein the 
fit was computed for a succession of values of one parameter, the value for 
the best fit was selected, and the procedure repeated for another parameter. 
About 10 such cycles of adjustment were used, each cycle representing one 
trial and error adjustment of every parameter. Following this, one final 
adjustment was made using a set of normal equations derived in a way given 
later in the discussion of errors. If we let Vx, be the measured yield and Y¢ 
the calculated yield, and o the standard error of the measurement, then the 
quantity minimized is 

1 9 
(5) A= D3. mr— Yer). 
Tr Oo, 
The summation is over all the measured points, including both the neutron 
and gamma-ray yields. The two yield curves were thus fitted simultaneously. 

The formula for the yield was simplified by factoring out the proton and 
neutron penetrabilities from the (pm) yield and the proton penetrability from 
the (py) yield and incorporating these in the weight factors. For computing 
the penetrability, a radius, a, of 4.774 10-" cm was used, obtained from the 
formula a = 1.4(A!“+1)X10-" cm (Gove 1958). The energies, measured 
vields, weights, and proton penetrabilities were placed in the computer in 
tabular form. The absolute yield of the reactions was not well known so that 
an adjustable normalization factor was placed in the neutron yield. Since the 
gamma-ray yield is homogeneous in the two gamma-ray widths such a factor 
is intrinsically contained in these and would be redundant if put in explicitly. 
Nine parameters were adjusted in all. These were the two resonant energies; 
six reduced widths representing the proton, neutron, and gamma channels 
for the two levels; and the normalization factor for the neutron yield. 

The variances and covariances of the adjusted quantities are given, as is 
well known, by the elements of the inverse of the matrix of the normal equa- 
tions. The normal equations are not obtained directly in the method we have 
used, but can be derived in the following way. Let us denote the nine para- 
meters by y,, 7 = 1, 2,...,9. A can be expanded as a Taylor series in terms 
of the y’s: 


aA 1 a°A 
A=A ~ = x bybyat ... 
ot d ay. by+5 » eae ay by 675+ 
The normal equations are the set given by 
9A 
doy, 


Stopping the series at the second-order term this gives the set 


a°A dA 
> =e y= : 
, O07; OY; 


(6) PS oe eo a 





“ST (BD mt 8 


= if 
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The matrix of the normal equations which we require is thus that whose rsth 
element is 

oA 

OVO. 





Once the minimum of A, which it may be noted corresponds to the vanishing 
of the right-hand sides of equations (6), has been found, these elements can 
be formed by varying the parameters first one at a time, giving 0?A/dy,?, and 
then two at a time, giving 0°A/dy,d7,;. By using both positive and negative 
increments of the parameters, both the first and second derivatives can be 
found if the true minimum has not quite been reached. All the coefficients of 
the normal equations (6) are then known and a further adjustment can be 
made by solving them. This was done at the end of the trial and error fitting 
and produced an appreciable improvement in the fit, demonstrating the 
accuracy of the normal equations. Inverting the normal matrix gave directly 
the standard errors of the parameters. 
RESULTS 

The fitted yields are shown for the (pm) case in Fig. 2 and for the (py) 
case in Fig. 3. The abscissae of these graphs are the center of mass energies, 
which are 14/15 times the proton energies in the laboratory system. In both 
these figures the yield of the 1.165-Mev resonance has been indicated. The 


regions of energy which were omitted from the fit are indicated by giving the 
theoretical yield as a dotted curve. The theoretical (pm) yield vanishes at 


* a 
9 C (pn) , 
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o 1 
« | 
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Fic. 2. Theoretical fit to the yield of the C'(pn)N* reaction. The fit was made only in 
the region covered by the solid line. The abscissae are the energy of the incident proton in 
the center of mass system. The yields are in arbitrary units. 
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1.07 Mev and increases somewhat below this, but does not exceed 1/2% of the 
observed yield down to where the study was started. This vanishing of the 
cross section has been noted by Wigner (1946). 


YIELO GAMMA RAYS) 


RELATIVE 


Ep (Mev) 


Fic. 3. Theoretical fit to the yield of the C4(py)N® reaction. The fit was made only in 
the region covered by the solid line. The abscissae are the energy of the incident proton in 
the center of mass system. The yields are in arbitrary units. 


The parameters for these curves are given in Table I. Eo is the resonant 
energy of the resonance. 7», Yn, and yy are respectively the reduced width 


TABLE I 


Resonant energies, Eo, reduced width amplitudes and partial and 

total widths. The partial and total widths for each resonance are 

calculated at the resonant energy. The energies and widths are for 
the center of mass system 


Narrow resonance Broad resonance 


Eo Mev 1 2302+ .0003 1. 4078+ 0035 
yp (Mev)!/2 0888+ 0035 6084+ 0048 
yn (Mev)"2 1482+ 0020 ~ 0472+ 0013 
+, (ev) 1 4424 044 3.1024 028 
ir hey 6.84.5 400.9463 
I” kev 34 64.9 4 04.2 
r, ev 4.24.7 19.244 
F kev 41.441.1 404.946 .3 


amplitudes for the proton, neutron, and ground state gamma-ray channels. 
r,, l,, and T, are the corresponding partial widths obtained from 


’ i oe 
Pena = 2 I cYa- 


r is the total width given by 


= = Deyn. 
; 
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With the exception of the gamma-ray partial width, these quantities are 
functions of energy by virtue of the penetration factor. The tabulated values 
are calculated at Eo for the level concerned. All of the energies refer to the 
center of mass system. 

The minimum value of A was 509. The expectation value of this quantity 
is (m—r) where m is the number of points and r the number of adjusted para- 
meters. There were 32 points in the (pm) curve, 46 in the (py) curve, and 9 
parameters, giving m — r = 69. A is thus sufficiently different from its expecta- 
tion value to suggest that the fit is not adequate. However, only the errors due 
to counting statistics were included in the calculation, and it is reasonable to 









expect that other errors will be present. 
With the exception of y,, the errors given in Table I have been calculated 


_ (509, -,\"" 
“a 


where N,~! is the diagonal element of the inverse of the normal matrix 
corresponding to the Ath parameter. The errors in y; are derived from the 
uncertainty in measuring the absolute y-ray vield. 






from the formula 








DISCUSSION 






A better physical insight into the formula can be obtained by expressing 
W.. as a sum of resonance-like amplitudes, although it is inconvenient to do 
this for computation. Wigner (1946) has shown how this can be done in 
general by diagonalizing the matrix e—&. In the present simple case W,,- 







has the form 


Wee = 1 Ne /D 


where N,.. is a complex function explicitly involving E to the first degree. 
D is a second-degree polynomial in E which can be resolved into its two fac- 
tors. W,.., can then be resolved into two resonance amplitudes by the partial 
fractions rule. Even in this simple case the expressions are somewhat com- 
plicated and will not be given explicitly. It is found that there is a relative 
phase between the two amplitudes which will vary from reaction to reaction. 
In a case such as the present one where only two dominant channels exist, 
the ratio of these two complex quantities is limited to a relatively small region 
of the complex plane by the requirement that reduced widths be real. 


If 















2 Ah Pe y; V2 
(7) Wyn = 1 \k.—E— (Gl, 5) Re E— (iP2/2)! 





where ¥; and ¥, are the complex amplitudes for the two levels and 
: , 
Wy 


9 





FE,-E- 






ils 
9 


E,;—-E- 





668 CANADIAN JOURNAL OF PHYSICS. VOL. 37, 1959 


are the factors of D, then the ratio Y1/ Y2 is limited to regions of the complex 
plane determined by the ratio ['u/T'22. Wp, refers here specifically to the (pm) 
process and the subscripts 1 and 2 to the lower and upper resonances. In Fig. 4 
is shown the complex plane of Y;/¥2. The region allowed by the requirement 
of real reduced widths is shaded. The crosses representing unrestricted fits 


c(pn) n'* 


A 


\ 
RATIO OF . < 


1.31 Mev: 150 Mev = thes 
AMPLITUDES St ro 


UNRESTRICTED FITS | 


Fic. 4. The complex plane representing the ratio of the amplitudes of the two interfering 
resonances in the C'(pn)N"™ reaction. The allowed region is that for which real reduced 
widths can be found for the Wigner two-level formula. 


were obtained in some earlier work using equation (7) directly. The adjusted 
quantities were | Y,|, | Yo|, arg ¥i/Y2, Ei, £3, and Ty. T2 is not well determined 
by the (pm) yield and a value taken from the (py) yield was used. The energy 
dependence of | Y;|, | Y2|, '¥, and [3 was included in an approximate way. To 
this approximation Y;/Y2 is independent of energy. The restricted fit was 
obtained under the requirements that it lie in the allowed region and be 
compatible with the (py) fit. 

A necessary and sufficient condition, for a definite E, J, and parity, that 
the distribution of products in outgoing channels be independent of the incom- 
ing channel is that the reduced widths of the levels* have the form (Newton 
1952) 

Yui = OB; 
where a, does not depend on the channel, 7, and 6; does not depend on the 
level, u. In the present case this would require 
— Wan; ay 


Y2p Y2n Y2y 


which is clearly not obtained. We thus do not have compound nucleus forma- 
tion in the region where the resonances overlap. Newton has shown further 
that where compound nucleus formation is obtained, then reaction cross 
sections vanish between resonancesf which evidently does not occur here. 


*The levels must overlap. Compound nucleus formation is always obtained at an isolated 
resonance. 

tThis should not be confused with the vanishing cross section that is always obtained, 
either between the rescnances or beyond them, when there are only two channels. 
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Although a knowledge of the absolute yield has not been used in the calcu- 
lation, the results do provide the absolute (pm) and (mp) cross sections. The 
reason for this apparent paradox is that we have been able to determine all 
of the neutron and proton partial widths from the shape of the yield curve 
through the interference which is present. o,, for the narrow level has been 
measured by Johnson and Barschall (1950) and ¢,, for the same level by 
Sanders (1956) and Macklin and Gibbons (1958). A comparison is made 
between these measurements and the values obtained from the theory and is 
given in Table II. 


TABLE II 


Comparison of measured and theoretical cross sections for the (np) 
and (pm) reactions for the level occurring at proton energy 1.230 Mev 
in the center of mass system 


Opn, millibarns 


Onp, Millibarns 


Measured 200 315 
Calculated 131 194 
Ratio, measured : calculated 1.53 1.62 


The agreement between the theory and experiment is poor, and lies outside 
the error of 30% given by Sanders. It is of interest that the one-level formula 
with yi, and y, from Table I gives oy, = 309 millibarns. This shows that the 
interference from the broad resonance is destructive in the region of the 
narrow resonance. An attempt to fit the absolute cross sections has not been 
made, but there is little doubt that the agreement with the shape of the yield 
curve would be materially spoiled in such an attempt. 

The fact that the neutron yield between 0.85 Mev and 1.0 Mev is sub- 
stantially greater than that calculated from the present fit may be significant. 
The yield curve of Macklin and Gibbons shows a typical S-wave threshold, so 
that the yield in this region must be attributed principally to 1/2+ states 
that have not been included in the calculation. If most of the reaction amplitude 
not specifically due to the 1.230-Mev resonance throughout the region fitted 
results from a non-resonant process rather than from the broad 1.408-Mev 
resonance, then y, for the latter may be considerably smaller than quoted. 
With a non-resonant contribution it is also likely that the absolute cross sec- 
tions can be fitted adequately. 

The (np) and (nyo) cross sections for thermal neutron energies computed 
from the present results are given in Table III. The experimental cross sections 
are those quoted by Bartholomew et a/. (1955). The large difference between 


TABLE III 


(np) cross section and (myo) cross section for the 
ground state gamma ray for thermal neutrons 


Calculated Measured 
0.13 barns 1.7 barns 
13 millibarns 


Onp oane 
ony 0.45 millibarns 
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the calculated o,, of Table III and the value of Bartholomew et al. results 
from a slight reduction in the neutron width of the broad resonance and 
from destructive interference from the narrow resonance which causes a six- 
fold reduction of the cross section. This cross section will have a large un- 
certainty due to this cancellation. The calculated o, is the same as the one 
calculated by Bartholomew et al. 

The cross sections for the elastic scattering processes C'(pp)C™ and 
N'4(nmn)N"™ can be predicted from the resonance parameters for the range of 
excitation energies studied provided the potential scattering phase shifts are 
known. No measurement of the C(pp)C"* process has been made. The process 
N'4(nn)N" has been studied by Johnson, Petree, and Adair (1951) and shows 
clearly the 1.23-Mev resonance, but no evidence of the broad 1.41-Mev one. 
The present results predict an extremely small neutron scattering amplitude 
for the broad resonance so that this process should be essentially a one-level 
one involving only the 1.23-Mev resonance. 

The signs of the reduced widths are arbitrary to a large extent and those 
given represent only one of several choices. Any change of sign which keeps 
the signs of the three quantities yipy2p, Yin¥2n, and Yiy72y unchanged also 
keeps the cross sections unchanged. 

The narrow and broad resonances have isotopic spins 1/2 and 3/2 respec- 
tively (Bartholomew et a/. 1955; Bartholomew ef a/. 1956). In the ordinary 
case no interference occurs in the total cross section between states of different 
total angular momentum. Analogously, it might be expected that interference 


will not be present between states of different total isotopic spin, contrary to 
what has been found here. However, the absence of interference in the ordinary 
case requires that the incident beam and the target be unpolarized and that 
polarization be undetected in the final nucleus and the emergent beam. Here 
the incident and emergent beams are completely polarized in isotopic spin 
space so that interference must be expected. 


NOTE ADDED IN PROOF 

A detailed comparison of the (pm) yields of Bartholomew et a/l., Sanders 
(1956), and Macklin and Gibbons (1958) has recently been made. In the 
range 1.19 Mev to 1.51 Mev (center of mass) differences have been found 
between the curves as follows: the energy scale of Macklin and Gibbons is 
18 kev less than that of the other two studies; the peak of the 1.23-Mev 
resonance is about 20% lower in the data of Bartholomew e¢ a/. than in that 
of Macklin and Gibbons; the yield near 1.45 Mev of Bartholomew et al. is 
about 40% higher than that of Macklin and Gibbons. Correspondence between 
the yields to better than 10% is otherwise obtained. An error of 18 kev may 
thus be present in the resonant energies found here. Uncertainties in the 
reduced width amplitudes revealed by these differences are difficult to estimate, 
but are probably less than those arising from other sources such as neglect 
of non-resonant contributions. 

Mr. Syurei Iwao (private communication) has recently extended the 
calculation by including a 1/2+ resonance near an excitation energy of 9 Mev 
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in the N® nucleus. This has successfully accounted for the yield below 1.19 
Mev and has reduced the disagreement between the measured and calculated 
cross sections at the peak of the 1.23-Mev resonance. 
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NOTES 





RADIATION FROM A SMALL LOOP IMMERSED IN A 
SEMI-INFINITE CONDUCTING MEDIUM 


JAMeEs R. Wait* 


It is well known that a circular loop antenna situated in the air and with 
its axis parallel to the ground will radiate a vertically polarized ground wave. 
Furthermore, if the electrical constants of the ground are known, the field 
strength of the ground wave can be readily calculated from convenient curves 
prepared by Norton (1941). He also discusses in detail the effect of raising 
the transmitting and receiving antenna above the ground. 

In this note the effect of burying one or both of the transmitting and 
receiving loops in the ground is considered. For purposes of analysis the 
ground is assumed to be represented by a semi-infinite homogeneous medium 
of conductivity oc, dielectric constant e. The air is assumed to be a semi-infinite 
insulator with dielectric constant ¢9. The permeability of the whole space is 
taken as u, equal to that of free space. At low radio frequencies the insulated 
loop can be represented as a magnetic dipole oriented in the direction of the 
loop axis. The fields vary according to the factor exp (wf). 

Choosing an (x, y, 2) co-ordinate system, with the insulating medium 
defined by z < 0 and the conducting medium defined by z > 0, the dipole 
is located at (0, 0, #) and oriented in the y direction. The electric and magnetic 
fields in the conductor can be represented in terms of a Hertz vector F as 


follows: 

(1) E = —curl F 

(2) luwH = —y?F+erad div F 
where 


2 


y? = lopw— euw?. 
A subscript 0 is added to the quantities in order that they shall pertain to the 
insulating region, so that 
2 2 
Yo = —€quw 
since 
KB = Ho. 
The formal solution of the problem can now be obtained in a straight- 


forward manner by matching cylindrical wave functions at the boundary 
z = 0 (Sommerfeld 1926). The formal result is given by Wait and Campbell 


(1953). 


*Present address: National Bureau of Standards, Boulder, Colorado. 
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0 (wyot+uoy’) (uo+u) 
= e'" /r, P2 = e''*/r» 

and where 

p= (x?+-y2)4 

ry = [p?+(2—h)?I}, ro = [p?+(2+h)*]}, 

Uo = (2+y2)}, u = (A?+7?)3, 
The constant C is proportional to the strength of the dipole and is given by 

C = tuwldA/4m 

where J is the total circulating current in the loop and dA is the area of the 
loop. In a previous paper (Wait and Campbell 1953) the integrals were evalu- 
ated in closed form with the restriction that |yop| « 1 and certain approximate 
results were given for other cases. It is the purpose of the present note to 
show that expressions for the radiation fields can be readily obtained by 


generalizing Sommerfeld’s results for dipoles at the interface. 
For a highly conducting ground, it is permissible to make the substitution 


u = y(1+d2/72)? & y(1+72/72)? & ¥ 


since the saddle point of the integrals is near 4 = yo when |y*| > |y3|. This 
approximation was used by Sommerfeld (1926) originally in the evaluation 
of the integrals for the case when the dipole is at the interface. The integrals 
Q and R can then be written 


Q ~ wn ee f (uotvi/y) ‘AJo(Ap)dd 
70 


ibicube 0 rag 2 a = ‘ 
Ree -2teen f eet at/a) = Gctod debe 


Integrals of this type can be expressed in terms of the error function as shown 
by Niessen (1948) for the case where |yop| >> 1. Thus R and Q are given by 


QO Way ep Fp) 


Re Qy7* i 6 TM 610° F(p) — F(p’)] 
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where 
pb = —v5e/(27?) and p! = —(yop/2)(y*/7}) 
and F(p) is an attenuation factor given by 
F(p) = 1—i(rp)* e~erfc. (ip) 


while for large values of p, F(p) = —1/(2p). 
The fields H, and H, are then given by 


(6) inwH, &% —2Cyi e761" pF (p) sin’ 
and 

(7) ipwH, & 2Cyi eT "np" F(p) sin @ cos > 
where 


@ = cos! y/p. 


‘These expressions, which are valid for |yop| > 1, were given in a previous paper 
(Wait and Campbell 1953) without proof. 

When (z+) tends to zero, the fields are the well-known expressions for 
the far-zone field intensity of a loop antenna situated on the surface. The 
effect of burying the transmitting loop and also the receiving loop is to simply 
modify the magnetic fields for the case of the loops on the surface by a factor 
exp[—y(z+)]. This conclusion is not modified for distances even when the 
restriction |yop| > 1 is violated as long as |yp| > 1 and p> (z+A). 

The additional attenuation G of the far-zone field, due to the immersion of 
the loops in the conducting medium, is then simply given by 


G = (z+h)/6 nepers. = 8.68 (2+h)/5 db 


where 6 is the “skin depth” of the radiation in the conducting medium given 
by 
5 = real part of 1/7 & [2/(cuw)]}. 


It is also of interest to note that an H, component is present but it varies 
as p~ in the radiation field and is not of any practical significance. 
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NOTES 
THE NOSE-ON RADAR CROSS SECTIONS OF FINITE CONES* 


H. Brysk, R. E. Hiatt, V. H. Weston, AnD K. M. SIEGEL 


Keys and Primich (1959) have completed the most thorough determination 
of the radar cross sections of perfectly conducting finite cones known to the 
authors. Certain features of these data are stated to disagree with theoretical 
predictions (Siegel 1958). The disturbing fact that some of the data of Keys 
and Primich (1959) do not fall off like 1/A* at the long wavelength end has 
led us to examine these results critically and to perform some experimental 
measurements of our own near the Rayleigh region. Part 1 explains briefly 
our nose-on Rayleigh predictions for the finite cone. Part 2 discusses the 
resonance region for the cone. Part 3 describes measurements we have made 
since obtaining a preprint of Keys and Primich (1959). Our results indicate 
no need to revise the theoretical discussion of Siegel (1958). At the same time, 
we find good agreement in the experimental data. 


1. RAYLEIGH REGION THEORY 


The reasoning behind the Rayleigh region theory, as given in Siegel (1958), 
is as follows: In a scattering experiment the probe measures results which 
depend upon polarization, wavelength, and properties of the body. Here we 
are discussing nose-on cross sections from a perfectly conducting body of 
revolution; the constitutive constants of the scatterer are thus simply fixed, 
and it is known in this case that the return is independent of polarization. A 
priori we expect the radar cross section in the Rayleigh region to have an 
inverse fourth power wavelength dependence. Since dimensionally the cross 
section is an area, we know this wavelength dependence must be multiplied 
by a factor of the dimension of length to the sixth power. When the wave- 
length is large in respect to the dimensions of the body, the body looks like a 
singularity to the wave. The return must then be insensitive to the fine 
structure of the body, because a geometric singularity on the body (e.g., a 
point, a line, a sharp interface) looks to the wave merely as a singular detail 
superimposed on a singularity. Thus the obvious first choice for a length 
to the sixth power is the volume squared. Accordingly, we write 


(1) o = Ck*(VF)?, 


where o is the nose-on backscattering radar cross section, k = 24%/X where 
d is the wavelength. C; is a pure number, V is the volume of the scatterer, and 
F is a dimensionless correction factor that takes into account any shape 
dependence of the cross section. In view of the expected insensitivity of the 
cross section to fine details of structure, it should be possible to express F 
fairly accurately as a function of a length-to-width ratio only. 


*Work carried out under contract with the Rome Air Development Center of the Air 
Research and Development Command, United States Air Force. 
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Rayleigh himself has solved exactly the scattering of an electromagnetic 
wave incident along one of the axes of a spheroid for the case in which the 
wavelength is large compared with all the dimensions of the spheroid. Thus we 
know C,F* for the spheroid. This expression tends to a constant value for 
very prolate spheroids, and to infinity for very oblate ones (since the disk has 
a finite cross section despite a vanishing volume). A very simple approxi- 
mation, which reproduces Rayleigh’s complicated formulae to within 1%, is 
(Siegel 1958): 

(2) F = 1+(sy)7 €*, 

(3) Ci = 4/n, 

where y is the ratio of the principal axis along the direction of incidence to 
that normal to it. 

The proposed extension of this result to shapes for which an exact solution 
is not available consists of retaining eq. (2), with y redefined as a constant 
times the ratio of the maximum dimensions along the direction of incidence 
and normal to it; the basis for the choice of the constant is that, in the limit 
as the shape is flattened, the disk cross section should emerge. In this way, a 
finite cone of base radius r and height A is assigned a cross section (Siegel 
1958): 

(4) o = (4/n)R VAL + (ny) 1eP, 
where 
V = (1/3)xrh, and y = h/4r. 


Keys and Primich compare cones of total angle 8° and 120° and point out 
that the theoretical difference in cross section on a volume-squared basis is 
28 db, whereas experimentally they obtained a difference of 4 to 5 db. If the 
F factor is taken into account, the theoretical cross sections differ only by 18 
db. There remains a discrepancy of the order of 14 db. This is, however, easily 
explained by the fact that Keys and Primich were using in their comparison 
cross-section values for the 8° cone at kr ~ .45, which is not in the Rayleigh 
region. For a cone with total angle @ and base radius 7 to be in the Rayleigh 
region, the following must hold: 

kr cot (0/2) < 1, 6 < 90° 
Rr 6 > 90° 
or, for an 8° cone, kr < tan 4° = .07. 

In our experimental results, we have yet to hit the Rayleigh region for our 
small angle cones, but we are approaching it. Our large angle cones are already 
there. Our experimental and theoretical results appear to be consistent. 


2. RESONANCE REGION THEORY 


In Siegel (1958), it was shown that, for small wavelengths, the main contri- 
bution to the nose-on cross section of a finite cone is from the ring singularity 
at the base. The physical optics results were computed by using wedge-like 
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local fields in evaluating the contribution of the ring singularity. The minor 
cone resonances (past the first) were predicted to lie at the same locations 
as for the ring. It was felt that the rise on the small wavelength side of the first 
resonance could be approximated by the ring. Keys and Primich have shown 
experimentally that the locations of the resonances (other than the first) 
agree indeed for rings and cones. 

We have pointed out that in the Rayleigh region the backscattering cross 
section for incidence along the axis of a perfectly conducting body of revolu- 
tion depends primarily upon the volume of the body, modified by a factor 
involving the length-to-width ratio. This implies that if the orientation of the 
body is reversed (front to rear) the Rayleigh return should not change sig- 
nificantly. Siegel has predicted that this should remain true as one enters 
the resonance region from the Rayleigh side—specifically that, if a cone in 
this region is turned around, the return should not be much affected. Honda 
(Honda, Silver, and Clapp 1956) has measured this front-to-rear ratio for 
several cones. In the small wavelength region, the reversed cone gives a return 
two orders of magnitude greater than nose-on. (Theoretically, this can be 
understood as due to the flat plate effect of the base in the reversed position.) 
For the smallest cone, on the other hand, which is just about at the first 
resonance (apex angle 15°, kr = 1.5), the rear contribution is larger by only 
a factor of 2+.5. This suggests that the location of the first resonance for the 
cone can be approximated by that of the resonance for the cone’s base, namely 
the disk. 

Comparing disk, sphere, and prolate spheroid (Siegel et a/. 1956) results, 
for incidence along the symmetry axis, it is evident that the location of the 
first resonant point, as a function of kr, where 7 is the maximum radius of 
revolution, decreases with elongation. Hence, the disk provides the upper 
bound for the location of the first resonance for spheroids. Since the disk is a 
limiting case for the cone too, we expect the disk answer also to be the upper 
bound for the cone. 

The first resonant point is the largest value of \ (or the smallest value of 
1/A) for which do/dy\ = 0. Hence, to exhibit the first resonant point, one 
must plot o vs. A (or 1/A), or some constant multiples of these quantities such 
as (o/mr*) vs. kr. On the other hand, if the ordinate (¢/A\*) is chosen (as was 
done in Keys and Primich (1959)), the maxima of the curve obtained are 
zeros of d(¢/\*)/dd, which do not necessarily coincide with zeros of do/dx. 
On a (o/mr’) vs. kr plot, Keys and Primich’s results for the cone show a first 
resonance at kr = 1.5. This compares with kr = 1.1 for the ring (Kouyoum- 
jian 1956; Weston 1957) and kr = 1.6 for the disk (Andrejewski 1952). 

These results will now allow us to predict the location of the resonance 
region for other bodies with greater precision than we previously could, 





3. PRELIMINARY EXPERIMENTAL RESULTS 


Measurements have been made at four different frequencies on the nose-on 
radar cross sections of small cones with apex half angles of 7.5°, 9.6°, 12°, 45°, 


and 60°. The four frequencies (3020, 3800, 9300, and 9900 Mc/s) were used in 
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order to provide more experimental points for a given number of models, and 
also to cross-check results obtained at different frequencies with different 
equipment. Measurements were also made on a cone with a half angle of 4° 
at frequencies of 8465 and 9900 Mc/s. Values of kr for the cones measured 


ranged from about 0.12 to 4.5. 
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Fic. 1. Nose-on cross sections of cones with various half angles, 6/2. O Experimental 
results; Rayleigh theory. 


The measurements were made on an indoor range, in a room lined with 
microwave-absorbing material. Low power, continuous wave oscillators, a 
microwave receiver with a sensitivity of about —86 dbm, and a hybrid tee 
waveguide system comprised the principal components of the equipment. 

All measurements were made at a distance equal to or greater than 2D?/), 
where D is the maximum dimension of the horn aperture or of the cone, 
whichever is greater. Calibration spheres were chosen so that their cross 
sections approximated those of the cones. Measurement procedures were as 
follows: With no target, the return from the background and the polyfoam 
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support cylinder was balanced out with the waveguide tuners. The return was 
then recorded for one or two spheres, for four or five cones, and then for the 
test spheres again as each was placed in position. If the system went out of 
balance, or if the second reading on the test spheres differed significantly from 
the first reading, the data would be discarded. The cones were not rotated, 
but repeated measurements, made after carefully repositioning the models 
on the support cylinder, were consistent. The points that appear on the 
accompanying curves represent the average of at least three, and usually as 
many as six or eight, separate readings. Before averaging, the values generally 
agreed to within +1/2 db. 

The agreement obtained between X-band and S-band data was very 
satisfactory. The agreement between our results and those given in Keys and 
Primich (1959), where the values of apex angle and kr in the two sets of data 
overlap, was most gratifying. 

It is interesting to note that, over the wide range of apex angles covered 
by the data, the location of the first maximum varies very slowly with cone 
angle. 
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Under this heading brief reports of important discoveries in physics may be published. These 
reports should not exceed 600 words and, for any issue, should be submitted not later than six weeks 
previous to the first day of the month of issue. No proof will be sent to the authors. 


Extension of Davydov’s Theory of Molecular Crystals to Vibronic Transitions 


The theory of the electronic states of molecular crystals was developed by Davydov (1948) 
in terms of exciton states. In these calculations (Davydov 1948; Winston 1951) only the 
electronic wave functions of the molecule were used, so that the resulting expression for the 
transition energies holds only for the pure electronic transition ve which does not occur in the 
observed spectrum. Instead, one observes transitions »,"," with vibrational quantum numbers 
v’, v’’ in the upper and lower states, respectively. It is therefore not correct to compare calcu- 
lated Davydov splittings with those observed in the spectrum, unless the effect of molecular 
vibrations on the calculated splitting is included. This is true even for the 0,0-transition, 
which is the main source for observing the Davydov splitting. Simpson (1957) was the first 
to show how vibrational effects can be taken into account in such calculations. One of the two 
limiting cases which he discussed is treated here in more detail. This case applies to Davydov 
splittings which are small compared with the molecular vibrations. 

Here the crystal energies can be calculated as first-order perturbations of the vibronic states 
of the free molecule in the same way presented by Fox and Schnepp (1955) for the pure elec- 
tronic states. In their formulae the wave functions yi, Yj of the ground and jth excited state 
of a molecule at the sth site in the crystal will be replaced by the product of these functions 
with the vibrational functions ¢}, ¢3. This yields the same eq. (6) of Fox and Schnepp (1955) 
for the excitation energy of a crystal state 


(1) E;,9(K) —Eo = (€; —€0) +D; g +M;,.(K). 


The first term on the right side is the excitation energy of the free molecule in its jth state 
including vibrational energies. The second term is independent of the wave vector K and 
only shifts the free molecule transition whereas the third term splits it into a band of states 
characterized by K. D and M are still given by eq. (7) and (8) of Fox and Schnepp (1955) 
except that the vibronic molecular wave functions are substituted for the electronic ones. The 
interaction operator V between any two molecules A and B can be written V = Vee+ Ven+ Van, 
corresponding to interaction between the electrons of A and B, the electrons of A and the 
nuclei of B and vice versa, and between the nuclei of A and B. Only Vee gives a non-zero con- 
tribution to the band-splitting term M, if the Born—Oppenheimer approximation holds. The 
resulting expressions for D and M can be reduced to one- and two-electron integrals through 
the use of the spinless electron densities p{", p{*, p{*}, which refer to the ground state, the jth 
excited state, and the transition 0 — 7 respectively. Thus the following expressions are obtained 


(2) bye= DD art f (oe? ay—oe"a)) 


| Vae(1) +f 2) Veo(1,2) ra a 


® Ae) =- > Somer j ps? (1) p6™ (2) Veo(1,2)dridrs 


x(feeresrao™) ; 


From this the following conclusion can be drawn. The band shift term (2) remains un- 
changed by including vibrations; it is only a more condensed form of eq. (7) of Fox and Schnepp 
(1955). The band-splitting term (3) is the product of two factors. The first is the same term 
as calculated for the pure electronic transition. The second factor gives the influence of the 
vibrational motion of the molecules. It is equal to the square of the vibrational overlap integral 
between the ground and excited state (Simpson 1957). Jt will lead to a different band splitting 
for each vibronic transition. On the other hand, the order of the different components—numbered 
by the subscript g—will be the same for all vibronic transitions. 
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The electronic transition moment for the transition from the crystal ground state to the 
excited states can be calculated in the same way, giving 


@) Ryne = dV TD Be) iax f oParrae”, 


The first factor gives the electronic transition moment, expressed in terms of the moments 
(ry) ;,0 of the free molecules in one unit cell. This is multiplied by the same vibrational overlap 
integral, which occurred in eq. (3). 
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Fox, D. and SCHNEPP, O. 1955. J. Chem. Phys. 23, 767. 

SIMPSON, W. T. and Peterson, D. L. 57. J. Chem. Phys. 26, 588. 
Winston, H. 1951. J. Chem. Phys. 19, 156. 


RECEIVED OCTOBER 27, 1958. W. A. BINGEL 
DEPARTMENT OF PHYSICS, 

DuKE UNIVERSITY, 

DuruHaM, NortTH CAROLINA. 





THE PHYSICAL SOCIETY 


MEMBERSHIP of the Society is open to all who are interested in Physics. 

FELLOWS pay an Entrance Fee of £1 1s. ($3.15) and an Annual Subscrip- 
tion of £2 2s. ($6.00). 

STUDENTS: A candidate for Studentship must be between the ages of 18 
and 26, and pays an Annual Subscription of 5s. ($0.75). 

MEETINGS: Fellows and Students may attend all Meetings of the Society 
including the annual Exhibition of Scientific Instruments and Apparatus. 


PUBLICATIONS include the Proceedings of the Physical Society, published 
monthly, and Reports on Progress in Physics, published annually. Volume XXI, 
1958, is now available (price £3 3s. ($9.45)). Members are entitled to receive 
any of the Publications at a reduced rate. 


Further information can be obtained from: 


THE PHYSICAL SOCIETY 
1, LowTHER GARDENS, PRINCE CoNnsorT ROAD 
Lonpon, S.W.7, ENGLAND 











NOTES TO CONTRIBUTORS 
Canadian Journal of Physics 


MANUSCRIPTS 


General.—Manuscripts, in English or French, should be typewritten, double spaced, on 
paper 83 X11 in. The original and one copy are to be submitted. Tables and captions 
for the figures should be placed at the end of the manuscript. Every sheet of the manuscript 
should be numbered. Style, arrangement, spelling, and abbreviations should conform to the 
usage of recent numbers of this journal. Greek letters or unusual signs should be written plainly 
or explained by marginal notes. Characters to be set in boldface type should be indicated 
by a wavy line below each character. Superscripts and subscripts must be legible and carefull 
placed. Manuscripts and illustrations should be carefully checked before they are suluaitead. 
Authors will be charged for unnecessary deviations from the usual format and for changes 
made in the proof that are considered excessive or unnecessary. 

Abstract.—An abstract of not more than about 200 words, indicating the scope of the work 
and the principal findings, is required, except in Notes. 

References.—References should be listed alphabetically by authors’ names, un- 
numbered, and typed after the text. The form of the citations should be that used in current 
issues of this journal; in references to papers in periodicals, titles should not be given and 
only initial page numbers are required. The names of periodicals should be abbreviated in 
the form given in the most recent List of Periodicals Abstracted by Chemical Abstracts. All 
citations should be checked with the original articles and each one referred to in the text by 
the authors’ names and the year. 

Tables.—Tables should be numbered in roman numerals and each table referred to in the 
text. Titles should always be given but should be brief; column headings should be brief and 
descriptive matter in the tables confined to a minimum. Vertical rules should not be used. 
Numerous small tables should be avoided. 


ILLUSTRATIONS 


General.—All figures (including each figure of the plates) should be numbered con- 
secutively from 1 up, in arabic numerals, and each figure referred to in the text. The author's 
name, title of the paper, and figure number should be written in the lower left corner of the 
sheets on which the illustrations appear. Captions should not be written on the illustrations. 

Line drawings.—Drawings should be carefully made with India ink on white drawing 
paper, blue tracing linen, or co-ordinate paper ruled in blue only; any co-ordinate lines that 
are to appear in the reproduction should be ruled in black ink. Paper ruled in green, yellow, 
or red should not be used. All lines must be of sufficient thickness to reproduce well. Decimal 
points, periods, and stippled dots must be solid black circles large enough to be reduced if 
necessary. Letters and numerals should be neatly made, preferably with a stencil (do NOT use 
typewriting) and be of such size that the smallest lettering will be not less than 1 mm high 
when the figure is reduced to a suitable size. Many drawings are made too large; originals 
should not be more than 2 or 3 times the size of the desired reproduction. Whenever possible 
two or more drawings should be grouped to reduce the number of cuts required. In such groups 
of drawings, or in large drawings, full use of the space available should be made; the ratio of 
height to width should conform to that of a journal page (4) X7} in.), but allowance must 
be made for the captions. The original drawings and one set of clear copies (e.g. small 
photographs) are to be submitted. 

Photographs.—Prints should be made on glossy paper, with strong contrasts. They 
should be trimmed so that essential features only are shown and mounted carefully, with 
rubber cement, on white cardboard, with no spice between those arranged in groups. In mount- 
ing, full use of the space available should be made. Photographs are to be submitted in 
duplicate; if they are to be reproduced in groups one set should be mounted, the duplicate 
set unmounted. 


REPRINTS 


A total of 50 reprints of each paper, without covers, are supplied free. Additional reprints, 
with or without covers, may be purchased at the time of publication. 

Charges for reprints are based on the number of printed pages, which may be calculated 
approximately by multiplying by 0.6 the number of manuscript pages (double-spaced type- 
written sheets, 8$X11 in.) and including the space occupied by illustrations. Prices and 
instructions for ordering reprints are sent out with the galley proof. 








Contents 


H. R. Fickel and R. H. Tomlinson—Thermal neutron absorption cross section 


a a i a ee oe CM Ie SO 
Bengt Kleman and Ulla Uhler—A *S—*S transition in CrH - - - - - 
K. W. Geiger—Absolute standardization of radioactive neutron sources. II. 

The use of the F’® (a,m) Na* reaction - - - - - = = = - 
M. Srirama Rao—Size of irregularities in the E region of the ionosphere - - 


A. G. Fenton, K. G. McCracken, D. C. Rose, and B. G. Wilson—Transient de- 
creases in cosmic ray intensity during the period October 1956 to January 





1958" <0 6. Ss ee Se ee ee ws ee Se ee 
W. J. Bratina—Investigation of deformation processes in Armco iron by means 
of internal friction at megacycle frequencies - - - - - = = 
Archibald W. Smith—The impedance, rectification, and electroluminescence 
of anodic oxide films on aluminum - - - - - + - - «= - 
J. P. Roalsvig, R. N. H. Haslam, and D. J. McKenzie—Photoneutron cross 
section of Ni® - - - - = = 5 2 2© 2£© 2£© 2£ 2 = « 


K. L. Chopra and T. S. Hutchison—Study of phase propagation in supercon- 
ducting aluminum by ultrasonic attenuation measurements - -_ - 


S. C. Loh—On toroidal functions - - - - - = = - -« «= - 
G. Herzberg and L. L. Howe—The Lyman bands of molecular hydrogen - - 
A. J. Ferguson and H. E. Gove—Analysis of two interfering 1/2+ levels in the 


disintegration of C‘ by protons - - - - - - = = - «= - 
Notes: 
James R. Wait—Radiation from a small loop immersed in a semi-infinite 
conducting medium - - - - = = = *© =*#© = -# « -« 
H. Brysk, R. E. Hiatt, V. H. Weston, and K. M. Siegel—The nose-on radar 
cross sections of finite comes - - - - - - - - - =. - 


Letters to the Editor: 
W. A. Bingel—Extension of Davydov’s theory of molecular crystals to 


vibronic transitions = set’ Get Gms fe) Gee ec Fe ener 


Printed in Canada by University of Toronto Press 








531 
537 


550 


557 


569 


579 


591 


607 


614 
619 
636 


660 


672 


675 








